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FOREWORD to the Second Edition 


This document is the second edition of the SIMWEST operating manual. The 
SIMWEST program described in the first edition was capable of modeling 
total wind energy storage systems. This edition also includes a descrip- 
tion of recent enhancements to the program which give it the capability to 
model solar photovoltaic systems. These enhancements were developed under 
NASA contract DEN3-42 "An Expanded System Simulation Model for Solar En- 
ergy Storage." The principal investigator for this contract was Dr. A. W. 
Warren. Co-investigators were Dr. Y. K. Chan and Dr. M. H. Dwarakanath. 
This program was conducted under the sponsorship of the Division of Energy 
Storage Systems, DOE, under the direction of Dr. G. C. Chang, and was 
administered by the NASA-Lewis Research Center Thermal and Mechanical 
Storage Section with Mr. L. H. Gordon and Mr. R. H. Beach as project 
managers. 
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1.0 INTRODUCTION 


Energy storage systems for the utilization of intermittent power sources 
have received increased study over the past few years. The analysis of 
storage requirements for optimal utilization of solar-derived energy sys- 
tems and the total cost of the resulting generator/storage system are often 
evaluated in such studies. The purpose of the SIMWEST (Simulation Model 
for Wind Energy Storage) program described in this document is to provide a 
tool for performing this needed analysis. It is a tool to aid in the 
design of a wind or solar-photovoltaic energy system for a given applica- 
tion and to allow the resulting system to be evaluated and verified through 
simulation. 

SIMWEST consists of a library of system components and a precompiler pro- 
gram which allows these components to be put together in building block 
form. The present library contains components for five types of energy 
storage systems. They are pumped hydro, battery, thermal, flywheel, and 
pneumatic. The SIMWEST program version described in this document is for 
use on the UNI VAC 1100 series of computers. 

The simulation program has proven to be efficient and versatile for per- 
forming parametric studies. It has a unique capability for simulating 
total wind/solar systems containing any one or combination of the above 
types of storage and at the same time has the flexibility and depth re- 
quired to perform thorough and meaningful parameter studies. 

1.1 SIMWEST OVERVIEW 

SIMWEST consists of two basic programs, and a library of generation, stor- 
age, environmental, and load components. The first program, the Model 
Generation Program, is a precompiler which generates computer models (in 
FORTRAN) of complex energy generation/storage systems, from user specifi- 
cations using SIMWEST library components. The second program utilizes the 
resulting computer model to perform cost and power utilization analysis. 
It handles input, output, integration of system dynamics, and iterates to 
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obtain convergence of implicit variables. The combination of these two 
programs provides a powerful tool for analyzing alternate generation and 
storage system designs. ^ 

Figure 1.1-1 shows the general organization of the SIMWEST program. In 
addition to the two programs described above, there is a third which 
performs file maintenance. It is used to incorporate user supplied data 
for new subsystem models. Although the program is shown as a number of 
subprograms, it can be executed as a single batch program by supplying the 
model description cards and the control cards describing the desired anal- 
ysis to be performed and the desired tabular and/or plotted output. , 

The SIMWEST model generation and simulation programs have a number of user 
oriented features which greatly enhance the value of the codes. Some of 
the more prominent features are shown in Table 1.1-1. These features and 
the supplemental components described in 1.2 enable the user to quickly 
build, debug, simulate and interpret alternative system designs. 

1.2 SIMWEST LIBRARY 

The SIMWEST library is listed in Table 1 . 2-1. It is made up of six types of 
components: environmental, generation, load, logical, storage and supple- 
mental. The two character mnemonic names are used to identify components 
in the users model. 

The degree of detail in the component models is based upon two design 
criteria. First, all models should contain sufficient detail to simulate 
all physical characteristics and constraints having significant impact on 
system cost effectiveness. Second, the models should be designed to mini- 
mize computer time and required user specification. It is assumed that a 
SIMWEST simulation might cover a time span of one year. Thus, from a 
computer run time and economic impact point of view a simulation step size 
of between 15 minutes and one hour was established as a design goal. 

a. 
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Figure 1.1-1 SIMWEST Program Organization 



Table 1.1-1 SIMWEST User Oriented Features 


MODEL GENERATION PROGRAM 

• Simplified Component Connections 

• Availability of all Input Parameters for Connection 

• Fortran Insertion Capability Between Components 

• Line Printer Schematic of User’s Model Provided 

• Automated Naming of Parameters and Variables 

• Built-in Diagnostic Capabilities 
SIMULATION PROGRAM 

• Free Field Data Inputs, Including Tables 

• Diagnostics on Data Inputs 

• Default Values Assigned to Unspecified Parameters 

• Optional Levels of Line Printer and Diagnostic Output 

• Multiple, Back-to-back Simulation Capability 

• Printer Plotter Output of Time Histories and Crossplots 
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As a result of the above design criteria, many physical components, such as 
the electrical components, were modeled mainly in terms of power flow and 
steady state response. This level of detail is consistent with a 15 minute 
time step and with the concept that important transients are on the time 
scale of demand curves or weather patterns, i.e., an hour or more, rather 
than on the time scale of electric motor transients of a few seconds. If 
short time transients were to be modeled, additional detail would be re- 
quired in the component models which would greatly increase the user's task 
of specifying the model. Further, the simulation time step would have to 
be reduced and computer runs would be much costlier. 

The environmental components listed in Table 1.2-1 simulate environmental 
conditions. In the present SIMWEST library a user can generate wind speed 
and ambient temperatures, or can use selected inputs from the recorded 
weather and insolation data on the Typical Meterological Year (TMY) tapes 
for one of 26 U.S. locations. These variables are generally used as inputs 
to physical components. 

The generation components consist of wind generation, solar-photovoltaic 
and utility routines. The wind turbine-generation components are fairly 
simple models for computing the power output of a conventional, horizontal 
axis wind machine given basic machine parameters. The solar-photovoltaic 
components are somewhat more sophisticated, especially in the collector 
thermal analysis, and have a number of modeling options which a user may 
employ, e.g., active or passive cooling. 

The storage components encompass such things as motors, generators, trans- 
missions, and flywheels. These components model actual physical hardware 
which might be used in a wind or solar energy system. The selection of the 
particular SIMWEST library set of storage components was based on the 
requirement that it be capable of modeling the five types of energy storage 
systems mentioned previously: thermal, flywheel, battery, pumped hydro 
and pneumatic. 

The load components in the SIMWEST library are used to simulate various 
types of power demand. They also monitor how well the system meets the 
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Table 1.2-1 SIMWEST Library Components 


ENVIRONMENTAL 


WIND WD 
AMBIENT TEMP TP 
TMY WEATHER TAPE ED 

WIND POWER GENERATION 

fURBI NE/GENER ATOR WP 
WIND TURBINE WT 
FIXED RATIO TRANSMISSION GR 
AC GENERATOR GE 

SOLAR POWER GENERATION 

SOLAR ORIENTATION (TRACKING) SO 
FLAT PLATE COLLECTOR FP 
FOCUSING LENS COLLECTOR FO 
PHOTOVOLTAIC ARRAY PV 


BATTERY STORAGE 


INVERTER IV 
RECTIFIER RE 
BATTERY BA 
ADMITTANCE AD 

FLYWHEEL STORAGE 

AC MOTOR MO 
VARIABLE RATIO TRANSMISSION TR 
FLYWHEEL/CLUTCH FL 

HYDRO STORAGE 

HYDRO PUMP PU 
HYDRO TURBINE HT 
HYDRO STORAGE HS 


PNEUMATIC STORAGE 


UTILITY GENERATION 


UTILITY UT 

LOGIC 

POWER DIVIDER PD 

POWER ACCUMULATOR PA 

PRIORITY INTERRUPT PI 

SWITCHES SW.SX 

SY.SZ 


LOAD 

ELECTRICAL LOAD LO 

THERMAL LOAD TL 


COMPRESSOR 

CO 

TURBINE 

TU 

ADIABATIC HEAT EXCHANGER 

HX.HY 

BURNER 

BN 

PNEUMATIC STORAGE 

CS 

THERMAL STORAGE 

STORAGE VESSEL 

TS 

SUPPLEMENTAL 

SATURATION 

SA 

RANDOM NUMBER GENERATOR 

RN 

TEST FUNCTIONS 

AF 

TABLE LOOKUPS 

FU.FV 

TRANSFER FUNCTIONS 

IT.LA.LL.TF 

ARITHMETIC ELEMENTS 

m,MB,MC 

COST MONITOR 

CM 

HISTOGRAM 

HG 

TAPE READ 

TA 

TIME CONVERSION 

TI 
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simulated demand and compute the value of the energy delivered to the load. 
Like the environmental components, these components may be computed from 
actual measurement data or from randomly generated data based on user 
furnished load profiles. 

The library’s logical components are the power dividers, power accumu- 
lators, switches and priority interrupts. Although physical hardware or 
logic devices could be built to serve the function of the logical com- 
ponents, they are not meant to represent any particular existing hardware. 
Instead, they are idealized components that allow the user flexibility in 
modeling a wide variety of system and control logic for operational evalua- 
tion of energy storage systems. In practice, the control function might be 
performed by a control room operator using a predefined control strategy or 
by use of a process computer. 

Finally, the supplemental components include such things as the tape read, 
the histogram and the cost monitor. These components serve to help the 
user run the simulation and analyze its results. , 

1.2.1 Storage Subsystems 

Figures 1.2-1 through 1.2-5 give example configurations of the five types 
of storage subsystems which can be modeled with the present SIMWEST 
library. For illustrative purposes the number of variables shown passed 
between components is limited. A description of the variables being passed 
is given in Table 1.2-2. 

A total energy system will generally be made up of elements from a number 
of different subsystems (see Figure 1.2-6). In addition, the SIMWEST 
program can be used for models which include networks of storage subsystems 
of the same type or a network of wind or solar generators. 

1.2.2 Logic Components 

The capability for modeling complex system control logic is provided by the 
power divider, power accumulator and priority interrupt components. Both 
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Figure 1.2-1 Pneumatic Storage Subsystem 
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Figure 1.2-3 Flywheel Storage 
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Figure 1.2-5 Thermal Storage 


Table 1.2-2 Partial List of Component Inputs and Outputs 

SYMBOLS 


p 

POWER 

RE 

POWER REQUEST 

MP 

MAXIMUM POWER 

RS 

ROTOR SPEED 

T 

TEMPERATURE 

TA 

AMBIENT TEMPERATURE 

M 

MASS FLOW RATE 

H 

RESERVOIR HEIGHT 

LO 

THERMAL LOAD DELIVERED 

WV 

WIND VELOCITY 

GR 

GEAR RATIO 

EF 

EFFICIENCY 

INT 

INTERRUPT FLAG 

PR 

PRESSURE 

PS 

PRIORITY SEQUENCE 

WY 

WEEK OF YEAR 

DW 

DAY OF WEEK 

TD 

TIME OF DAY 

SP 

SURPLUS POWER 
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the divider and accumulator operate on a priority basis. The priority 

interrupt is used by other system components to change the priority setting t 

of the divider and accumulator. ^ } 

The power divider has one input power port and four output power ports (not 
all output ports need be used for a given simulation). The divider also 
has an input request associated with each of its output ports. A power 
request originates with a component which is directly or indirectly con- 
nected to an output port. The user specifies priorities of either 0, 1, 2, 

3, or 4 to be associated with each of the output ports. If the input power 
exceeds that requested of the port with highest priority (priority 1) then 
the excess power goes to the port with the next priority. This process 
continues until either all power is distributed or all requests of non-zero 
priority ports are met. A port with zero (0) priority does not receive 
power. Such ports are included to model backup or switch operated com- 
ponents. In these situations, the connected component would change the 
zero priority setting of the power divider by use of a priority interrupt. 

Two or more ports may be assigned the same priority in which case the user 
may specify weights to be associated with each port. Then if there is not 
enough power available to satisfy all requests of equal priority, the power 
is divided between them in proportion to the user specified weights. 

The power accumulator is similar to the divider except that instead of 
distributing power from a single input port among four output ports, it 
accumulates power from four input ports and sends it out through a single 
output port. The power accumulator accepts power requests from the down- 
stream component and allocates requests to each of its input ports in order 
to service the downstream component. 

An example illustration of the use of power dividers and power accumulators 
is given in Figure 1.2-6. It is seen that power from the turbine/generator 
is distributed with highest priority (priority 1) going to the power accum- 
ulator that services load 1. Since the power accumulator servicing load 1 
has its priority 1 input port connected to the power divider, it will try 
first to satisfy load 1 from the turbine/ generator and then from the 
utility. If the power divider satisfies load 1 and there is power left v 
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over, it will be used to satisfy the request from the battery. Finally, if 
the battery is full or if its charging rate is met, then the excess power 
goes to the flywheel. The battery also has a priority zero connection to 
the utility. If the battery remains in a discharge state for more than a 
specified amount of time, it can change the utility priority (from 0 to 1) 
to receive needed power. 

Also in Figure 1.2-6, we see that load 2 prefers to draw power from the 
flywheel before turning to the battery. This configuration tends to keep 
the flywheel as discharged as possible, using it primarily as a means to 
absorb large influxes of power. 

1.3 SIMWEST OUTPUT 

There are three basic forms of SIMWEST output to facilitate the analysis of 
wind and solar energy storage systems; line printer plots, histograms of 
system variables and time sequenced output of variable values. Each 
SIMWEST library component is associated with a number of output variables. 
Prior to simulating a given system the user may select any of these outputs 
for plotting or tabular output. For example, he may want to plot the 
energy of pneumatic storage as a function of time and/or as a function of 
temperature. If the user wants a time sequenced listing of all variable 
values, he may specify the time step between printouts. The listing of all 
variables has proven to be a useful tool in understanding the performance 
of the system under consideration and a valuable aid in validating the 
system design. 

SIMWEST also provides a special output which computes life cycle and level- 
led energy costs per kwh. This output is produced by the cost monitor 
component and is illustrated in Figure 1.3-1. The levelized energy costs 
are based on energy delivered to the loads during a simulation and fore- 
casted to a full years' system operation. This output permits direct 
comparison of capital and energy costs for alternative system configura- 
tions, enabling a user to perform economic trade studies and system sizing. 
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• 
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Figure 1.3-1 Cost Monitor Output For Fresnel Lens Model 
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1.4 TESTING 

Reference [l] describes two simulation studies which were used to test the 
original SIMWEST program. Reference [6] describes the NASA-Lewis approved 
simulation studies for the expanded SIMWEST program. These studies 
provide an excellent test and illustration of the program's capability to 
model complex wind/solar energy systems. 

Prior to performing the simulation studies and throughout its development 
the SIMWEST program was systematically tested. First components were 
grouped into simple systems and simulations were performed. During these 
simulations system parameters were driven so as to force the individual 
components through every normal program path and to assure that all com- 
ponent outputs assume a wide range of values. The number of components and 
the number of ways they can be connected makes it impossible to exercise 
every combination. However, the subsystem groupings that were used were 
representative of the expected program usage. Sections 8 and 9 describe 
some of the test cases for the wind and solar-photovoltaic generation 
components. 

In terms of computer efficiency, it was found during the testing that the 
program exceeded original expectations. Even on very complex systems, 
such as represented by the NASA-Lewis test case, convergence of logic 
variables was quite rapid. Convergence generally took place in less than 
six iterations per simulation time step. As an example, the year simula- 
tions used in the NASA defined parameter study of reference [1] took less 
than 420 CPU seconds on the CDC 6600. For comparison, the CPU time on the 
UNIVAC 1100/40 is approximately two to three times as great as that on the 
6600, and CPU time on the Cyber i/5 is a factor of two to three times 
smaller than that of the 6600. 

1.5 PROGRAM USAGE 

While the user need not be a SIMWEST expert or software specialist to make 
efficient use of the program, he should thoroughly think through and be 
familiar with the characteristics of the system he wants to simulate. 
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Component models, if not carefully specified, may perform in unexpected 
ways. If the systems logic is not well thought out, the resulting system 
may be significantly out of balance and subsystems may not be fully uti- 
lized. The test case described in reference [6] illustrates the process of 
sizing and logic adjustment to satisfy system performance objectives. 

A number of useful procedures were developed during the simulation 
studies. First it was found that when simulating a complex system, it is 
best to separately develop and test subsystem portions of the model. This 
allows problems or unexpected results to be isolated and understood prior 
to the introduction of the more complex characteristics associated with 
the total system. 

It was found during the simulations that the use of Fortran statements in 
the model definition is very useful for creating special input to system 
components and for defining special outputs to be plotted or statistics to 
be printed. For example, Fortran statements enable the user to generate 
and interpret trade study data by computing component input parameters 
from user specified system parameters. The use of Fortran statements is 
simple and should be encouraged early in SIMWEST applications. 

Computer simulation costs may be minimized by appropriate tradeoffs be- 
tween run time and simulation accuracy. Run time is most directly affected 
by the integration step size, the total simulation length, and the average 
number of iterations through the model at each time step. For long dura- 
tion runs, an hour step size is usually acceptable. Models having smaller 
time constants than the step size may be approximated by implicit steady 
state conditions and solved by iteration through the model. If a model 
requires many iterations for convergence then it may be useful to isolate 
the source of instability in order to modify or simplify that portion of 
the system model. It has been generally found in the simulation studies 
that use of a few seasonal weekly simulations is adequate to predict long 
term performance for system trade studies and design optimization. Based 
"he results of [6], four to six week long simulations are recommended 
this purpose. 
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When making a year simulation run, it is best to break it into twelve 
monthly simulations. Thus, measures of performance such as plots, histo- 
grams and performance statistics are available on a monthly basis. In 
addition to giving better visibility of the system performance, this helps 
limit the job core size. The twelve monthly simulations can be submitted 
as a single run with the results of a given month acting as initial 
conditions for the next month. The user only needs to submit new data 
cards for data which changes from one month to the next. 
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desired, the independent and dependent axis scale ranges can also be speci- 
fied. The independent scale range is specified by the word XRANGE followed 
by the minimum and maximum values for this scale. The dependent scale 
similarly is specified by the word YRANGE. If scale ranges are not speci- 
fied, values will be used that span the given data. 

SI MANUAL SCALES 

SI AUTO SCALES (Default Condition) 

The SI MANUAL SCALES command allows the plotted output requested by the 
DISPLAY commands to be plotted on manual scales specified by the YRANGE and 
XRANGE commands. The SI AUTO SCALES command can be used to return plotting 
to the automatic scaling mode. Auto scales are selected so that they span 
each plotted quantity. The auto scale option is the default used until 
manual scales are requested. The PRINTER PLOTS command is also required to 
obtain plots. 

Example 3.5-1: 

SI MANUAL SCALES, PRINTER PLOTS 
DISPLAY1 

WV2WD, VS, TIME, YRANGE =10,40 

PI PD, VS, TIME, YRANGE = 0,1000 

P2 PD, VS, TIME, YRANGE = 0,1000 

DISPLAY2 

P2 IV, VS, TIME 

RE2BA, VS, TIME 

RE1L0, VS, TIME 

0ISPLAY3 

PI PD, VS, P2 PD, YRANGE = 0,1000, XRANGE = 0,1000 


NG PAGE NOT FILMED 
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TITLE 


The TITLE command allows a title to be placed on all plotted output. Up to 
74 characters may follow the delimiter that follows the TITLE command. The 
TITLE command may be changed before each analysis. Once defined, the title 
remains in effect until a new title is entered. 

Example 3.5-2: 

TITLE = BATTERY TEST MODEL 

3.6 ITERATION AND DIAGNOSTIC CONTROL 

There are three built-in parameters in any SIMWEST model : CYCLES, DLINES 
and RESET. These parameters are specified similar to component parameters 
using the PARAMETER VALUES command. 

CYCLES controls the number of iterations through the model to obtain steady 
state. If CYCLES < 0. then only one pass is made through the model. If 
CYCLES is a positive integer then the maximum number of iterations through 
the model is equal to CYCLES + 1. If cycles is positive, but not an 
integer, then the maximum number of iterations is equal to the smallest 
integer value exceeding cycles. A maximum of 20 iterations are permitted 
per time step. Most of the models tested require no more than six itera- 
tions per time step to attain steady state. A complex model with cascaded 
logic components may require more. 

Each of the model output variables are monitored each pass for convergence. 
If all of the outputs are converged within 3% of their previous values, 
then one final pass is made through the model. Otherwise, all variables 
exceeding 5 % of their previous value are printed out after the last itera- 
tion. 

Since output statistics are only updated the last iteration, some of the 
variables printed indicating nonconvergence are just statistics, and as 
such should be ignored. 


50 


BCS 40180-2 Rev. 



DUNES controls the amount of convergence-related printout to be con- 
trolled as well as the amount of diagnostic printout put out by the library 
component. If DUNES >0 then the total number of diagnostic printouts is 
no greater than DUNES. Figure 3.6 shows a typical section of diagnostic 
printout using DUNES >0. If DUNES <0 then only library component diag- 
nostics are printed with no greater than - DUNES of output. Typically, 
DUNES = 50 is sufficient to catch most simulation errors per run. 


TS STORAGE TEMPERATURE 50, i** OUTSIDE MINIMUM *0,000 AND MAXIMUM 212,000 

-18 STORAGE TEMPERATURE - 5**TJl OUTSIDE MlNIrUR *0,000 -- AND MAXIMUM 212,000 ' 

TINE* 00, SO 

P2 NT nOnCON Vf P6ENCE j— OLD -VALUE ■ Jl,*tl - NE* VALUE* SO.SOO — “ 

P2 CE NONCONVERGENCE, OLD VALUE* JO, *10 N£« VALUE* 20,0*0 

PL 6E NOnCOnvERGENCE, OLD VALUE* 1,275 NE* VALUE* 1,2 11 

NS RESERVOIR VOLUME 77210,004 DROPED fcELO* NINIMUN 00000,000 


-TO-OTORAliE TEMPERATURE 50,**4 OUTSIDE MINIMUM ■ *0,000 — AND MAXIMUM 2 12,000 

TS STORAGE TEMPERATURE 50,*** OUTSIDE MINIMUM *0,000 AND MAMIMUN 212,000 

TS STORAGE TEMPERATURE 5*. *S* 0UTS10I MINIMUM *0,000 AND maximum 212,000 


FIGURE 3.6 TYPICAL DIAGNOSTIC OUTPUT 

RESET controls the initialization value for the random number generators 
if several simulations are run back to back. If RESET >0 (Default) then 
the same random numbers are used for each simulation. If RESET<0 then the 
random numbers at the start of each simulation are obtained from the last 
value at the end of the previous simulation. 

3.7 DEFINE COMMANDS 

DEFINE STATES 
DEFINE RATES 
DEFINE PARAMETERS 
DEFINE VARIABLES 
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These program commands may be used to define the alphanumeric names that 
will be used to refer to states, rates, parameters, and variables. All 
system models formed by the Model Generation program have model-related 
names generated for all states, variables, and parameters in the model. 
State variable derivatives, (Rates), are generated as Rl, R2, ... for all 
models. Rl, R2, ... refer to the rates of the first, second, ... states 
respectively. If it is desired to replace these machine generated names 
with other names, the DEFINE command may be used to substitute any eight 
character names of the analyst's choosing. These names are associated with 
the corresponding numeric quantities located in the labeled commons /CX/, 
/CXDOT/, /CP/, and /CV/. The appropriate location for each quantity is 
printed out along with the quantity name prior to each simulation. Each of 
these commands is followed by phrases containing the location numeric fol- 
lowed by an alphanumeric name with one to eight characters, the first of 
which must be alphabetic. 

Example 3.7: 

DEFINE STATES 

1 = PRESSURE, 2 = STROKE, 5 = VELOCITY, 7 = ANGLE 

DEFINE PARAMETERS 

5 = MASS, 35 = DCT AREA 

DEFINE VARIABLES, 1 = T OUTLET, 2 = LIQ H20 

Note that the program commands, numeric values and alphanumeric names must 
be separated by delimiters which are: [,] , equals [=] , left parenthesis 
[(], right parenthesis [)], or three or more consecutive spaces. 

3.8 EXAMPLE OUTPUT 

Figure 3.8 shows a sample of the output print format generated using PRINT 
CONTROL = 3. This sample is taken from the Wind Turbine and File Read run 
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4.0 JOB CONTROL PROCEDURES 


In this section, we describe job control procedures for running and maintain- 
ing the SIMWEST programs. For the convenience of the user, a number of proce- 
dure files have been set up which simplify the user control cards required. 

In Section 4.1, we describe the control cards for executing the model genera- 
tion and analysis programs. Section 4.2 describes the procedures to maintain 
the programs and update the component library. 


4.1 MODEL GENERATION AND ANALYSIS EXECUTION 

Figure 4.1-1 shows an overview of the program structure to execute a simula- 
tion run. The program FILOAD is only executed when the component I ibrary 
Is updated, and is thus described in the next section. The user input data 
for the model generation program is put on a file called EASYCARDS. A proce- 
dure file called XQTEASY is then used to generate the model Fortran and com- 
pile this model. Similarly, the user input data for the analysis program 
is put on a file called NONSIMCARDS, and a file called XQTANALYSIS maps the 
relocatable elements into absolute file elements, and executes both the simu- 
lation and printer plot programs. 


A job control stream to execute these programs in a batch environment is given 
by: 


6RUN ... 

€DELETE,C EASYCARDS. 
WSG,UP EASYCARDS. 
ODATA, IL EASYCARDS. 


INPUT DATA DECK 
FOR MODEL 


OEND 

0*SG,A XQTEASY. 
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♦SCRTCH9 IS FORTRAN SOURCE CODE OUTPUT 




FIGURE 4.1-1 SIMWEST PROGRAM EXECUTION STRUCTURE 



@ADD,PL XQTEASY. 
@DELETE,C NONSIMCARDS. 
@ASG,UP NONSIMCARDS. 
@DATA,IL NONSIMCARDS. 


INPUT DATA DECK 
FOR ANALYSIS 


ASG,T 2.U9B. , Reel No.* 
@END 

@ASG,A XQTANALYSIS. 
@ADD,PL XQTANALYSIS. 
@FIN 


The job control procedures XQTEASY and XQTANALYSIS are shown in Figures 
4.1-2 and 4.1-3. If a user is creating data inputs from a terminal, then 
it may be somewhat simpler to create new job control procedures similar to 
XQTEASY and XQTANALYSIS, but substituting his data input file names for 
EASYCARDS and NONSIMCARDS, respectively. If the same model is used for a 

series of runs, then only the analysis program is required for execution. 

✓ 

However, it is safer and also relatively inexpensive to execute both pro- 
grams when using the above job stream. Whenever the file read component is 
desired, the user must either substitute his file for FI or F2, or add the 
following job cards to XQTANALYSIS: 

@ASG,A MYFILE. 

@USE M, MYFILE. 

where MYFILE is the user time history file and M is a unit number between 
13 and 18. (See 7.38 for a discussion of the tape/file read component.) 

4.2 PROGRAM MAINTENANCE AND LIBRARY UPDATES 

Whenever the component library is updated, the user must compile the Fortran 
code and run the FILOAD program to furnish the model generation program com- 


*Used whenever TMY environmental tape data is to be input. 
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©HDG SlMWEST MODEL GENERATION 
OASG, AX MGABS . 

OASG,A M18. 

OUSE 18,M18. 

©ASG,T M7. 

OUSE 7,M7. 

©ASG,T SCRTCH8. 

OUSE 8,SCRTCH8. 

©DELETE, C SCRTCH9. 

OASG, UP SCRTCH9. 

OUSE 9,SCRTCH9. 

OASG,T SCRTCH10. 

OUSE 10, SCRTCH10. 

OASG,T SCRTCHll. 

OUSE 11, SCRTCHll. 

OASG,T SCRTCH12. 

OUSE 12,SCRTCH12. 

OASG,A EASYCARDS. 

OUSE 5, EASYCARDS. 

OXQT MGABS. EASY 
OASG,AX ASRO. 

OASG,AX ASS! . 

OADD,PL 9. 

OFREE 18. ,7. ,8. ,9., 1C. ,11. ,12. 


FIGURE 4.1-2 XQTEASY JOB CONTROL FILE 


OHDG SIMWEST ANALYSIS 
OASG,AX MAPANALYSIS. 
OADD,PL MAPANALYS IS. 
OASG ,AX ASABS. 

@ASG,AX FI. 

OUSE 11, FI. 

OASG, AX F2. 

OUSE 12, F2. 

OASG,T SCRTCH25. 

OUSE 25,SCRTCH25. 
©DELETE, C SCRTCH26. 
OASG, UP SCRTCH26. 

OUSE 26,SCRTCH26. 

OASG, AX NONSlMCARDS. 
OUSE 5, NONSlMCARDS. 
OXQT ASABS.NONS IM 
OXQT ASABS. NSMPPT 
OFREE 11., 12., 25., 26. 


FIGURE 4.1-3 XQTANA LYSIS, JOB CONTROL FILE 
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. nnn PRIMARY and xxx SECONDARY INDEPENDENT VARIABLE POINTS EXCEEDS THE 
zzz WORD STORAGE LIMIT FOR THE FOLLOWING TABLE. SOME DATA WILL BE 
LOST. 

The maximum amount of data allowed for each table is given in the 
Input Requirements List produced by the Model Generation program. 
Check that given data falls within this limit or for data card errors. 

5.2 DIAGNOSTIC MESSAGES FOR LIBRARY COMPONENTS 

A diagnostic message associated to a component is printed when a variable 
gets out of bounds during analysis. Adjustment of component parameters may 
be necessary. 

In component alphabetical order, these diagnostic messages are: 

AD: INPUT POWER xxxx TOO HIGH RELATIVE TO ADMITTANCE xxxx AND RATED 

VOLTAGE xxx 

ADMITTANCE POWER LOSS xxxx EXCEEDS INPUT POWER xxxx 

BA: POWER REQUEST xxxx EXCEEDS BATTERY CAPABILITY. CHECK VC, VO, AND RT. 

BN: BN INLET AIR MASS FLOW RATE xxxx GREATER THAN MAXIMUM ALLOWABLE xxxx 

CO: MAX ITERATIONS FOR COMPRESSOR EFFICIENCY. NP, XNP, RS = xxxx, xxxx, 
xxxx 

CS: CS STORAGE TEMPERATURE xxxx GREATER THAN ALLOWABLE xxxx 

CS MASS OF AIR IN STORAGE xxxx BELOW MINIMUM ALLOWABLE xxxx 
CS MASS OF AIR IN STORAGE xxxx EXCEEDS MAXIMUM ALLOWABLE xxxx 

ED: INPUT ERROR, DAY OF YEAR DY IS OUT OF RANGE 

TAPE INPUT ERROR OR EOF 

FL: FLYWHEEL POWER LOSS xxxx EXCEEDS CHARGING POWER xxxx 

FLYWHEEL LOSS xxxx EXCEEDS DISCHARGING POWER xxxx 
FLYWHEEL CLUTCH LOSS xxxx EXCEEDS MAXIMUM INPUT POWER xxxx 
FLYWHEEL CLUTCH LOSS xxxx EXCEEDS DELIVERABLE POWER xxxx 
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FLYWHEEL KINETIC ENERGY xxxx EXCEEDS CAPACITY xxxx 

FLYWHEEL KINETIC ENERGY xxxx FALLS BELOW MINIMUM REQUIREMENT xxxx 

GE: GENERATOR OUTPUT EXCEEDS RATED POWER 

HS: HS INLET MASS FLOW RATE xxxx OR OUTLET MASS FLOW RATE xxxx IS GREATER 

THAN MAXIMUM xxxx 

HS RESERVOIR VOLUME xxxx EXCEEDED MAXIMUM ALLOWABLE xxxx 
HS RESERVOIR VOLUME xxxx DROPPED BELOW MINIMUM xxxx 

HT: HT TURBINE CHARACTERISTIC PARAMETER OUT OF RANGE 

HT INLET MASS FLOW RATE xxxx GREATER THAN MAXIMUM DESIGN VALUE 

HX: HX EXIT TEMPERATURE xxxx GREATER THAN MAXIMUM ALLOWABLE xxxx 

IV: IV POWER LOSS xxxx EXCEEDS INPUT POWER xxxx CHECK RATED DC VOLTAGE VDC 

MB: WARNING-DIVISOR IN MB EQUALS 0., HAS BEEN SET = 1. 

MO: MOTOR INPUT POWER xxxx .GT. RATED INPUT POWER xxxx 

MOTOR SLIP xxxx EXCEEDS RATED POWER SLIP xxxx 
STATOR RESISTANCE xxxx OR DAMPING xxxx TOO HIGH FOR MOTOR 

PV: WARNING: INSOLATION OR TEMPERATURE AT CELL EXCEED RANGE 

RE: RE POWER LOSS xxxx EXCEEDS INPUT POWER xxxx 

RE, AC INPUT POWER xxxx TOO LARGE IN RELATION TO TRANSFORMER REACTANCE 
xxxx AND RATED AC VOLTAGE xxxx 

TA: FILE DATA OUT OF RANGE. INITIAL VALUE = xxxx ON UNIT xx 

TIME POINT PAST TABLE RANGE. LAST VALUE = xxxx ON UNIT xx 
READ ERROR OR END OF FILE ON UNIT xx 
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7.0 LIBRARY COMPONENT DESCRIPTIONS 


This section describes the mathematical algorithms and input/output struc- 
ture of the SIMWEST library components. Each component writeup contains a 
brief textual description of the algorithms, a mathematical expression 
summarizing its function, a list of input and output variables, a descrip- 
tion of the calculation sequence and logic used in the model, and the model 
code. A figure is provided which shows the nominal input and output 
connections, and the state variables of each component. 

There are a number of features and conventions in the component descrip- 
tions which require some elaboration. These are briefly summarized below. 

7a. INPUT/OUTPUT NAME LISTS 

A potentially confusing factor is the way port numbers on input parameters 
and output variables are designated. On the model generation input cards 
the name of the physical quantity and the port number are separated by a 
comma. For example, the power variable with port designation 1 is denoted 
P,l. To emphasize the distinction between the physical quantities and port 
numbers they are listed separately in the name lists of the component 
writeups. For example, P 1 in the name list denotes the power variable 
(or parameter) with port designation 1 even though in other parts of the 
text it may simply be denoted PI. 

Another convention in the name lists is that the alphabetic symbol 'O' is 
shown as 0 to distinguish this symbol from a zero. Elsewhere in the text 
symbols such as V0 may be referred to as VO. 

7b. INPUT PARAMETER SPECIFICATION 

All input parameters are associated with default values. Many of the 
parameters have default values denoted in the parameter description by the 
letter D. For example, in the Battery component the default value for 
terminal resistance, RT, is D = .001 ohms. All input parameters for which 
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a default value is not so specified have a default value of .99999. 
Default values are intended to enable users to put models together quickly 
by specifying a minimum of input data. Users need only specify detailed 
parameter values for those components of current interest. One must be 
careful using this approach since the operating characteristics and effi- 
ciency of a lOkw rated device may, for example, be quite different than for 
a lOOkw device. 

Any user-specified input parameter can be driven by one or two dimension 
table lookups using the FU and FV components. This enables the user to 
build more detailed models using time or other output variables to drive 
the tables. For example, if one needs to specify cost of peak load 
generation to the utility component as a function of peak load request, 
then one adds FU as an input to UT and specifies load request as an input 
connection to FU. The desired function table for FU is specified in the 
simulation input. 

It may be noted that not all of the components have maintenance or operat- 
ing cost inputs. Thus, whenever these costs are important, one can aggre- 
gate such costs and input lumped costs to the model. For example, the 
maintenance cost of the hydro storage system may include maintenance costs 
for the pump and turbine. 

7c. COMPONENT LOGIC 

In constructing SIMWEST components, we have adopted several conventions to 
aid communication with the logic components. All physical components 
distributing power are given two input parameters EF and MP (port 1) and 
two output variables EF and MP (port 2). The output EF is the product 
efficiency of all components in the distribution subsystem up to and in- 
cluding the given component, and MP is the maximum power deliverable at the 
output of the component. Each storage component has in addition a power 
request input denoted RE (port 1), a power request output denoted RE (port 
2), and a priority interrupt flag denoted INT. 
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Figure 7.0 shows the logic and physical variable connections for power flow 
in and out of a hydro reservoir. Power flows from the power divider to the 
pump at a rate not to exceed the request RE from HS. The HS request is 
computed by dividing the input maximum power by the input (or pump) effi- 
ciency EF. Hence, the maximum power flowing to HS cannot exceed RE*EF = 
MP. Similarly, the input request to HS is computed by the PA component so 
as not to exceed the maximum input power MP divided by EF (turbine effi- 
ciency). Hence, the power that flows to PA cannot exceed RE*EF = input 
maximum power. 

When the hydro reservoir is empty, the interrupt flag is turned on and the 
priority sequence is changed so that the reservoir is given access to power 
flowing into the divider. 

7d. UNITS 

Most of the SIMWEST components are coded in English units. However, SI or 
metric units were used to code the solar-photovoltaic components: ED, SO, 
FP, FO, and PV. This is generally not a problem since there are at most 
only a few interconnection variables between the solar-photovoltaic gener- 
ation components and other SIMWEST components, and units conversions are 
easily handled using an MA arithmetic component. (See for example the 
Fresnel Lens Model, section 9.3.) 
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FIGURE 7.0 SAMPLE CONNECTIONS FOR LOGIC COMPONENTS 
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Section 7.5 

Replaces pages 117 - 122 of the original document. 
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SYSTEM COST INPUTS 

• CAPITAL COSTS 

• MAINTENANCE COSTS 


UTILITY INPUTS 

• UTILITY ENERGY DELIVERED 

• VALUE OF UTILITY ENERGY 

• SURPLUS ENERGY SUPPLIED 



• SOUR ENERGY DELIVERED 

• VALUE 9F ENERGY DELIVERED 


• TOTAL LOAD DEMAND 


This component sums the capital, operating and maintenance costs of all 
system components. The total yearly cost TC is then computed using a fixed 
charge rate factor which represents depreciation, cost of money, insurance 
and taxes. 

a 

The total energy delivered to the loads plus surplus energy is then summed 
and yearly energy delivered TED computed. Cost of operation in mills is 


1 This component must be placed last in the model generation input file, 
i.e., just prior to the END OF MODEL command. 
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then given by 


CM 

System cost/kwh = TC * 1000. /TED 

Similarly, the value of energy delivered to the loads is summed minus the 
utility energy value and including the value of surplus energy, and fac- 
tored to give yearly energy value delivered VED. Energy value in mills is 
given by 

Load value/kwh = VED * 1000. /TED. 

Cost per value delivered is the ratio of the above two equations. 

In addition to the above cost calculations, percent of total load supplied 
by storage PCW, percent of load supplied by utilities PCU, and percent of 
energy surplused to the utilities PCS is computed. The total cost in mills 
to meet the load is then given by 

Load cost/kwh = (system cost/kwh * PCW + utility cost/kwh * PCU)/100., 

where 


Utility cost/kwh = value of utility energy * 1000. /uti lity energy 

delivered. 


Inputs 

Parameter/Port 

Description 

Units 

CR 

Capital charge rate 

%/year 

LE 

System life expectancy 

years 
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CM 


Common Block 

Inputs 

Description 

Units 

CC 

Total yearly capital costs 

$ 

CM 

Total yearly maintenance costs 

$ 

CO 

Operating and fuel costs over TMAX 

$ 

TMAX 

Simulation time interval 

hr 

VDE 

Value of energy delivered (including surplus) 

$ 

TDE 

Solar energy delivered (including surplus) 

kwh 

TLD 

Total load demand 

kwh 

UTV 

Value of utility energy 

$ 

UTD 

Utility energy supplied 

kwh 

SPD 

Surplus energy supplied 

kwh 

Outputs* 


Total yearly costs (TC) 

$ 


Yearly energy delivered (TED) 

kwh 


Cost of energy per kwh 

mills 


Yearly value delivered (VED) 

$ 


Cost per value delivered 
Percent of load supplied by 

- 


Storage (PCW) 

- 


Utility (PCU) 

- 


Surplus energy load factor (PCS) 

- 


Total load cost per kwh 

mills 


I 


f 


Printout only occurs when simulation is completed, 
variable symbol is required. 


Thus no output 
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no 

o 


SUBROUTINE CM 


ENTRY POINT 000213 


STORASE USED: C00EI1) 000226; OATAIOI 000332; BLANK C0MH0N(2I 000000 

COHNON BLOCKSt 

0003 COST 000011 
C004 CIHPL 000001 
LOOS C1IME 000001 
0006 CSIhUL 000010 


EXTERNAL REFERENCES I BLOCK , NAME I 

C0J7 NWOUS 
L010 NIOiS 
C.C11 NERR3S 


STORASE ASSIGNMENT (BLOCK* TYPE* RELATIVE LOCATION* NAME) 


0001 


0G002C 

100L 

0000 


000016 

200F 

0000 


00003 S 

300F 

0000 


000106 

40 OF 

0300 


000211 

500F 

0033 

R 

OOQLOO 

cc 

GPuO 

R 

000003 

CCY 

0003 

R 

000001 

CMA 

0003 

R 

000002 

CO 

0000 

R 

000002 

COY 

ccoo 

R 

CCOulS 

CPHVH 

cooo 

R 

C00C11 

CPV 

0006 


000003 

0UM 

0000 

R 

000005 

EOE 

0004 

I 

000000 

IMPL 

CCQO 


CLC315 

INJPS 

ooco 

I 

000006 

I VOE 

0000 

I 

000001 

LLE 

ooco 

R 

rooji 2 

PCD 

0003 

R 

000314 

PCS 

COOO 

R 

0LGU13 

PCU 

L0C3 

R 

03001C 

SPD 

0003 

R 

0CCD04 

TDE 

0005 

R 

000130 

TIME 

0003 

R 

030305 

TLD 

L0u6 

R 

0L0li07 

TMAX 

COOO 

R 

oocooo 

TMAX 1 

COCO 

R 

000004 

TOY 

0000 

R 

000 JO 7 

TOYN 

0003 

R 

000307 

UTO 

0003 

R 

0C0C06 

UTV 

C003 

R 

000003 

VOE 

0000 

R 

oonoio 

VDEN 










SS 

Q 


1 


P 


o 2 
c S 
> G) 
F nn 



00100 

1* 

COST 

000000 


aoioi 

2* 


SUBROUTINE CM(0UHM,FCR*LET 

OD&OOO 

CO 

00101 

3* 

C 


ooooco 

o 

C/1 

00101 

4* 

C 

PURPOSE SUMMARIZE WIND ENER6E STORAGE COSTS AND LCVELIZED 

000000 

J0131 

5* 

c 

ENERGY COSTS PER KWH* 

030000 

o 

00101 

6* 

c 


000300 

CO 

o 

00101 

7* 

c 

WRITTEN BY A*W. WARREN VERSION 1* MAY 1977 

000000 

001D1 

6* 

c 


QOGOOO 

1 

no 

00101 

9* 

c 

INPUT PARAMETERS 

OOuOGO 

UOIOI 

1G* 

c 

0C3P3D 

70 

to 

CC131 

11* 

c 

FCR - FIXED CHARGE RATE FACTOR INCLUDING DEPRECIATION* 

D3COOO 

< 

00101 

12* 

c 

MONEY COST* INSURANCE, AND TAXES* PER YEAR 

ooonao 

• 

CC101 

13* 

c 

LE - SYSTEM LIFE EXPECTANCY , YEARS 

o^opao 


00101 

14* 

c 

TMAX - SIMULATION TIME, HR 

00GP00 


0D101 

15* 

c 

CC - TOTAL YEARLY CAPITAL COSTS, S 

000000 


U3101 

16* 

c 

CM - TOTAL YEAPLY MAINTENANCE COSTS, S 

ooonoo 


0 0131 

17* 

c 

CO - TOTAL OPERATING AND FUEL COSTS OVER TMAX, % 

O r, DDOO 


uOlOl 

ie* 

c 

VDE f VALUE OF ENERGY DELIVERED OVER TMAX, S 

OOuPGO 


bClOl 

19* 

c 

TOE - TOTAL ENERGY DELIVERED OVFR TMAX, KWH 

oourao 


00101 

20* 

c 

TLD - TOTAL LOAD DEMAND OVER TMAX, KWH 

cncoao 
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00101 

21* 

C 

UTV - VALUE OF UTILITY ENERGY SUPPLIED LESS SURPLUS VALUE 

00101 

22* 

C 

UTD - TOTAL UTILITY ENERGY DELIVERED, KWH 

00101 

23* 

C 

SPD - TOTAL SURPLUS ENERGY SUPPLIED TO UTILITY, i 

00101 

24* 

C 


□ 01 03 

25* 


COMMON /COST/ CC, CPA, CO, VDE.TDE.TLD, UTV , UTD ,SPD 

OS? -~4 

26* 


COUPON /CIMPL/IHPL /CTIME/ TIME /CSIHUL/ 0UM!7),TMAX 

00105 

27* 


REAL LE 

coins 

28* 

C 

INITIALIZATION 

00105 

29* 

C 


0C106 

30* 


IF1IMPL.GT.3IG0 TO 100 

00110 

11* 


DUHH=0.0 

join 

32* 


CC = 0 . 

□ Cl 12 

33* 


CPA = 0. 

00113 

34* 


CO = 0. 

031 14 

35* 


VDE= 0. 

00115 

36* 


TDE= 0. 

00116 

37* 


TLD- 0. 

00117 

38* 


DTV=ii. 

001 20 

39* 


UTO-Q. 

00121 

40* 


SPD=C. 

□ C122 

41* 


THAXi= TMAX*. 99999 

00122 

42* 

C 


CC122 

43* 


100 IFIT1PE.LT.TMAXI IRE TURN 

00125 

44* 


IFllPPL.LE.llRETURN 

30125 

45* 

c 

, 

00125 

46* 

c 

COST SUMMARY OUTPUT 

U012S 

47* 

c 


00127 

48* 


LLE = LE 

□ 0130 

49* 


6RITE 1 6 ,200 1LLE 

00133 

50* 


200 FORMAT <1H1,35X,39H SOLAR/WIND ENERGY STORAGE COST SUMMARY // 

00133 

51* 


1 1H ,4GX,I2,17H YEAR LIFE CYCLE 1 

00133 

52* 

c 


03134 

S3* 


COY = C0*8762./TMAX 

□ 0135 

S4* 


CCY = CC*LE*FCR*.01 

00136 

S5* 


TOY : COY + CPA ♦ CCY 

30137 

56* 


WFITE16.30O>CCY,CMA,COV,TOV 

CC145 

57* • 


300 FORMAT!//// 30X.22H0 YEARLY SYSTEM COSTS/ lH«,29X,lH*/ 1H-.42X, 

30145 

58* 


1 I2HCAPITAL COST ,12X.F8.U,2H % / IN ,42X, 17H (INCLUDING FIXED • 

0014b 

59* 


2 SHChARGES) / 1H0,42X,16HFIXED 0 ♦ M COST, 8X,F8.0,2H S /1HQ , 

00145 

60* 


3 42X,21H0P£PATINL ♦ FUEL COST, 3X,F8.0,2H % / 1H0,42X,5HT0TAL , 

□ 0145 

61* 


4 19X,F3.C,2H 5 1 

□ 0145 

62* 

c 


□ 0146 

63* 


EOE = TOE * 8762. /TMAX 

00147 

64* 


1 VOE : VDE * 8763. /TMAX 

00150 

65* 


TOYN = TOY* IQPO./ EOE 

00151 

66* 


VOEN = VDE* 1000,/ TOE 

03152 

67* 


CPV = TOYN / VDEN 

00152 

68* 

c 


00153 

69* 


UPlTE(6,40OIE0E,T0YN,IV0E,VDEN,CPV 

0C152 

7C* 


400 FORMAT!//// 30X,16Hn ENERGY DELIVERED / 1H*,29X,1M* / 1H-, 

00162 

71* 


1 42X.16HENEPGY DELIVERED, 7X,F9.3,4M KWH / 1HP.33X.5GOH*) / 

00162 

72* 


1 1H ,33X,1M*,48X, 1H* / 

00162 

73* 


2 1H ,33X,1H*, 8X , 19HENEPG Y COST PER KWH, 7X,F6.1,9H MILLS * / 

□ 3162 

74* 


2 1H ,33X,1H*,49X,1H* / 1H , 

□ 0162 

75* 


3 33X,1CJ(5M*****) / 1H0,42X,?SHVALUE OF ENERGY OEL IVERED , 17 , 

00162 

76* 


4 2H S / 1H , 42*, i2HI VALUE OF FUEL SAWED) / 1H D ,42 X , 2 JHENERG Y VALUE 

30162 

77* 


5 PER KkH, 6X,F 6,1 ,6H MILLS / 1 HO, 42 X , 24HC0ST PER VALUE DELIVERED, 


oooooo 

030C00 

ononco 

oooooo 

03000? 

onuroo 

oooooo 

oooroo 

oooooo 

OlliOOO 

000002 

00DP03 

DPCD04 

oocros 

000006 

000007 

000010 

OOtiPil 

CPCP12 

000^13 

000014 

300014 

000020 

000025 

000026 

0 r, 0026 

C0u026 

000040 

000047 

000055 

000355 

000055 

Q00D55 

oroo6i 

0 , 'CCfc6 

C 10102 
000102 
000102 
000102 
000102 
000102 
000102 
^00106 
0PC120 
000124 
300130 
000130 
000132 
OOC 144 
000144 
000144 
000144 
0^0144 
000144 
000144 
0 n C 144 
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10162 

78* 


6 2X.F6.2) 

000199 

10162 

79* 

C 


0*50199 

>0161 

00* 


PC0 = IT0E-SPD»*100./TL0 

000199 

1816*1 

81* 


PCU = UTO*1UO./TLO 

0*50150 

1016$ 

82* 


PCS= SPDMOO./TLD 

0 *58159 

10166 

83* 


CPKWH: IT0VN*(T0E-SP0> ♦ UTV*1003.)/TLD 

030160 

J0167 

86* 


6RITE I6,50C»PCO »PCU(PCS»CPNUH 

0*50167 

•0175 

as* 


500 FORMAT!//// ICXfllHO LOAD rACTOR / IH«,29X, 

3 3u2?0 

.0175 

86* 


1 1H* / 1H-.42X, 

0PC2C3 

•C175 

8 7* 


1 26HPEPCENT OF LOAD SUPPLIED , F6.1, 2M / 1H ,92X* 28NBY TOTAL S 

0*50200 

•C175 

88* 


20LA.R SYSTEM / 1H0,*?X,2*HPERCEKT OF LOAD SUPPLIE0 ( 2X.F6«1 / 

0*10200 

Cl 75 

89* 


2 1H t 91X,UH PY UTILITY / 

0*50200 

0175 

9C* 


3 1hO ( 92X ( 26HPcRCENT OF SOLXR ENERGY « F6.1 / 

000230 

0175 

91* 


3 lH ,*»2X»9hSURPLUS£D / 

030230 

0175 

92* 


9 iHj,42X # i JHCOST TO MEET LOAD , 3X.F6.1.6H MILLS/ 

0*10200 

Cl 75 

93* 


5 1H t * 2 X , 29H( SOLAR ♦ UTILITY) / 1 H 1 > 

0**0?30 

0175 

9** 

C 


0*'C230 

0176 

95* 


RETURN 

a**c?oo 

•0177 

96* 


END 

0*50225 
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ED 

I 

V 7.7A ENVIRONMENTAL DATA (TMY TAPE) 

Day-of-Year (DY) 

Time-of-Day (TD) 

This component reads data values from the Typical Meteorological Year 
(TMY) tapes or data with a similar format structure such as the University 
, of Wisconsin insolation and environmental data tape or the SOLMET tapes. 

Only one ED component is allowed per model. (Unit 2 is reserved for the 
input tape.) The file structure assumes hourly recorded data with one 
record or card image per hour of data. Twenty-four hourly records are read 
into core at a time and linear interpolation is used to obtain the output 
values at the current simulation time. The component TI is used to supply 
the time inputs DY and TD. Standard outputs with the TMY tape are direct 
and global solar insolation, dry bulb temperature, and vyind speed. For 
non-standard outputs or non-TMY format tapes the user may specify the input 
format to read one to eight data variables. The following limitations 
apply in this case: 

( 
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1) 


ED 

Time information is decoded in integer month (1-12), day (1-31), and 
hour (0-24) format. 

2) Output variables are decoded in F or E format, even if recorded in 
integer format. 

3) Where data is missing, fill in with 9‘s is assumed. The code checks 
for certain 9 fill values, namely 99., 999., 9999., and 99999. If any 
one of these values is read, then the corresponding data input is 
replaced with 0. or the previous value, depending on the sign of IND. 
(However, one must use FN.O format N=2,3,4,5 for this option and a 
scale multiplier if necessary to obtain the desired exponent.) 

Inputs/Port Description Units 

NST Number of tape blocks to skip at start* 

NX Number of output variables (default = 4, max = 8) 

IND Indicator function: 

0 = no read 

±1 * standard format and units (default) 

±2 = user-specified format and units 

IND>0 sets missing data = 0 

INDO sets missing data = previous value 


T 

For the TMY tapes we may compute NST from the station number (NSTAT) shown 
in table 7.7A and the start day (DSTART): 

NST = (NFILE-1 )*365 + DSTART-1 
where NFILE = /NSTAT if NSTAT < 13 

\NSTAT-13 otherwise 
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Inputs/Port 

(cont 'd) Description 

TS* Time shift of data (default = -0.5) 

TD Time of day (0-24) 

DY Day of year (1-365) 

Ml Units multiplier for XI (default = 1) 


M8 Units multiplier for X8 (default = 1) 

A1 Units addition factor for XI (default = 0) 

• • 

• * 

• • 

A8 Units addition factor for X8 (default = 0) 



Units 

hours 

hours 


Compensation term since solar radiation data is an integrated total over 
the observation interval. 
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Outputs/Port 


ED 

Description Units 

XI 1st output variable 

p 

(IND =±1: beam radiation in w/m 
X2 2nd output variable 

o 

(IND =±1: global horizontal radiation in w/m ) 

X3 3rd output variable 

(IND =±1: dry bulb temperature in °C) 

X‘ 4th output variable 

(IND = ±1: wind speed in m/s) 


X8 8th output variable 


Format Specification 

A user-specified format may be input in order to select non-standard en- 
vironmental outputs or to read a tape other than the TMY insolation tape. 
The following sequence of data cards is recommended for insertion in the 
model generation input following the MODEL DESCRIPTION command: 

FORTRAN STATEMENTS 
DIMENSION FMT ( 12 ) 

COMMON/READER/ N, FMT 
DATA FMT/72H. . . ) 

1 /N/NN/ 

where the format specification contains up to 71 characters inserted after 
* 72H ' and followed by and NN = the number of characters per data 
record. 


HOD 
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ED 

The format specification must conform to the following rules: 

1) The first two words read are station and year identifying informa- 
tion. These words must be either A format or nH format with up to six 
characters for station and two characters NN for year 19NN. 

2) The next three words are two-digit integers containing month (1-12), 
day (1-31), and hour (0-24) information. 

3) The next one to eight words specify the location of the output vari- 
ables XI... X8 and must be given in F or E format. 

NOTE: The tab or column spacing control T may be used to read data from 

files which are not ordered as in 1) to 3), e.g., T71, A5, Tl, A2,...). 

For example, the standard TMY tape format specification is 

Station Yr-Mo-Dy-Hr Beam Rad. Global Rad. Temp Wind 
A5, A2,3I2,11X, F4.0 ,26X, F4.0.45X, F4.1,7X, F4.1) 

and N = 132. 

The general format for variables on the TMY tape is summarized in Figure 

7.7A. 
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-SOL A R 

E-A P 

I A T 

.1 r- s 

O R S E R V 

UT 
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W6A.N 


SO LAX 

LST 

ETR 

RAC l AT 10 / VALVES KJ/n* 

Sj 

STM 


TIME 


TIME 

KJ/» 2 

D 

D 
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OUTSaE 
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# 

BC 
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DY 

HEm5T 
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I 
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xxxxx 
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xxxx 

XXXX 

XXXX 

1XXXX 

1XXXX 

1XXXX 

1XXXX 

1XXXX 

1XXXX 

1XXXX 

1XXXX 

1XXXX 

xx 1 

002 

003 




004 

101 

102 

103 

104 

103 

106 

107 

10S 

109 

110 

111 


S I 

r r 

l C E 

M E 

r E 0 R 0 

L 0 C I C A 

L 0 

9 S 

E R 

V A 1 

nn 

r 

0 

c 

SKY 

VSBY 
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PRESSURE 
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VINO 

1 

s 

> 

E 
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IcPs 
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T 

0 

N 
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STA- 

DRY 

DEW- 
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s 

0 

p 

0 

T 
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LEVEL 

TION 
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PT. 

I 
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u 

X 
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A 
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XX 

xxxx 

IXXXX 

xxxx 

XXXXXXXX 

XXXXX 

XXXXX 

XXXX 

XXXX 

XXX 

XXXX 

XX 

XX 

X. 

*11 


TAPE 

RECORD 


FIELD NUMBER 

POSITIONS 

DESCRIPTION 

002 

01-05 

WBAN STATION NUMBER 

003 

06-15 

SOLAR TIME ( YR ,M0 , DAY , HOUR .MINUTE ) 

004 

16-19 

LOCAL STANDARD TIME {HR AND MINUTE) 

101 

20-23 

EXTRATERRESTRIAL RADIATION 

102 

24-28 

DIRECT RADIATION 

103 

29-33 

DIFFUSE RADIATION 

104 

34-38 

NET RADIATION 

105 

39-43 

GLOBAL RADIATION ON A TILTED SURFACE 

106 

44-48 

GLOBAL RADIATION ON A HORIZONTAL SURFACE- 
OBSERVED DATA 

107 

49-53 

GLOBAL RADIATION ON A HORIZONTAL SURFACE- 
ENGINEERING CORRECTED DATA 

108 

S4-58 

GLOBAL RADIATION ON A HORIZONTAL SURFACE- 
STANDARD YEAR CORRECTED DATA 

109,110 

59-68 

ADDITIONAL RADIATION MEASUREMENTS 

111 

69-70 

MINUTES OF SUNSHINE 

201 

71-72 

TIME OF COLLATERAL SURFACE OBSERVATION 
(LST) 

202 

73-76 

CEILING HEIGHT (OEKAMETERS) 

203 

77-81 

SKY CONDITION 

204 

82-85 

VISIBILITY (HECTOMETERS) 

205 

86-93 

HEATHER 

206 

94-103 

PRESSURE (KILOPASCALS) 

207 

104-111 

TEMPERATURE (DEGREES CELSIUS TO TENTHS) 

208 

112-118 

HIND (SPEED IN METERS PER SECOND TO TENTHS) 

209 

119-122 

CLOUDS 

210 

123 

SNOW COVER INDICATOR 


140F 


Figure 7.7A TMY Tape Format 
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A complete description of the available data, and the meaning of the 
recorded outputs, is contained in the SOLMET user's manual [3]. The TMY 
tape was derived from SOLMET tapes of the 26 stations with rehabilitated 
solar radiation data, and has the same format as the SOLMET tapes except 
that tape deck number and detailed cloud data have been omitted. Table 
7.7A shows the identity and location of the 26 stations on the TMY tape. 

Calculation Sequence 
If IND = 0 Return 

1) INITIALIZATION (first pass only) 

• Set defaults and initialize LTD = -1 

• Skip NST blocks to position the file 

• Read first data block and write out identification informa- 
tion. (Error exit to 6)) 

t Go to 4) 

2) Table Interpolation for Output (DY = DYF) 

• If OY > DYF go to 3) 

§ If DYF >DY go to 5) 

• If LTD = TD return (LTD = last time C(I,J) was accessed) 

§ X(I) = TBLU1 (TD, TO, C(1 , I ) ,0,24)*M( I )+A( I ) I = 1,...NX 

§ LTD * TD 

• Return 
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3) Read One or More Data Blocks (DY > DYF ) 

• Read DY-DYF data blocks. (Error exit or EOF exit to 6)) 

4) Decode Using Specified Format 

• Decode day-of-year (DYF) and time information (TO) and put 
output variables in array C( I » J) 1=1,24 and J=1,NX. Check 
for missing data values in C(I,J). 

• Go to 2) 

5) Backspace the File (DYF > DY) 

• Backspace and read first data block 

• Decode day-of-year (DYF) 

• Go to 4) if DYFiDY. Otherwise print diagnostic and stop. 

6) Read Error or EOF Encountered 

t Print diagnostic and stop. 


140H 


BCS 40180-2 Rev. 


ED 



TABLE 

7.7A TMY TAPE STATIONS AND LOCATION 


STATION 

NUMBER 

WBAN 

IDENTIFIER 

STATION 

LATITUDE 

LONGITUDE 

1 

3927 

Fort Worth, Texas 

32°50' 

97°03 ' 

2 

3937 

Lake Charles, Louisiana 

30°07 ' 

93°13' 

3 

3945 

Columbia, Missouri 

38°49 ' 

92°13' 

4 

12832 

Apalachicola, Florida 

29°44' 

84°59 1 

5 

12839 

Miami, Florida 

25°48' 

80°16 1 

6 

12919 

Brownsville, Texas 

25°54' 

97°26‘ 

7 

13880 

Charleston, South Carolina 

32°54' 

80°02 ' 

8 

13897 

Nashville, Tennessee 

36°07 1 

86°41' 

9 

13985 

Dodge City, Kansas 

37°46' 

99°58' 

10 

14607 

Caribou, Maine 

46°52' 

68°01' 

11 

14837 

Madison, Wisconsin 

43°08 ' 

89°20' 

12 

23044 

El Paso, Texas 

31°48' 

106 °24' 

13 

23050 

Albuquerque, New Mexico 

35°03’ 

106°37‘ 

14 

23154 

Ely, Nevada 

39°17 ' 

114°51' 

15 

23183 

Phoenix, Arizona 

33°26 ' 

112°01 1 

16 

23273 

Santa Maria 

34°54* 

120°27 ' 

17 

24011 

Bismarck, North Dakota 

46°46 1 

100°45 ' 

18 

24143 

Great Falls, Montana 

47°29' 

111°22 * 

19 

24225 

Medford, Oregon 

42°22* 

122°52‘ 

20 

24233 

Seattl e-Tacoma, Washington 

47°27 1 

122°18' 

21 

93193 

Fresno, California 

36°46 1 

119°43' 

22 

93729 

Cape Hatteras, North Carolina 


75°33‘ 

23 

93734 

Washington, D.C. 

38°59‘ 

77°28' 

24 

94701 

Boston, Massachusetts 

42°22' 

71°03' 

25 

94728 

New York, New York 

40°47 1 

73°58' 

26 

94918 

North Omaha, Nebraska 

41°22' 

96°01 ' 
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ENTRY POINT 000644 


o 

C_. 


SUBROUTINE ED 


STORAGE USEOt CODE (1) 030704; DATA(O) 001624; BLANK COMMON(2> 000000 

, COMMON BLOCKS t 

0003 READER 00091S 

UC04 CIMPL 3000G3 


EXTERNAL REFERENCES (BLOCK, NAME! 

GC35 NTRAN 
0036 NDCODS 
uCQ7 TBLUl 
0010 NI03S 
OGll NI02S 

0012 NhDUS 

0013 NS70PS 

0014 NERR3S 


STORAGE ASSIGNMENT (BLOCK, TYPE, RELATIVE LOCATION, NAME) 



OOJ1 


0L024G 

10QL 

OOC1 


000101 

147G 

0001 


000125 

1636 

0001 


000141 

1736 

0001 


000138 

2L 

C0U1 


CG024C 

2C0L 

0001 


0001SC 

201G 

0001 


000270 

240G 

0001 


C00365 

2566 

0001 


000424 

2706 

wool 


0C0312 

30UL 

OOCO 


001506 

306F 

0001 


000345 

400L 

cool 


030512 

500L 

0001 


000605 

507L 

0CG3 


00151 6 

5C8F 

LOCI 


C0C614 

6CCL 

occc 


001530 

608F 

coco 

R 

000330 

A 

0000 

R 

000324 

AA 

LOCO 

R 

001351 

B 

COCO 

R 

000024 

C 

0000 

R 

C31504 

CIJ 

cooo 

R 

C01361 

CL 

3000 

R 

OC 1421 

DM 

COCO 

R 

0Lli.CC 

DYF 

UPG3 

R 

300001 

FMT 

oroo 

R 

PC 143 5 

FMTP 

coco 

R 

ocojie 

FRMT 

0003 

I 

021472 

I 

LCjO 

I 

0C1344 

IB 

00C4 


COCJOl 

ICNT 

ocar 

I 

001501 

10 

CP04 

I 

roroco 

1MPL 

0000 


0C1S71 

IN JPS 

crjo 

I 

OC Is 71 

I NX 

coao 

I 

P01475 

I REWiN 

0004 


P0DD02 

ITEST 

cooo 

I 

031503 

11 

ocoo 

I 

001473 

J 

coco 

I 

001502 

31 

0000 

I 

001505 

L 

oc-oc 

R 

001470 

LTD 

0033 

I 

001477 

LI 

0303 

R 

001463 

M 

DC 33 

I 

DC OLuC 

N 

UOCO 

1 

001474 

NO 

0000 

I 

r 0147b 

N 1 

ccoo 

R 

001452 

OFFSET 

0007 

R 

OCOOOO 

TBLUl 


cOCO R 0L1371 TO 


00100 

1* 

CEO 



000030 

00101 

2* 


SUBROUTINE 

E0(X,X2,X3,X4,X5,X6,X7,X8,NST ,NX,IN0,TS,T0,0Y, 

003000 

00101 

3* 


1M1,M2,M3,M4,MS,M6,M7,M8,A1,A2,A3,A4,A5,A6,A7,A8} 

000000 

00101 

4* 

C 



030000 

ooiai 

5* 

C 

PURPOSE 

THIS COMPONENT READS THE TYPICAL METEOROLOGICAL 

00CM3C0 

00101 

6* 

C 


YEAR TAPE WITH A STRUCTURE SIMILAR TO THE 

03G030 

com 

7* 

C 


SOLMET DATA TAPE. USER MAY SPECIFY FORMAT FOR NON- 

oponoo 

0 3101 

8* 

C 


STANDARD TAPES 

0 70030 

CC101 

9* 

C 



ooccoo 

00131 

10* 

C 

WRITTEN 

8Y Y.K.ChAN, 1C-5-7P, VERSION l 

03C030 

00101 

11* 

c 



O'UiOOO 

uOlOl 

12* 

c 

METHOD 

TWENTY FOUR HOURLY RECORDS ARE READ INTO CORE 


uciai 

13* 

c 


AT A TIME AND LINEAR INTERPOLATION IS USEO TO 

COG3O0 



/ 
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0C140 

71* 



A(3):A3 


000060 

00141 

72* 



A (4 )=A4 


000062 

00142 

73* 



A 15 )-A5 


000064 

OC143 

74* 



A (6 1=A6 


CD0966 

00144 

75* 



AI7>=A7 


000070 

uC14S 

76* 



Aldl=A8 


G00072 

CC146 

77* 



00 11 1=1, INX 


oocini 

0 0151 

78* 



IF ( * ( I ) ,E0 • *99999) Ht 1 1 -1 • 



00153 

79* 


11 

1FJM1I.E0.. 999991 A ll» = 0. 


000105 

00153 

8C* 

C 




DOC IPS 

00153 

81* 

C 


SET DEFAULT TAPE RECORO FORMAT 

TO STANDARO 

003105 

uCISfc 

82* 



IFtABSlINOl.GT.l.ClIGO TO 2 


OOull 3 

3016U 

33* 



Mil 1 = 1 »/3 ,6 


000117 

00161 

84* 



11121 = 1 ./3.6 


003121 

0C162 

85* 



00 3 1=1,12 


000125 

00163 

86* 


3 

FMT(1»=FRMTIII 


OOC 125 

00167 

87* 



N=1 32 


OOC 127 

00170 

88* 


2 

CONTINUE 


OOC 132 

0C171 

89* 



FMTPO >=GFFSETtl> 


O n Gl 3? 

uC'l 7 2 

90* 



DO 5 1=1,12 


00C141 

00175 

91* 


5 

FMTPII*1»=FMTIII 


000141 

30177 

92* 



LTG=-l. 


■ 3001*3 

C32 0C 

93* 



DO 1C 3= 1 , I NX 


GOC153 

00203 

94* 


20 

CUJi:C'. 


0^155 

00203 

95* 

c 


. 


000150 

GO* 33 

96* 

c 


POSITION THE FILE 


000150 

00233 

97* 

c 




000159 

00205 

98* 



CALL N Tft AN 12,10) 


00G1S1 

002 06 

99* 



fcO=NST«.CCl 


0OJ155 

GS237 

ICO* 



CALL NTRANI2,7 , NO, 221 


C 00 1 66 

G02 37 

1GI* 

c 




09C166 

0 02 07 

102* 

c 


READ FIRST DATA BLOCK 


0.00)66 

u0207 

103* 

c 




00U166 

C 02 10 

1C4* 



IREblN=C 


00G174 

01211 

lu5* 



M = 4*N 


09&175 

Ou2 12 

lufc* 



CALL NTRANI?,2,N1,AA,L1»22I 


C9U290 

30213 

1C7* - 



IFtLl.LT.OIGO TO 630 


oooPin 

00215 

lue* 



DECODE (N,FMTP P AA 111 HB,B 


00C213 

0C221 

1C9* 



bRITE(C, 398)1 Bill, IP(2> 


0OQ227 

U022S 

11C* 


308 

F0RHAT(1HC,?X,15HE0C STATION ID=, A5 ,10X , 7MYEAR 

19, A2) 

09J236 

UC226 

111* 



00 TO 40u 


Q0C'?36 

U32 26 

112* 

c 




000236 

00226 

113* 

c 




000236 

C0227 

114* 


100 

CONTINUE 


00G24D 

OC227 

115* 

c 




00u?40 

03230 

116* 


200 

CONTINUE 


000240 

0 02 3 0 

117* 

c 




003240 

j C2 SO 

118* 

c 


INTERPOLATION FOR OUTPUT IF CURRENT 

OAV OF YEAR HAS 

000240 

GC23C 

119* 

c 


BEEN LOCATED 


090240 

002 3 w 

120* 

c 




00C240 

C0231 

121* 



IF(Or.&T.(OVF«.lllGO TO 3CJ 


G0C240 

uC2 3 3 

122* 



IFIDY.LT.IDTF-.lllGO TO 500 


00i.->44 

□ 0235 

123* 



IF(L10.EC.T01RETURN 


093254 

CC237 

124* 



DO 2C1 1 = 1, INX 


CP1279 

□ 0242 

125* 


201 

X(I >=TFLU1(TD,T0,C( 1 , I 1 , 0 , 24I*M II 1*A|I) 


QOO270 

00244 

126* 



LTO=TO 


OOC 0Q4 

0C?45 

127* 



RETURN 


0O03D6 
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00245 

12B* 

C 





Q3G336 

302*6 

129* 


.300 

continue 



Q7C312 

302*6 

130* 

c 





000312 

□ 02*6 

131* 

c 


IF CURRENT DAY OF YEAR HAS NOT BEEN 10CATE0, READ 

HORE TAPE 

000312 

30246 

132* 

c 





09C712 

002*7 

133* 



10 = DY -DYF -.9 



01G712 

0 32 SC 

134* 



CALL NTBAN) 2, 7, ID, 22) 



030723 

00251 

135* 



CALL NTPSNI 2, 2, Nl, AA, L 1, 22) 



C3G7 II 

002 5 2 

136* 



1F< LI .a. 0) GO TO 600 



oicm 

CT2;5 2 

137* 

c 





0HC341 

3C252 

136* 

c 


DECODE DAT* AND TINE OF DAY 



030341 

0Q2S2 

139* 

c 





Q n 034 I 

00254 

1 4 C* 


4 00 

CONTINUE 



0 3D 34 * 

0C255 

141* 



CO 4C? 1=1,24 



030345 

0 0260 

142* 



Jl = HI-l)*N)/6 



030771 

0C261 

143* 



11=11-1 >*N-6*J1 



010374 

00262 

1 4 4 ♦ 



FPTP ( 1 ) =OFF SET 1 1 1* 1 ) 



0 n u377 

n 02 6 2 

145* 



DECODE (N.FMTP, A A l Jl*l) IIB,B 



03C-402 

00267 

146* 



DO 401 J= 1 , 1KX 



010424 

CC272 

147* 



C(I,U>=PIJ) 



030424 

00273 

143* 



cij=ci i ,j » 



030425 

00274 

149* 



lFMCIu.L0.59. 1 .OR.fCIJ.EO ,999 • )«0R . (Cl J .E0.9999. ).0R. 


030426 

u5274 

ISC* 



1 1 Cl J. EL. 999 99. ))C(I,J)=CL!J) 



Ci 1Q426 

00276 

151* 



1F( IND.LT.2.)CLC J)=CI I ,J) 



030451 

00310 

152* 


401 

CONTINUE 



D30470 

00302 

153* 



T 0 ( I ) = IE l 5 ) *TS 



03D470 

00303 

154* 


402 

CONTINUE 



1 03u5Q0 

OC3DS 

155* 



L = IBID 



03U500 

00306 

156* 



DYF = 16«4) *0(1)0 



032502 

00337 

157* 



CO TO 2C'C 



310510 

003 07 

156* 

c 





030510 

0 03 lu 

159* 


soc 

CONTINUE 



015512 

OC310 

16G* 

c 





G0U512 

CC31C 

161* 

c 


IF DAY OF YEAR ON TAPE IS PAST 

CURRENT OAY OF YEAR, 

J"051? 

0031 J 

16 2* 

c 


BACKSPACE TAPE. 



012512 

03311 

163* 



IFIIREWIN.GT,0)GO TO 507 



Q~041? 

00313 

164* ' 



2 0 = C Y-DYF-I * 1 



0iU«14 

3C314 

165* 



CALL NTB*N)2,7,ID,2?) 



03C526 

0331a 

166* 



I RE U I Nil 



030534 

C 03 1 6 

167* 



CALL NTRA N 12,2, N1,AA,L 1,22) 



0 n u e 26 

00317 

ue* 



IFILi.LT.OGO TO 600 



030546 

00321 

169* 



FPTPtl )=0FFSET)1) 



01Q551 

00322 

170* 



DECODE (N,FnTP ,AA 11 ) )IB ,B 



030553 

003 26 

171* 



L = IPI3) 



0-10567 

0 03? 7 

172* 



DYF = IP ( 4 ) *DM(L> 



01L571 

03330 

173* 



IF(DYF .LT.tDY+.l ))G0 TO HOC 



0°C577 

03332 

174* 


507 

4PITE»6,S' - 6) 



03C635 

00334 

175* 


SOB 

F0k‘'AT(1HC,47HEDC INPUT epror, DAY OF 

YEAR OY IS OUT OF 

RANGE ) 

C306H 

03335 

176* 



STOP 



CiDfll 

00336 

177* 

c 





C10611 

00335 

176* 

c 


IF ERROR IN READ, PRINT DIAGNOSTICS 


030611 

3 0336 

179* 


60C 

URITLCfc ,606) 



030614 

003*0 

1 5 0* 


606 

F ORB AT )lhO,27HEDC TAPE INPUT ERROR OR 

EOF ) 


030623 

00341 

161* 



STOP 



CiOL 20 

00342 

1p2* 



END 



Q10703 
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Section 7.8A - FO and 7.8B - FP 

Delete pages 151 and 152 of the original document 
and insert revision pages 151 - 152N between 
pages 150 and 153. 


preceding page blank not filmed 


xx-17 


f RECEDING PAGE BLANK NOT FILMED 



00277 

185* 

C 


C00715 

00 

C0277 

186* 

c 

PRIORITY INTERRUPT 

000715 

o 

t /> 

00277 

1&7* 

c 


000715 

V» 

OC 301 

188* 


CClrEl-ECE 

C00739 

o 

C0302 

189* 


ECO=EO«EDE 

000737 

w 
1— » 

00 303 

19C* 


ITI (KE.GT.ECO).ANO.(INT.EO.l. llINTzO. 

C00712 

00 

00305 

191* 


1FC (KE.LT.ECn.8ND. (I NT.EO. -1 . > IINTZO . 

000760 

w 

1 

00337 

192* 



000776 

ro 

00307 

193* 


IFIKr.tr. C0)INT=1. 

000776 


0031 S 

191* 


IFIKE.G7.E1 )IN7~-1. 

001001 


00313 

19S* 


IF( («r .6 7 . ECO). ANO.IKF.lt. EC1 1>INT=0. 

001012 


00315 

19 b * 


If 1 IHPL.LE. 1 (RETURN 

001031 


C0315 

197* 

c 


001331 

£• 

00315 

198* 

c 

STATISTICS 

001031 

c 

0C315 

199* 

c 


001031 


00317 

200* 


HTrAHAXl (PE »KE 1 

CO 1 CIO 

s 

C0323 

201* 


HPC-AHAX1 (NPC,KEOI 

C01C96 

£ 

0C321 

202* 


HPD-APAX 1 (HPD.-KEO 1 

001051 

c 

00322 

203* 


SPC-SPC*TINC*P1 

001062 


00323 

231* 


SP0=SP0*IINC*P2 

001066 

> 

00323 

2C5* 

c 


001066 

O 

00321 

206 * 


IF(TI^.LT.TNAX1»»ETURN 

C01C72 


00 32b 

207* 


CCI=CCI*CC 

001101 

>%N 

OC 327 

208* 


CMl=CMI*C« 

001101 

w 

00327 

2390 

c 


CO 11 09 


00330 

210* 



C011 07 

z 

00339 

211* 


RETURN 

001107 

* 

CC331 

212* 


END 

001332 

5: 






o 






H 






*5 
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7.8A FRESNEL LENS SOLAR COLLECTOR 


Solar Insolation (ST)—**) 
Air Velocity or Wind (WD)— *4 
Ambient Temperature (TA)— 
Fluid Inlet Temperature (TFI>- 



Cell Temperature (TC) 

Fluid Flow Rate (FMD) 

Outlet Fluid Temperature (T2) 
Thermal Power Obtained (PI) 


The Fresnel lens collector model performs a thermal analysis for a con- 
centrating photovoltaic array which tracks the sun. The array may be 
cooled passively or by forced air or fluid. Fins may be used on the back to 
increase convective heat transfer to the environment. Figures 7.8A-1 and 
7.8A-2 show the physical construction of the array and the equivalent 
thermal network for the focusing collector. The purpose of the model is to 
compute the cell temperature TC, and the fluid pump rate FMD when fluid 
cooling is used. The analysis is based on a similar thermal model in 
SOLCEL [4]. 
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BASIC EQUATIONS 

1) Energy absorbed by the collector per unit area 



QH = ST*TAU*(ABC-EFF) 

where 

ST = direct beam solar insolation 

TAU = lens transmittance 

ABC = cell absorptance 

EFF = nominal cell efficiency 


2) Heat balance equations for the thermal network of 7.8A-1: 


^h ' Qt0P + %0T 

Q T0P ° H T0P (TS - T ET0P) = H l (TS-TL) 

Q B0T = H B0T (TS ' T EB0T ) * QpiN + Q FLU 
Q FIN = H FIN (TS-TA > * H,(TS-TI) 

q flu * h flu< ts - t flu> 


3) The temperature variation in the insulating bond between the cell and 
the heat sink is given by a radial conduction equation for r>a: 


.2lV’ T B 2 t „ 
^T + — -«rr B .o. 


with 9r specified at the cell radius r=a and at the equivalent lens 
radius r=b. This equation may be solved using modified Bessel func- 
tions to compute Tg at r=a given the overall heat transfer coefficient 
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and equivalent temperature of the collector minus bonding. Thus the 
cell, bonding, and collector thermal diagram reduces to 

QH*AL - 

• - ■WWW ■ • 

TC TB TE 


where 

2 

AL = lens area = n b 
TE = ^ H T0P* T ET0P +H B0T* T EB0T^ / ^ H T0P +H B0T^ 
Input Specification Notes 

Minimum input parameters to specify FO are 


CM0 

= 

Cooling mode option 

TF0 

= 

Outlet fluid temperature (CM0=2) 

NT 

= 

Number of cooling tubes (CM0=2) 

HI 

= 

Thermal conductivity/thickness of back insulation (CM0=2) 

AL 

= 

Area of lens 

NL 

= 

Number of lenses 

CL 

= 

Collector length 

CW 

= 

Collector width 

RC 

= 

Radius of solar cell 

FIR 

= 

Cooling fin/collector area ratio (CM0-O) 


The user should check inputs for consistency with those used in the 
photovoltaic model PV. For example 

? ? 

FO collector area = CL*CW >AL*NL >PV array area 

FO concentration ratio = AL/(tt*RC ) * PV concentration ratio 

? ? 

FO cell area = rr*RC >PV array area/number of cells 
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Inputs/Port 

Description 

Units 

ST 

Direct beam solar insolation 

w/m 

WO 

Air or wind velocity (default = 0.) 

m/s 

TA 

Ambient temperature 

°C 

TFI 

Inlet fluid temperature 

°C 

TF0 

Specified outlet fluid temperature 

°c 

CM0 

Cooling mode (default = 0.) 

- 


0 - natural air cooling 

1 = forced air cooling 

2 = fluid cooling 


AL 

Lens area 

m 2 

TAU 

Lens transmittance (default = 1.) 

- 

ABC 

Cell absorptance (default = .95) 

- 

EFF 

Nominal cell efficiency (default = .12) 

- 

SPA 

Lens to heatsink space (default = .025) 

m 

EL 

Emittance of lens (default = .9) 

- 

ES 

Heatsink emittance (default * .5) 

- 

El 

Emittance of the back surface (default = .5) 

- 

CW 

Collector width 

m 

CL 

Collector length 

m 

NL 

Number of lenses on collector 

- 

RC 

Radius of solar cells (default = .025) 

m 

ABL 

Absorptance of the lens (default = .05) 

- 

SPT 

Specific heat of coolant (default = 4184) 

j/kg-K 

HI 

Conductivity/thickness of the back insulation 

w/m 2 -K 


(default = 10 9 for no insulation) 
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Inputs/Port 

(cont'd) Description Units 


FIR 

NT 

MFM 

DT 

C0S 

THS 

DEN 

C0C 

HC 

CC 

CM 

C0P 


Cooling fin to flat plate area ratio 
(default = 1 for no fin) 

Number of cooling tubes 

Maximum fluid flow rate kg/s 

Diameter of cooling tubes (default = .015) m 

Conductivity of heatsink (default = 202) w/m-K 

Heatsink plate thickness (default = .003) m 

Coolant density (default = 980.) kg/m 3 

Conductivity of the coolant (default = .657) w/m 

p 

Conductivity/thickness of the cell insulator w/m-K 

Q 

(default = 10 for no insulation) 

p 

Capital cost per unit collector area per year $/m 
Maintenance cost per year $ 

Cost of operating power $/kwh 


Outputs/Port 


Description 


TC 

TS 

FMD 


T 1 

T 2 

PH 


P 1 


Cell temperature 
Heatsink temperature 
Fluid flow rate 
Inlet fluid temperature 
Outlet fluid temperature 
Collector energy absorbed 
Thermal energy collected 


Units 

°C 

°C 

kg/s 

°C 

°C 

kw 

kw 
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Outputs/Port 
(cont *d) 


FO 

Description Units 


REA 

^^.^Jeyjolds number (air cooling) 

- 

^Te'f* 

tfeyno^ number (fluid cooling) 

- 

LTI 

Last time at which the collector 
calculations were performed 

hr 

0P 

Operating Power used (state) 

kwh 


CALCULATION SEQUENCE 

RL = (AL/tt)* 5 

1) Solar Power Absorbed by the Collector 

QH = ST*TAU(ABC-EFF) 

PH = QH*AL*NL/1000. 

If QH<0.1 set TC = TA, FMD = PI = 0P = 0 and return 
If LTI = TIME and ITFI - T1K.1, return 
LTI = TIME 

2) Convert TA,TF0,TFI to °K 

3) Initial Temperature and Flow Rate Estimates 

TS = TA + QH/20 
TL = (TS + T0)*.5 
TF = (TFI + TF0)*.5 
TI = TL 

FMD = IFLU = 0 

If CM0 = 2 and TF0 >TFI, IFLU = 1 
If IFLU = 1, 


1 51 G 
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CALCULATION SEQUENCE (cont'd) 

RO = NT*SPT*(TF0-TFI)/(AL*NL) 

FMD = MIN(0.5*QH/RO,MFM) 
o Iterate 4) to 8) three times: 


Appendix 

(2)-<3) 

Ibid, (8) 


Ibid, (8) 


5) Fin Factor and HFIN Heat Transfer Coefficient 

HC = CNVC(TI,TA,WD,CL) Ibid,(2)-(3) 

HR = RADC(TI,TA,EI,1.) Ibid, (8) 

FAC = 4.318 - 4.3375*EXP(-.26795*FIR) (First pass) 

HFIN = (1/HI + 1/(HC*FAC+HR))" 1 

6) HFLU Heat Transfer Coefficient to Fluid and REF 

HFLU = 0. 


4) HTOP Heat Transfer Coefficient and TETOP 
T sky = .0552*TA 1,5 


HC1' 


REA, 


= CNVC(TL,TA,WD,CL) 


HR1 = R ADC ( TL, TSKY, EL, l.)*(TL -TSKY) 

(TL-TA) 

HI = HC1 + HR1 
TM = ,5*(TL+TS) 

HC2 = (7.25x10' 5 *TM+4.325x10" 3 )/SPA 

HR 2 = RADC(TS,TL,ES,EL) 

HL = HC2 + HR2 
HTOP = (1/HI + l/HL)** 1 
TETOP = TA+ST*(ABL+(1-TAU)*TAU*ABC)/H1 
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CALCULATION SEQUENCE (cont'd) 



If IFLU = 0 go to (7) 

HFLu\ 

= FLUC(NT,DT,CW,C0S,THS, FMD,DEN,TF,C0C) Ibid,(5)-(6) 

REF / 

7) HBOT Heat Transfer Coefficient and Temperature TEBOT 

HBOT = HFIN + HFLU 
TEBOT = ( HF I N*TA+HFLU*TF ) /HBOT 

8) Temperature and Flow Rate Updates 

H = HTOP + HBOT 
TE = (HTOP*TETOP+HBOT*TEBOT)/H 
TS = TE + QH/H 
TL = TS - HTOP* ( TS-TETQP ) /HL 
TI = TS - HFIN(TS-TA)/HI 
QFLU = HFLU(TS-TF) 

FMD = 0. 

If QFLU >0, FMD = QFLU/RO 
If QFLU >MFM*RO, 

FMD = MFM 
RA = QFLU/MFM 

TF = TF I +RA*AL*NL* . 5/ ( SPT*NT ) 

9) Check for QFLU<0 

If QFLU <0 set IFLU = 0 and repeat (4)-(8) once 

10) Cell Temperature 
ALPH = H/(C0S*TH$) 

X = SQRT(ALPH)*RC 
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CALCULATION SEQUENCE (cont'd) 



Y = SQRT(ALPH)*RL 
BETA = QH*AL/(2tt *C0S*THS*X) 

A = BETA*I1(Y)/(K1(X)*I1(Y)-K1(Y)*I1(X) ) 

B = BETA*K1(Y)/(K1(X)*I1(Y)-K1(Y)*I1(X) ) 

TB = A*K0(X)+B*I0(X)+TE 
TC = TB+QH*AL/(ff *RC 2 *HC) 
where I0,I1,K0,K1 are modified Bessel functions. 

11) Output Calculation 
T2 = 2*TF-TFI 

Convert TC,TS,T1,T2,TA,TFI,TF0 to °C 
PI = QFLU*AL*NL/1000 . 

TKP = 5.E-4*CL*CW 

1 ° 

0P=TKP+ I .0742*(CW*CL)- 2835 *WD- 567 

y .85xl0" n *FMD 2 * 855 -DT< ’ 4 * 702 ^*NT*CL 


if CM0 = 0 

if CM0 = 1 and WD>0 
if CM0 = 2 and FMD >0 


REFERENCES FOR FO 

1. J. K. Linn, "Photovoltaic System Analysis Program- SOLCEL," Sandia 
Laboratories Report SAND77-1268, 1977. 

2. E. L. Burgess and M. W. Edenburn, "One Kilowatt Photovoltaic Sub- 
system Using Fresnel Lens Concentrators," Paper 11.6, IEEE Photo- 
voltaic Specialists Conference, Baton Rouge, November 1976. 


BCS 40180-2 Rev. 


151 J 



0C101 

63* 

c 

07 -DIAMETER OF COOLING TUBES ,M, (DCF AULT =.01 5) 

050000 


00101 

64* 

c 

COS -CONDUCTIVITY OF MEAT SINK ,N/M-K , (DCF AULT- 2Q2 1 

000000 

tn 

ooici 

65* 

c 

THS -HEATSINK PLATE THICK NESS , M , CDEFAUL T= .003 I 

000030 

~r 

30101 

66* 

c 

DEN -COOLANT DENSITY, KG/M3, 

CDEF AULT =980 I 

000000 


30101 

67* 

c 

COC -CONDUCTIVITY OF THE COOLANT, V/M-K , CDEFAULT=. 6571 

oouooo 


00101 

68* 

c 

HC -CONOUCTIVITV/THICKNtSS 

OF THE Cl LL INSULATOR, 

QOC'000 


30101 

69* 

c 

V/M2-K, IDEFAULT = U**9 

FOR N: INSULATION) 

003000 


30101 

70* 

c 

CC -CAPITAL COST PER UNIT 

COLLECTOR AREA PER YEAR, 5/M2 OOODOO 


30101 

71* 

c 

CM -MAINTENANCE COST PER Y 

EAR 

oouooo 


0011 1 

72* 

c 

COP -COST OF OPERATING POWER ,11/KwH 

O^COQT 


0 3101 

73* 

c 



003000 


001 *13 

74* 


COMMON /ClhPL/IMPL .ICNT.ITEST 


003000 


00104 

75* 


COMMON /CTIME/TIME /CSIMUL/DUM l?» ,!KAX 

aooroo 


3010s 

76* 


COMMON /C0ST/CC»P,C4A,CP0 


a*»3oao 


00106 

77* 


REAL NL,NT,NFM,LTI 


OOcDOO 


0010? 

78* 


DOUBLE PRECISION MMB S 1 0, MHBSI 1 ,MMd SKO,MHBSK l 

Q9QD30 


001 10 

79* 


DOUBLE PRECISION X,T 


oooroo 


L01 1 u 

80* 

c 



003D00 


301 tO 

61* 

c 

INITIALIZATION 


003 n 00 


aono 

8 2* 

c 



ooeoca 


00111 

63* 


IFCIMPL.GT.01G0 TO lDO 


ooooco 


33113 

84* 


IFIWD.EQ..99\>9 3) WD = ?, 


000002 


33116 

6 5* 


IFICM0.E0..9'3999)CM0 = 0« 


303006 


00117 

86* 


IF! At«£0..99i99) AL=.u9 


C00012 


QC121 

87* 


1FCTAU.E0 .. 99999 |TAU=1 . 


O' 1 301 7 


0C123 

88* 


1FI ABC .EC, .99999) ABCr.eS 

o o 

, 0000.24 


30126 

69* 


IFlFFF.EC..99999)EFF=.12 

“ 2 

000031 


30127 

90* 


IF ISP A. EC.. 999 99 ISP A = .325 

Tj O 

OOCOJ6 


00131 

91* 


IF CEL .E0..99999)EL=.9 

3 § 

000043 


30133 

92* 


IFt£S»£C..99999ILS=«5 

o f 

0OC050 


30136 

93* 


IF IE1 .E0..99999)EI= .5 

» /= 

OOOPSS 


•301 3 7 

94* 


IF IRC. EO. .999 99 )RC = .025 


000062 


30141 

95* 


1FIABL .EC.. 99999 )ABL=.OS 

g 

c > 

003067 


00143 

96* 


IFISPT.EC..99999)SPT=4184 

a* & 

000074 


C014 5 

97* 


1FIMI.E0. .99999)hl=! .E9 

£T m 

003131 


00147 

98* 


1FIF1R.EC..999V9)FIP=1. 

5 __ 

motes 


30151 

99* ■ 


IFIDTcEC. .999 99 >DT = .^15 

< </) 

U00113 


00153 

ICO* 


IF'.COS . EC.. 99999 ICO S = 202 


000120 

m 

3 Cl 56 

lcl* 


IFITHS.EC..99999ITHS=.003 


000125 

w 

O 

00157 

102* 


1FI0EN.EC..99999IDEN=983 


0 JO 132 

oo 

00161 

103* 


1FI COC .EC. .999 99 I COC = .6 57 


003137 

-P» 

00163 

lu 4* 


1FI HC.EO.. 99999 IHC=1 ,E9 


006144 

o 

LC165 

1C 5* 


RL=SCRTIAL/3.l4la926) 


000151 

CD 

00166 

lufc* 


F AC=4 .3 16 -4 ,3375*EXPI -«2679S*FIR) 


073160 

0 

1 

00167 

lu7* 


TMAX1=TMAX*. 99999 


00J171 

INJ 

03167 

106* 

c 



000171 


L'Q17c 

lu9* 


100 CONTINUE 

. 

000175 

fD 

33173 

11C* 

c 



000175 

< 

3017C 

111* 

c 

SOLAR POWER ABSORBED BY THE COLLECTOR 

000175 


0017C 

112* 

c 



000175 


00171 

113* 


CH=ST*T»U*C»BC-EFFI 


000175 


30172 

114* 


PH=CH*AL*NL/1CD0. 


OOC.702 


00173 

1 15* 


I F 1 OH jGT. C.ll CO TO 201 


•3 00 20 6 


0 0175 

116* 


TS = TA 


□03012 


00176 

117* 


TC=TA 


300214 


00177 

118* 


0P3=3. 


000215 


002 00 

119* 


FKD=0. 


006216 
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00201 

120 



P1=0. 


01C21T 

00202 

121* 



60 TO 920 


000220 

00203 

122* 


201 

IFHLTI.EQ. TIME). AND. tABS<TFI-Tl».LT. ,11160 TO 920 


000222 

00235 

123* 



lti=tihe 


030243 

0020b 

124* 

C 




030241 

□ 0235 

125* 

c 


CONVERT TA , 7F0 , TF I FROM CENTIGRADE TO KELVIN 


010240 

03235 

126* 

c 




010243 

0C216 

127* 



TA=TA*273 


010242 

00237 

128* 



TF0=TF0*27J 


010245 

602 III 

1*9* 



TF1-IF 1*273 


01-250 

032 I U 

130* 

c 




010250 

00210 

131* 

c 


INITIAL TEMPERATURE AND FLOW RATE ESTIMATES 


010250 

002 1 U 

132* 

c 




C1U250 

302 1 1 

133* 



TS=TA+GH/20. 


01G253 

30212 

134* 



TL=»TS*TA1*.5 


010257 

00213 

135* 



TF=tTFl*TF01*.5 


DOC 262 

0021<« 

136* 



TI-TL 


01C265 

00215 

137* 



IFLU=0. 


013266 

00215 

138* 



F MO-O • 


010267 

00217 

139* 



IFI < AB: 100-2. 1 .LT..1) .AND.«TF0.GT.TFI»UFLU=1 


030273 

00221 

140* 



1FIIFLU.NE.IIC0 TO 301 


300712 

GC223 

141* 



fcO = NT*SPT*< TFO-TFI »/CAt*NL 1 


01-715 

00224 

142* 



FHD-hFN 


QIC 326 

00225 

143* 



IF(RO.OT.C.IFMD=AMIN1I .5*OH/RO,MFNI 


01C733 

03227 

144* 


301 

CONTINUE 


C13344 

U0227 

145* 

c 




010744 

0 3*27 

146* 

c 


ITERATE HEAT COEFFICIENT CALCULATION THREE TINES 


01-744 

□ C227 

147* 

c 




010344 

0023 J 

146* 



LO0P = 9 


0 ?, a744 

00231 

149* 


400 

CONTINUE 


0*0745 

032 7 2 

1SG* 

c 




030345 

00231 

151* 

c 

HTOP t HEAT TRANSFER COEFFICIENT, AND TETOP.TOP EQUIVALENT 

TEMPERATURE 

013745 

00232 

152* 



TSKY-.G552*ITA**1.51 


010345 

002 33 

153* 



CALL CNVCIHC1,REA,TL,TA,WD,CL» 


010352 

CC23H 

154* 



CALL RADC«HP1,TL,TSKY.EL, 1.) 


020362 

□ C235 

155* 



hri~hri*<tl-tskvi/itl“Tai 


010771 

00235 

' 156* • 



H1=HC1.HK 1 


010431 

03237 

157* 



TP:.5* ITL*TS1 


J31J422 

002*40 

158* 



h C2 “17.25*1 •E-5*TN.4, 32 5E -31/ SPA 


010435 

002*41 

159* 



CALL RADC<hR2,TS,Tl,ES,EL » 


013412 

002 42 

16C* 



HL-HC2 *HR 2 


010421 

00243 

161* 



HTOPAl./C 1./H1M./HL1 


Q1C424 

03244 

162* 



TETGP:TA«ST*< A6L«< 1-TAU)*TAU*ABC)/H1 


31C434 

00244 

163* 

c 




310434 

03244 

164* 

c 


HEAT TRANSFER COEFFICIENT HF IN 


010*34 

00245 

165* 



CALL CNVCIhC2,RE J TI,TA,w0,CL» 


-*C445 

00245 

166* 



CALL RACCIhP,TI,TA,El,l.) 


010455 

0C247 

167* 


„ 

HFIN = 1 ./I 1 ,/Hl*l ./<HC2*FAC*HR» > 


01C464 

0C24 7 

166* 

c 




C10464 

00247 

169* 

c 


FLUID HEAT TRANSFER COEFFICIENT HFLU AND REYNOLDS NUMBER 

REF 

010464 

OG25C 

170* 



HFLU=0. 


01C477 

002 51 

171* 



IFlIFLU.E0.2lG0 TO 7u0 


CIO 500 

0C253 

172* 



CALL flucjhflu,ref,nt,dt,cw,cos,ths,fmd,oen,tf,cocj 


030502 

00253 

173* 

c 




010532 

00253 

174* 

c 


equivalent eottom temperature ano heat transftr coefficient 

0*0532 

00254 

175* 


700 

CONTINUE 


01-520 

00255 

176* 



HBOT-HFIN.HFLU 


010520 
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00256 

177* 



TEBOT=tHFIN*TA*HFLU*TFSXHBOT 

C30522 

00256 

17S* 

C 



070522 

0025 6 

179* 

C 


UPDATE TEMPERATURE ANO FLOW RATE 

070522 

00257 

180* 



H=HTOP*«BOT 

3 n 0531 

00260 

lei* 



TEilhTOP*TETQP*HBOT*TEBOT»/M 

onesso 

00261 

182* 



T5-TE ♦C'H/M 

07&S*1 

00262 

18 3* 



TL=TS-HTOP*ITS-T£TOP>/HL 

07u5*5 

GC263 

164* 



YI=TS-hr 1N*<TS-TA)/HI 

00C552 

00264 

165* 



lflu=hflu*«ts-tf» 

0*10560 

00264 

166* 

c 


6R1TEI6 ( 108)HFIN ( HB0T«TEB0? tHTC:*, TETOPfH »TEtTS»TL»TI «OFLU«RO 

0*70560 

00264 

167* 

£4 

108 

FORMAT 1 1H t *FOC * f 8E10.2,/,5X,8EI3.2» 


0026s 

168* 



F**D = -w. 

070564 

00266 

169* 



IFiQFLU.LC.I'.IGO TO 800 

004555 

002 70 

190* 



1F(GFLU.GT.(MFH*R0> )G0 TO 799 

370*67 

00272 

191* 



F^D-OFtU/RO 

07057* 

30271 

192* 



60 TO 60u 

070177 

OC274 

193* 


799 

FFO=HfM 

036601 

00275 

194* 



hA=CFLU/«FH 

000602 

00276 

195* 



TF=7FI*RA*AL*NL*.5/ISPT*NT> 

0766-35 

00277 

196* 


800 

CONTINUE 

3*10616 

00277 

197* 

c 



070616 

001 DC 

198* 



L0QP=LG0P^1 

03«tl6 

00111 

199* 



IF< LOOP *LE» 2160 TO 400 

me 4 ?o 

001-31 

i'ut* 

c 



070*20 

0C301 

201* 

c 


CHECK FOR EFFECTIVE FLUIO C00LIN6 

0"0623 

00131 

2c>2* 

c 



070620 

03133 

2u3* 



IFICFLU.GE.O.IGO TO 900 

07U623 

0 0105 

20** 



1FLU=0. 1 

036626 

uOlUe 

2C5* 



60 TO 470 

076627 

00337 

206* 


900 

CONTINUE 

37*2631 

00307 

2C7* 

c 



00C6J1 

30307 

208* 

c 


CELL TEMPERATURE 

000631 

0C137 

2u9* 

c 


' 

000631 

00310 

210* 



ALPH=h/<COS*TKS» 

370631 

003 U 

2U* 



X=S8RTCALPH)*RC 

UCC635 

0 3312 

212* 



Y=SCRTULPhJ*RL 

070645 

uCl 1 3 

213* - 



6E7A=0H*AL/<2.*3.141S9*C0S*ThS*X1 

070650 

00114 

214* 



8I1Y=mkPSIH1,T,IER) 

070663 

00315 

215* 



fcKAXXKMPSK l(lfXtlER) 

00C672 

00316 

216* 



tIiX=MHPSIia,X f IER» 

COuTOl 

00317 

217* 



bMV;MHPSKHl f Y,IERI 

070710 

02320 

218* 



bK jXzPf'BSK. Il *X, IER) 

070717 

0C321 

219* 



feIJX:p^t>SIJIl.X,l£R» 

376726 

03322 

22C* 



A=dETt*PllY/IBHlX*e[ 1Y-BK 1 Y*B I IX > 

07C73S 

0 0323 

22 l* 



b=8FTA*Bh J Y/(PK!X*0HY~BKIY*B1 1X1 

070746 

00324 

222* 



IB=/i*Pnnx-*6*Bl6X*TE 

070752 

00325 

223* 



TC=TB*L'H*AL/<3.141S9*RC*RC*HCI 

070757 

00325 

224* 

c 



071.757 

00125 

225* 

c 


OUTPUT CALCULATION 

070757 

CO 325 

226* 

c 



070757 

00326 

227* 



TC=TC-273 

07U770 

00327 

228* 



TS=TS-273 

070772 

0033 0 

229* 



Tl=TFI-273 

030775 

00331 

23P* 



72=2 »*TF -IF 1-273 

071?30 

00132 

231* 



T A = T A-273 

mi.-aS 

00331 

232* 



TFI=TFI-273 

031-713 

03334 

233* 



TFO = TF 0-273 

071311 


i 
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00335 

234* 


P120FLU*AL*NL/1000, 

001014 

00336 

23S* 


*n=o. 

0H1"?1 

00337 

236* 


lF1ABSICHO-l.).LE..l>Rei:.O742*|ICW*CL)**.2035)*WO**.S67 

001022 

00341 

237* 


IFlFMO.LE.O.IGO TO 909 

001047 

00343 

238* 


lF(CH0.GT.l.l»REl=7.65£-lI*CFHD**2.*551*«DI**f-4.7C2l»*NT*Cl 

otiose 

00345 

235* 

909 

CONTINUE 

C"io75 

003«o 

24C* 


TKP=5.E-4*CL*CW 

G01075 

00347 

241* 


1F»I0P.NE.uI0P0=TKP*RE1 

oonoo 

00351 

242* 

92C 

1F<T1*E.LT,TMAX1)RETURN 

001106 

00353 

243* 


IFUPPl .ir«?)RCT0f»* 

C^l 114 

coiss 

244* 


CCAF=CCAP*CC*»L*NL 

001121 

0035b 

245* 


CPArCHA*CH 

071130 

00357 

246* 


CP0=CPC*C0P*0P 

011133 

(J 036u 

247* 


RET URN 

001137 

0036 1 

248* 


two 

001515 


n 

O 


7.8B FLAT PLATE SOLAR COLLECTOR 


Solar Insolation (ST)L_-J 

Collector Tilt (TLT) 

Fluid Inlet Temp. (TFI 
Air Velocity or Wind (WD 
Ambient Temperature (TA). 



Cell Temperature (TC) 

Fluid Flow Rate (FMD) 

Outlet Fluid Temperature (T2) 
Fluid Cooling Power (PI) 


The flat plate component performs a thermal analysis on a nonconcentrating 
photovoltaic array. Three types of cooling may be used: 

• Front surface cooling using natural or forced air. 

• Back surface cooling using natural or forced air with or without 
a finned back surface. 

• Fluid cooling using tubes on the back and N glass covers (N = 

0 , 1 , 2 , 3 ). 

Figures 7.8B-1 and 7.8B-2 show the physical construction of the array and 
the equivalent thermal network for the flat plate component. The purpose 
of the analysis is to compute the cell temperature TC and the fluid pump 
rate FMD when fluid cooling is used. The analysis is based on the flat 
plate thermal model in SOLCEL [4], except that an empirical equation due to 
Klein is used to compute the top loss coefficient for 1 to 3 glass covers. 


152 


BCS 40180-2 Rev. 




FP 


TEMPERATURES 



T AMB 


.GLASS COVER(S)(OPTIONAL) 

■AIR SPACE 
•SOLAR CELL 

-INSULATION 
•MOUNTING PLATE 
■INSULATION (OPTIONAL) 
.FLUID COOLING (OPTIONAL) 


Figure 7.8B-1 Physical Diagram of Flat Plate Collector 



Figure 7.8B-2 Equivalent Thermal Network for Flat Plate Collector 
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BASIC EQUATIONS 


z t rr " ~ — - - - 

for the network of Figure 7.8B-2. 

Q H - st * t n< ab - eff) = q TOp + Q 


q T0P * H T0p( T CEU - W 

q bot ■ h bot (t ceu - r EB0T ) - H C (r CELL . t plate ) 
= q fin + q flu 

Q FIN ■ H FIN< W - T AM e ) - Hj(T pute . t insul ) 
Q FLU - FMD*P(TF0 - TFI) = H clI1 (T n r , 

flu u plate 'fluid' 


Where H T0P» h B0T» h C-“ denote transfer coefficients, and 
T N = trar, smittance of the N-covers 
AB - collector cell absorptance 
EFF = nominal cell efficiency 


T ebot . equivalent both- tea*. (. TflHe «,th „o fl u1d cooling) 

T ' = f ' U ’ d SP6CifiC heat/unit <*" area * No . of cooling tubes 
FLUID - average fluid temperature = (TF0 + TFI ) /2. 


Input Sp ecification Notes 

Minimum input parameters to specify FO are 


CM0 

TF0 

NG 

HI 

CW 

CL 

NT 

FIR 


Cooling mode option 

Number o'f^assTo^s^ (CM0=2) 

Col?e U ct t or i ifdt t h hiCkneSS ° f the back '"Elation 
Collector length 
Number of cooling tubes (CM0=2) 

Cooling fin/col lector area ratio (CM0=O) 


The user should check the consistency of these inout<; (o 

area) with those of the trackinn r P * {e * 9 *’ collect or 

component PV. 9 C ° mP ° nent 50 and the Photovoltaic 
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Inputs/Port 

Description 

Units 

ST 

Global solar insolation 

2 

w/m 

TLT 

Collector tilt 

Deg 

WD 

Air or wind velocity (default = 0.) 

m/s 

TA 

Ambient drybulb temperature 

°C 

TFI 

Inlet fluid temperature 

°C 

TF0 

Specified outlet fluid temperature 

°c 

MFM 

Maximum fluid flow rate 

kg/s 

RE 

Tracking power request 

kw 

CM0 

Cooling mode (default = 0.) 

0 = natural air cooling 

1 = forced air cooling 

2 = fluid cooling 


NG 

Number of glass covers (default = 0.) 

- 

TN 

Transmittance of the N-covers 

- 

AB 

Collector cell absorptance (default = .9) 

- 

EFF 

Nominal cell efficiency (default = .12) 

- 

EC 

Emittance of cell (default = 0.5) 

- 

EG 

Emittance of the glass covers (default = .9) 

- 

EP 

Emittance of the back surface (default = .9) 

- 

CW 

Collector width 

m 

CL 

Collector length 

m 

SPT 

Specific heat of coolant (default = 4184.) 

j/kg-K 

HI 

Conductivity/thickness of the back insulation 
(default = 10 9 for no insulation) 

w/m 2 K 


152C 


BCS 40180-2 Rev. 


Inputs/Port 

(cont'd) 

Description 

Units 

FIR 

Cooling fin to flat plate area ratio 
(default = 1. for no fin) 

- 

NT 

Number of cooling tubes (default = 1) 

- 

or 

Diameter of cooling tubes (default = .015) 

m 

C0P 

Conductivity of mounting plate (default = 202. 

) w/m-K 

THP 

Mounting plate thickness (default = .003) 

m 

DEN 

Coolant density (default = 980.) 

kg/m 3 

C0C 

Conductivity of the coolant (default = .657) 

w/m-K 

HC 

Conductivity/thickness for cell insulation 
(default = 10^ for no insulation) 

w/m^-K 

CC 

Capital cost per unit area per year 

%/m 2 

CM 

Maintenance cost per year 

$ 

CP0 

Cost of operating power 

$/kwh 

Outputs/Port 

Description 

■in its 

TC 

Cell temperature 

°C 

TP 

Mounting plate temperature 

°C 

FMD 

Fluid flow rate 

kg/s 

T1 

Inlet fluid temperature 

°C 

T2 

Outlet fluid temperature 

°c 

PH 

Collector energy absorbed 

kw 

PI 

Thermal energy collected 

kw 

0P 

Operating power used (state) 

kwh 

REA 

Reynolds number (air cooling) 

- 

REF 

Reynolds number (fluid cooling) 

- 

LTI 

Last time at which the flat plate array 
calculations were performed (used internally) 

hr 

152D 

BCS 

40180-2 Rev. 



CALCULATION SEQUENCE 

1) Solar power absorbed by the collector 

QH = ST*TN*{AB - EFF) 

PH = QH*CL*CW/1000 

If QH<0.1 set TC = TA, FMD = PI = OP = 0 and return 
If LTI = TIME and |TFI - Tl| < .1, return 
LTI = TIME 

2) Convert TA, TF0, TFI to °K 

3) Initial temperature and flow rate estimates 

TC = TA + QH/20 
TI = (TC + T A ) * . 5 
TF = (TFI+TF0)*.5 
TP = TI 
FMD = 0 
IFLU = 0 

If CM0 = 2 and TF0 >TFI, IFLU = 1 
If IFLU = 1, 

RO = NT*SPT*(TF0 - TFI)/CW*CL 
FMD = MIN(MFM,0.8*QH/R0) 
o Iterate 4) to 8) three times: 

4) HTOP heat transfer coefficient and REA 
TSKY = .0552*TA 1,5 
HCl\ 

Y CNVC(TC, TA, WD, CL) See (2)-{3) in Appendix 

REA/ 

If NG = 0, 
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CALCULATIONS (cont'd) * a 

HR1 = RADC(TC,TSKY,EC,1. )*(TC-TSKY) Ibid, (8) 

TCTAJ ‘ 

HTOP = HC1 + HR1 
If NG > 0, 

HTOP = HTGLAS(N,TA,TC,HC1,EC,EG,TLT) Ibid, (7) 

5) Fin factor FAC and HFIN heat transfer coefficient 

HC2 = CNVC(TI, TA, WD, CL) Ibid, (3) 

HR 2 = RADC(TI, TA, EP, I.) Ibid, (8) 

FAC = 4.318 - 4.3375*exp(- .26795*FIR) (first pass) 


HFIN = (1/HI + 1/ (HC2*FAC + HR2))" 1 

6) HFLU heat transfer coefficient to fluid and REF 

HFLU = 0. 

If IFLU = 0 go to 7) 

HFLU\ 

= FLUC ( NT, DT, DW, C0P , THP , FMD, DEN , TF , Cf)C ) Ibid, (5)-(6) 

REF/ 

7) HBOT heat transfer coefficient and equivalent temperature TEBOT 

HBOT = (1/HC + 1/ ( HFIN + HFLU))’ 1 
TEBOT = (HFIN*TA + HFLU*TF)/(HFIN + HFLU) 

8) Temperature and flow rate updates 

TC = (QH + HTOP*TA + HBOT *TEBOT ) / ( HTOP + HBOT) 

TP = TC - HBOT*(TC-TEBOT)/HC 
TI = TP - HFIN*(TP - TA)/HI 
QFLU = HFLU* (TP - TF) 
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CALCULATIONS (cont'd) 

FMD =/u. if QFLUSO. 

^QFLU/RO if QFLU >0. 

If QFLU > MFM*RO, 

FMD = MFM 
RA = QFLU/MFM 

TF = TFI + RA*CL*CW*.5/(SPT*NT) 

9) Check for QFLU < 0 

If QFLU <0 set IFLU = 0 and repeat 4) to 8) once 

10) Output calculations 

T2 = 2*TF - TFI 

Convert TC,TP,T1,T2,TA,TFI,TF0 to °C 
PI = QFLU*CL*CW/1000 
If CM0 = 0 
0P = RE 

If CM0 = 1 and WD>0 

0P = RE+.0742*(CW*CL)- 2835 *WD‘ 567 
If CM0 = 2 and FMD>0 

0P = RE + 7.85 x 10" li *FMD 2,855 *DT^“ 4,702 ^*NT*CL 
REFERENCES FOR FP 

1. S. A. Klein, M.S. Thesis, "The Effects of Thermal Capacitance Upon the 
Performance of Flat Plate Solar Collectors," University of Wisconsin, 
1973. 

2. J. A. Duffie and W. A. Beckman, Solar Energy Thermal Processes. 
Wiley, 1974. 

3. F. Kreith, Principles of Heat Transfer , 3rd Edition, International 
Textbook Company, 1973. 
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SUBROUTINE FP 


ENTRY POINT 000741 


STORAGE USEO: C0DE(1) 001270; DATA (01 OoOlll; BLANK C0NN0NC2) 000000 
COMMON BLOCKS! 

0033 CIMPL 0C0001 

0004 CTIME 000001 

0005 CSIMUL 000310 
uOC6 COST Q00003 


EXTERNAL REFERENCES (BLOCK. NAME) 


LOOT 

CNVC 

0010 

RADC 

0011 

HTGLAS 

0012 

FLUC 

00.13 

EXP 

D 014 

XPRR 

0015 

NERR3S 


oo 

c~> 

oo 

-p» 

o 

00 

0 

1 

ro 

XJ 

n> 

< 


STORAGE ASSIGNMENT (BLOCK. TYPE. RELATIVE LOCATION. NAME) 


0031 


0wQi<46 

100L 

0001 


000173. 

2P1L 

0001 


000315 

3C1L 

U001 


000316 

4Q0L 

0001 


000360 

401L 

0031 


CCD j72 

4C2L 

0001 


000494 

7PCL 

0031 


000530 

799L 

0001 


000545 

BPCL 

0001 


000560 

90 OL 

3001 


0u06fc7 

9C9L 

C0C1 


000674 

92DL 

0036 

R 

corooo 

CCAP 

0006 

R 

000301 

CMA 

0006 

R 

030002 

CPOS 

liOQ5 


000003 

OUM 

DC CO 

R 

000001 

FAC 

COCO 

R 

0C0D21 

HBOT 

uooo 

R 

000311 

HC1 

0000 

R 

000014 

HC2 

0000 

R 

Cl.001 7 

HFIN 

COOC 

R 

00C020 

HFLU 

COCO 

R 

000012 

HR1 

COGO 

R 

000016 

HR2 

0011 

R 

OOQOOO 

HiGLAS 

0003 

R 

OCCOll 

HTOP 

CiOOC 

I 

"00005 

IFLU 

0033 

I 

000000 

IMPL 

0000 


00006 3 

INJPS 

COOP 

I 

000007 

LOOP 

0C3D 

R 

00CD23 

QFLU 

COCO 

R 

000002 

CH 

0000 

R 

C00024 

RA 

3000 

R 

000315 

REM 

coco 

R 

000325 

RE 1 

0033 

R 

CCOCQo 

RO 

0031* 

R 

3C0C22 

TEbOT 

COCO 

R 

000004 

TF 

cooo 

R 

000303 

TI 

C004 

R 

0C3QD3 

TIME 

0035 

R 

003007 

TMAX 

DODO 

R 

030033 

TMAX 1 

OOOC 

R 

C0001C 

TSKY 










acioo 

1* 

CFP 



OOOOCO 

00101 

2* 



SUBROUTINE FP<TC,TP,FMD,T1,T2,PH,P1 , OP.OPO .IOP.REA .REF.LTI. 

030000 

00101 

3 * 


1 

st.tlt.wd.ta.tfi.tfo.mfm.re.cmo.ng.tn.ab. 

030000 

00101 

4* 


2 

eff.ec.eg.ep.cw.cl.spt.hi.fir.nt.dt.cop.thp.den. 

oocooo 

UG101 

5* 


3 

COC.HC.CC.CM.CPO) 

000000 

00131 

6* 

C 



000000 

00131 

7* 

C 


PURPOSE THIS COMPONENT PERFORMS A THERMAL ANALYSIS 

ooorco 

UC131 

8* 

C 


ON A NONCONCENTRATING PHOTOVOLTAIC ARRAY. 

oocooo 

00101 

9* 

C 


THPEE TYPES OF COOLING MAY BE USEDC 

000000 

00101 

10* 

C 


FRONT SUPFACE COOLING USING NATURAL OR FORCED AIR 

03D n 09 

03131 

11* 

c 


BACK SURFACE COOLING USING NATUR4L GR FORCED AIR 

030330 

UC131 

12* 

c 


WITH OR WITHOUT FINE. 

C"OPDO 

03131 

13* 

c 


FLUID COOLING USING TUBES ON THE BACK AND NG 

COD3UO 

0C131 

14* 

c 


GLASS COVERS ( NG=0 , 1 ,2 » 3 ) • 

D30G00 



!■ ' *> 

r 
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00101 

15* 

C 




000000 

C0101 

16* 

c 

WRITTEN 

BY Y.K.CHAN, 11-6-78* VERSION 1 


OOCPOO 

00101 

17* 

c 




oocooo 

CC101 

16* 

c 




ononoo 

00101 

19* 

c 

METHOD 

BASED ON THE FLAT PLATE THERMAL MODEL IN SOLCEL 

• 

OH^OO 

00101 

20* 

c 


EXCEPT THAT AN EMPIRICAL EQUATION DUE TO KLEIN IS USED 

030030 

00101 

21* 

c 


TO COMPUTE THE TOP LOSS COEFFICIENT FOR 1 TO 

l 

□ "•CHOD 

00101 

22* 

c 


GLASS COVERS 


oouooa 

00131 

23* 

c 

calling 

SEOUENCE 


Doo^on 

00101 

21* 

c 

OUTPUTS 


C 30000 

00101 

25* 

c 

ST 

-GLOBAL SOLAR INSOLATION, W/M2 


DOO-’OO 

ocitu 

26* 

c 

TC 

-CELL TEMPERATURE ,C 


coeooo 

00101 

27* 

c 

F MO 

-FLUID FLOW RATE, KG/S 


030000 

00131 

28* 

c 

T 1 

-INLET FLUID TEMPERA TURE.C 


O’lCODO 

0 0131 

29* 

c 

T 2 

-OUTLET FLUID TFMPER ATURE ,C 


DOCOOO 

00101 

30* 

c 

PH 

-COLLECTOP ENERGY ABSORBED, KW 


DOODCO 

00131 

31* 

c 

PI 

-THERMAL ENERGY COLLECTED, KW 


013300 

00131 

32* 

c 

OP 

-OPERATING POWER USED ISTA TE ) ,K WH 


030000 

00101 

33* 

c 

REA 

-REYNOLDS NUMfeEPtAIP COOLING) 


oonooo 

00101 

31* 

c 

RLF 

-PEYNOLOS NUMBERiFLUID COOLING) 


Q1CC30 

00101 

35* 

c 

LTI 

-LAST TIME AT WHICH THE FLAT PLATE ARRAY 


03C000 

ocioi 

36* 

c 


CALCULATIONS WERE PERFORMED (USED INTERNALLY) 

033000 

□ 0101 

37* 

c 

INPUTS 


DDCUO 

0010a 

3 8* 

c 

TLT 

-COLLECTOR TILT , DEGREES 

5 O 

COuDOO 

ccim 

39* 

c 

WO 

-AIR OR WIND VELOCITY.M/S, (DEFAULTED.) 


03D330 

00101 

1C* 

c 

TA 

-AMBIENT PRY3ULB TEMPERATURE, C 

S) r- 

D3D0D0 

□ 0131 

11* 

c 

TFI 

-SPECIFIED INLET FLUIO TEMPER A TURE , C 

03G300 

00101 

12* 

c 

tfo 

-SPECIFIED OUTLET FLUID TEMPER ATURE ,C 


C f ’uPQO 

00131 

13* 

c 

MFM 

-MAXIMUM FLUID FLOW RATE, KG/S 

<r > 
5s» o 

c nr 

O3JO0O 

ocm 

14* 

c 

RE 

-TRACKING POWER KECUEST.KW 

033000 

00131 

15* 

c 

CMO 

-COOLING MODE (DEFAULTED) 

03C330 

00101 

16* 

c 


P-NATURAL AIR COOLING 


G3DOCO 

00131 

<)7* 

c 


1=FCRCED AIR COOLING 

000000 

00101 

H8 * 

c 


2 -FLU ID COOLING 


003030 

00101 

19* 

c 

NG 

-NUMBER OF GLASS C0VERS(0EFAULT=0) 


oocooo 

00101 

50* 

c 

TN 

-TRANSMITTANCE OF THE NG GLASS COVERS 


03C000 

00101 

51* - 

c 

AB 

-COLLECTOP CELL ABSORPTANCE (DEFAULT: ,9) 


oocooa 

□ 0131 

52* 

c 

EFF 

-NOMINAL CELL E FF IC IENC Y ( DEFAULT: . 1 ? ) 


036000 

00101 

53* 

c 

EC 

-EKITTANCE OF CELLI0EFAUlT=.5 ) 


030000 

00101 

51* 

c 

EG 

-EMITTANCE OF GLASS COVE«S(DEF AL'LT=,9) 


030000 

00131 

55* 

c 

EP 

-gmittancf of the back SURFACE IDEFAJLT:,9> 


030000 

□ 0131 

56* 

c 

CU 

“COLLECTOP WIDTH, M 


oocnoo 

001 'll 

57* 

c 

CL 

-COLLECTOR LENGTH, M 


0 30000 

00101 

58* 

c 

SPT 

-SPECIFIC HEAT OF COOLANT ,J/KG -K , ( DEFAULTER 181 ) 


030000 

00101 

59* 

c 

HI 

-CONDUCTIVITV/THICKNESS OF THE BACK INSULATION, W/M2-K, 

oocneo 

00101 

60* 

c 


(DEFAULT:!. EV FOR NO INSULATION) 


030000 

00131 

61* 

c 

FIR . 

-COOLING FIN TO FLAT PLATE AREA RATI0(DEFAULT=1, 

FOR NO FIN) 

000^00 

□ 0101 

62* 

c 

NT 

-NUMBER OF COOLING T UBE S ( DEFAULT; 1 ) 


G OC303 

00131 

63* 

c 

DT 

'DIAMETER OF COOLING TUBE S , M, ( D EF AULT=. 31 5 ) 


330000 

30101 

61* 

c 

COP 

“CONDUCTIVITY OF MOUNTING PL A T E , W /M -K , 4 DEF AULT=2Q2) 

□on ooo 

00101 

65* 

c 

THP 

- M OUNTING PLATE THICKNESS ,M,( DEFAULT:, a r '3> 


OOliOOO 

JC1Q1 

66* 

c 

DLN 

-COOLANT DENSITY, KG/M3, (0EFAULT = 983) 


OOOOuO 

00101 

67* 

c 

COC 

-CONDUCTIVITY OF COO L AN T , W/M-K , ( DEF AUL T :,6&7 


□ounoo 

□ 0101 

6 B* 

c 

HC 

“CONDUCTIVITY/THICKNESS FOP. CELL INSULATION ,W/M 2“ 

K, 

030OC0 

00101 

69* 

c 


(DEFAULT=1.E9 FOP NO INSULATION) 


030000 

00131 

70* 

c 

CC 

“CAPITAL 30ST PER UNIT ARFA PER VFAR,*/M2 


DOUiCD 

ocm 

71* 

c 

CM 

-MAINTENANCE COST PfR YEAR,* 


030000 
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00101 

- 72* 

C CPO -COST OF OPERATING POWER* S/KMh 


000000 

00101 

73* 

C 


oogooo 

0C103 

74* 

COMMON /CIMPL/IMPL 


010100 

0010(i 

“ 75* 

COMMON /CTImE/TIME /CSIMUL/DUM <71 ,TMAX 


OIOOOO 

00105 

76* 

COMMON /C0ST/CCAP,CMA,CPOS 


000000 

UC1C6 

77* 

REAL L TI* MFM* NG.hl 


o^cnoo 

0 0106 

78* 

C 


0101C0 

00106 

79* 

C INITIALIZATION 


000000 

00106 

, 80* 

C 


0^0000 

00137 

81* 

lFClMPL.bT.t 1 GO TO 1U0 


Q1LP00 

00*111 

82* 

1FIW0.EC. ,99999 )t.D=3. 


000002 

uniij 

83* 

IFtCMO.EO ..99999 ICMO^O. 


OUOPOfe 

0 0 1 1 5 

84* 

IF C NG.EO.. 99999 )NG=0. 


OOUP12 

00117 

85* 

1F< Ab. EC. .999991 AB = .9 


000016 

□ 0121 

86* 

lF<EFF,Eti..999 99!E r Fi.l2 


000023 

03123 

8 7,* 

IFiEC.EC..99999)tC=.5 


000030 

0012a 

8 8*. 

IF t EG .EO . • 99999 ) EGS .9 


000035 

00127 

8 9* 

lF(EP#tO«»99999)CP-»9 


000042 

□ 0131 

90* 

1FISPT.EC..99999»SPT=41E(| 


000047 

00123 

91* 

IFJH1 .EC.. 999991 l(I = l.E9 


000^54 

GC135 

92* 

IF<FIR.E{f..99999)FIR=l* 


000061 

00137 

93* 

IF ( NT . EC ■ .99999 ) NT- 1 


0 ■’ 0 6 6 

3 C 1 4 1 

94* 

IFlDT.EO.. 99999) 0T=. CIS 


01C173 

0 214 3 

95* 

1FCCCP.EC ..99999 »C0P=2a2. 


D10I00 

□0145 

96* 

IFITHP.E0..99999)THP=.Q03 


□10135 

□ 014 7 

97* 

IF I DEN .EC . .99 999 )DEN=9 8C 


, 0H0112 

00151 

98* 

iriCOC. EC.. 99999 )COC=. 657 


000117 

0015 3 

99* 

1F1PC.E0.. 99999 )hC=l.E9 

4 - 

000124 

UC1SS 

luO* 

TMAXUTMAX*. 99999 


0001 3 1 

□ CIS b 

lui* 

FAC=4.21E-i.3375*EXP«-.26795*FIR» 


000134 

00157 

132* 

IOC CONTINUE 


O^O 146 

00157 

103* 

C 


□00146 

0015 7 

104* 

C SOLAP POWEP ADSORBED BY COLLECTOR 


003146 

00157 

105* 

C 


010146 

00160 

S U 6* 

<IHxST*TN*Ufj-EFF) 


00U146 

□ 0161 

lu7* 

PH=Ch*CL*CW/1001. 


0^0 153 

□ 0162 

108* 

1 F < CM .GT. 0 • 1 >. GO TO 2dl 


000157 

□ 0161 

109* 

1P=TA 


000163 

□ 0165 

110* 

CP0=0. 


000165 

00166 

111* 

TC=TA 


070166 

00167 

112* 

FMD=u. 


000167 

□ 01 73 

113* 

PUD. 


000170 

00171 

114* 

bO TO 920 


000171 

□ Cl 7 2 

115* 

201 IFI (LTI.EO. TIME) .AND. <ABS<TFI-T1» .LT..1) JGO TO 920 


01C173 

□ □174 

116* 

LTUTIME 


000211 

00174 

117* 

C 


010211 

□ 0174 

118* 

C CONVERT TA.TFO.TFI TO KELVIN 


000211 

00174 

119* 

C 


010211 

00175 

120* 

1 A-T A*273 


010213 

w 01 76 

121* 

TFO=TFO+273 


010216 

□ Cl 77 

122* 

TFI-TFI*273 


010221 

GC177 

122* 

C 


070221 

0C177 

12 4* 

C INITIAL TEMPERATURE AND flow RATE ESTIMATES 


010221 

00177 

125* 

c 


C10721 

0 2220 

126* 

TC=TA*CH/2u. 


0 r ‘lj22 4 

00201 

127* 

TI=ITC*TA)*.5 


000231 

00202 

128* 

TF~ITFI+TF0>*.5 


010233 


“n 

"U 


/ 

v. 
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00203 

129* 



TP=TI 

000236 

00204 

130* 



FMD=0. 

000237 

00205 

131* 



IFLU=0 

000240 

0C2 0b 

132* 



lFHAf<S(CM0-2.> .LTc.l) .AND.<TF0.6T.TFIH 1FLU = ! 

000241 

UOilO 

133* 



IFUFLU.NE.HGO 70 301 

000263 

CC212 

134* 



KO=NT*$PT*iTFO-Tm/(CW*CLJ 

000266 

0 0213 

135* 



FMU=MFM 

000277 

00214 

136* 



1F<P0.GT.0» )FHD=AMIN1IMFM, ,8*0H/R01 

0^0301 

0021b 

137* 


301 

CONTINUE 

000715 

00216 

138* 

C 



030315 

00216 

139* 

c 


ITERATE HEAT TRANSFER COEFFICIENT CALCULATION THREE TIMES 

000315 

OC216 

140* 

c 



OOt. 31 5 

00217 

141* 



LOOP=n 

000715 

00220 

142* 


40G 

CONTINUE 

C 0 u 7 1 6 

00220 

143* 

c 



000716 

00220 

144* 

c 


HTOP, HEAT TRANSFER COEFFICIENT and REA, REYNOLDS NUMBER 

000316 

00220 

145* 

c 



00G316 

00220 

146* 

c 



000316 

00221 

147* 



TSKY=.Q552*«TA**1.5! 

COO 31 6 

00222 

148* 



CALL CNVC JHCl,REA,TC,TA,WOfCL) 

CnC'32 3 

QC222 

149* 



IFJNG.GT.0.1G0 TO 401 

000333 

00225 

150* 



CALL RADCf HRl,TC,TSKr,EC f t.» 

000336 

00226 

151* 



hpi:hri*itc-tskyi/ctc-ta> 

000345 

0022 7 

152* 



HTGP=HC1*HR1 

000354 

002JC 

153* 



GO TO 402 

000356 

00231 

154* 


401 

hT0P=HTGLAS«NG,TA,TC,HC1,EC,EG,TLT» 

000360 

C-?2J2 

155* 


402 

CONTINUE 

000772 

0023 2 

156* 

c 



C0U372 

00232 

157* 

c 


HFIN HEAT TRANSFER COEFFICIENT 

000372 

00232 

158* 

c 



000372 

OC233 

159* 



CALL CNVC(HC2,REN1,TI,TA,U0,CL) 

000372 

QC23t 

lbO* 



CALL RADC(HR2,TI,TA,EP,1.1 

Q0U401 

0 0*35 

161* 



HFIN=1 ./<1./HI*1 »/tHC2*FAC*HR2J 1 

000410 

00235 

162* 

c 



000410 

00235 

lb 3* 

c 


HFLU, HEAT TRANSFER COEFFICIENT TO FLUID AND REF, REYNOLDS NUMBER 

00041C 

00235 

164* 

c 



000410 

00236 

165* ‘ 



HFLU-C • 

000423 

00277 

166* 



IFlIFLU.EOoTJGO TO 700 

000424 

CC241 

167* 



CALL FLUC(HFLU,REF,NT,DT,CW,COP,THP,FMD,DEN,TF,COCI 

000426 

0024 1 

168* 

c 



000426 

00241 

169* 

c 


EQUIVALENT BOTTOM TEMPERATURE TE90T AND HEAT TRANSFER COEFFICIENT MBOT 

000426 

00241 

170* 

c 



000426 

0 0242 

171* 


700 

CONTINUE 

Q 0 0 4 4 4 

3C*43 

172* 



HeOT = l . / 1 1,/HC + l ,./IHFIN4HFLUt 1 

OOC444 

00244 

173* 



TEbOT= IHFIn*TA*HFLU*TF >/» HFIN ♦HFLU) 

D00455 

00244 

174* 

c 



000455 

00244 

175* 

c 


UPDATE TEMPERATURE AND F LON RATE 

000455 

00244 

176* 

c 



000455 

3C245 

177* 



TC=<QH«HTOP*TA*HBOT*TEBOT )/ IHT0P*HB0T1 

000464 

00246 

178* 



TP=TC-HBOT*«TC-TEBOT »/HC 

O^OuVS 

00247 

179* 



TI=TP-HF1N*(TP-TA)/HI 

000502 

CG25U 

160* 



CFLU-HFLU*t TP-TF 1 

C7CJ507 

00251 

181* 



FMD = L . 

000513 

00252 

182* 



1FI0FLU.LE. 1. >G0 TO 800 

000514 

00254 

163+ 



IF<CFLU,GT. CMFM*fiO))GO TO 799 

000516 

00256 

164* 



FMQZCFLU/RO 

006523 

00257 

185* 



GO TO ecc 

QOC526 
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00263 

23b* 

799 FP*D=MF« 

900530 

0026 1 

187* 

R A-OFLU/HFM 

000531 

OC262 

IBS* 

TFsTFI*RA*CL*6W*.5/t$Pt*NT* 

000534 

00263 

139* 

800 CONTINUE 

000545 

00263 

190* 

C 

0Tj545 

00264 

191* 

L00P=L00P*1 

QOC545 

0026$ 

192* 

1FILOOP.LE.2IGO TO 400 

000547 

0-3265 

193* 

C 

000547 

00265 

194* 

C CHECK FOR EFFECTIVE FLUID COOLING 

0^0547 

00265 

ft 

195* 

C 

0G&54? 

Q0267 

196* 

1FI0FLU.GE.0.JG0 TO 900 

01355? 

Cl C2 7 1 

197* 

IFLU-0 

0PC555 

03272 

198* 

GO TO 400 

030556 

3 02 73 

199* 

900 CONTINUE 

OOCSfcC 1 

uC273 

2 jG* 

C 

0 0 C 5 6 0 

GG27J 

2« 1 * 

C OUTPUT CALCULATION 

ooqseo 

03273 

2u2* 

C 

D 00560 

00274 

203* 

TC=TC-272. 

J0&569 

0D275 

2D 4* 

TP=TP-273. 

036562 

U 02 7 1> 

205* 

Tl=TH-273. 

000565 

022 77 

206* 

T2=2.*TF-TFI-273. 

00u570 

30303 

207* 

7A=TA-272» 

030575 

JG3Hi 

2J8* 

TFlrTF 1-273. 

000630 

00332 

2 J9* 

TFC=TFO-273 0 

one 501 

03303 

210* 

Pl=CFLU*CL*C«/10C0. 

a oo <c-4 

30304 

211* 

REi-w. 

o^ceu 

0C3 05 

212* 

IF1 AfcS ICH0-J.1.LE..2 )RE1=. 0742*1 CCU*CL)**.2835>*t VD**.S67> 

003612 

30307 

213* 

IFJFPD.LE.O.IGO TO 939 

0*’06 37 

30311 

214* 

IF« 48 S 1 CmO-2. 1 .LE.. USE 1 = 7. 85E -11*1 FHD**2. 855 )* IDT** 1-4.7021 1 

003642 

03211 

215* 

1 *NT*CL 

900642 

0 0313 

216* 

909 CONTINUE 

C00667 

00314 

217* 

I FJ.IOP .NE.U10PD = RE*RE1 

QPD667 

30316 

218* 

920 IFCTIME.LT. THAXURETURN 

030674 

03320 

219* 

1FIImPL.LT. 21RETURN 

9* 1 C732 

00322 

220* 

CCAPrCCAP*CC*CL*CW 

ooo 7 u 

33323 

221* 

CM A -CM A ♦ C M 

030716 

03324 

222* 

CPOS=CPOS*CPO*OP 

0^0721 

00325 

223* 

RETURN 

0 r '0725 

00326 

224* 

END 

001267 


7.9 ONE DIMENSION TABLE LOOKUP 





T ab I es Description 

FTA Tabular values of function 


Input quantity 

ABS(AN) < 0.5 for equispaced interpolation 
(AN < 0 prevents extrapolation) 

Outputs 

Var i ab I e/Port 

F0 Output quantity 


I nputs 

Parameter/ Port 

FIN 

AN 


Calculation Sequence 

F0 = FTA(FlN) 


NOTE: A maximum of 18 points is allowed in the table. 


PRECEDING 


page blank not filmed 
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0 



0 





Revis ion Paaes 

l *ft l 


OS 










Sections 7.27 - 7.29 

{/ ^ s Replaces pages 243 to 266 of the original document. 

c . * . 



( 


PRECEDING PAGE BLANK NOT FILMED 


xx- 19 


zmW'T&r* 1 








09 








r> 








CO 









C0173 

79* 


4CS FORMAT! 1HQ»19H STATOR RESISTANCE ,F12.3,12H OR 0AI 

IPIN6 * 


000172 

o 

00173 

60* 


XF12.3,23H TOO HIGH FOR MOTOR 1 



000172 

00 

0017*1 

81* 


1FI1MPL.E0.2ITCNT?ICNT*1 



000172 

o 

00174 

82* 

£ 




000172 

r\s 

00174 

83* 

c 

EFFICIENCY AND MAXIMUM OUTPUT ROVER 



000172 


001741 

84* 

c 




00C172 


00176 

85* 


409 CONTINUE 



000201 


00177 

86* 


P2r0. 



000201 


IC200 

87* 


EF2=CFI 



00C201 


00201 

88* 


MP2TAMIN1 IMP! «RAP) 



000203 


00202 

89* 


60 TO 5CQ 



00021 1 

t* ' 

00203 

90* 


400 EF2=m*P2/Pl 



000213 

o 

002041 

91* 


HP2-AM1N 1 IH PI V R API *P2/P1 



000216 

E 

00205 

92* 


IF IRS .NE .0. )T0=P2*737.6/0ME6A 



000226 

C 

C020& 

93* 

c 




000226 

5 

QC207 

94* 


500 iniMPL.LE.llRETURN 



00023S 

o 

00207 

9 5* 

c 




CC0235 


00207 

96* 

c 

’STATISTICS 



000235 


PG207 

97* 

c 




C0023S 

Q 

00211 

98* 


MT-ANAxl (TOvMTI 



000243 

t=i 

00212 

99* 


MPNSAMAX1 IP2/RAP t HPNI 



060251 


CC213 

100* 


SP£SP«P2*7INC 



0CC260 

r 

00213 

101* 

c 




000260 

> 

00214 

102* 


iriTICE *L? •TMAX11RETURN 



000264 

z 

00216 

103* 


CCI-CCI*CC 



600273 

5* 

00217 

104* 


CHI9CN1«CM 



000276 


00217 

las* 

c 




000276 

Z 

00220 

106* 


RETURN 



000301 


00221 

107* 


END 



000433 


PA 

7.27 POWER ACCUMULATOR 


■►MAX. OUTPUT POWER (MPO) 

► POWER REQUESTS 

(RE1 ,RE2,RE3,RE4) 

► POWER OUTPUT (PO) 


This component sums power from four input ports and allocates power re- 
quests to each port's source of power generation. An input power request 
is allocated according to user-supplied weights within the ports of high- 
est priority. If an input power request (load) exceeds the maximum power 
that can be delivered by the ports of highest priority, then the remaining 
load is allocated to the next priority ports. (See 1.2.2 and 7c for 
further discussion.) 


POWER REQUEST (REO) — — ► 

INPUT POWER 
(PI ,P2,P3,P4) 

PRIORITY SEQUENCE 
(PS1 ,PS2,PS3,PS4) 


MAX. POWER 
(MP1 ,MP2,MP3,MP4) 


INPUT EFFICIENCY _ 
(EF1 ,EF2,EF3,EF4) 
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PA 

1 

Inputs 


Parameter/Port 

Description 

Units 

RE- 

0 

Load request , 

kw 

EF 

1,2, 3,4 

Input efficiency from port i 

- 

P 

1, 2,3,4 

Input power from port i (default =0.) 

// 

kw 

PS 

1,2,3, 4 

Priority sequence {default = 1*2, 3, 4) 

- 

F 

1, 2,3,4 

Allocation weight (for equal priorities) 

“ 

MP 

1,2, 3, 4 

Maximum power (default =0.) 

kw 

Outputs 



Variable/Port 



MP 

0 

Maximum deliverable power (2MP(i)) 

kw 

RE 

1,2, 3, 4 

Power request for port i 

kw 

P 

0 

Power output 

kw 

SP 


Supplemental power request to meet load 
(Power deficit = REq - 2 MP^-) 

kw 

Statistics 



SRE 


Sum of energy requested 

kwh 

PC 

1,2, 3,4 

Percent of cumulative load request 
delivered by port i 

% 

1 No 

capital costs 

assigned since this is an allocation component 

, not 


a physical device. 
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PA FAIR SHARE ALLOCATION 
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ro 

03 


SUBROUTINE PA 


ENTRY POINT 0C0514 


STORAGE USED S CODEtl) 000746; 


DATA (0) 000107; 


BLANK COMMON (2) 000000 


COMMON BLOCKS: 


0003 CIKPL 0C00Q1 

0C34 CSIMUL 000010 


Vi 


EXTERNAL REFERENCES (BLOCK, NAME) 
0005 NERR3S 


STORAGE ASSIGNMENT (fcLOCK, TYPE, RELATIVE LOCATION, NAME) 


0001 


000372 

100DL 

0001 


000375 

20001 

0001 


000231 

2346 

0001 


000240 

2416 

0001 


000260 

25*6 

CC01 


CGC3P2 

273G 

CDC1 


000321 

305G 

□ 001 


030357 

325G 

0001 


000103 

4GL 

0001 


000275 

402L 

OP 01 


□00413 

SOOL 

0001 


0002 75 

feOCL 

00U1 


000312 

703L 

C001 


00C155 

8 PL * 

0001 


000347 

BOOL 

QC01 


CLQ366 

900L 

GO 04 


corona 

OUM 

0000 

R 

000030 

FR 

GOOD 

R 

000342 

FRU 

0000 

I 

000037 

I 

0003 

I 

OCOLUC 

IMPL 

0000 


000054 

INJPS 

ooco 

I 

000041 

K 

CDOO 

R 

000334 

LL 

0000, 

R 

030035 

LOLD 

0000 

R 

000014 

HP 

0000 

R 

030304 

PR 

0003 

R 

030000 

R 

ODOO 

R 

D0D324 

SMP 

0000 

R 

030044 

SRI 

PC 03 

P 

000043 

SRO 

0000 

R 

•500020 

SW 

0004 

R 

CGC006 

TINC 

□ COO 

R 

030336 

TINC 1 

000 4 


030007 

TM AX 

00-0 

R 

000010 

U 

0000 

R 

0CCU40 

XI 









r- _ ' . 






00100 

l* 

CPA 



00101 

2* 


SUBROUTINE PA1MP0, 


CC101 

3* 


1 Rl, R2, R3 , R4, 


00101 

4 * 


2 PO,SP, 


30101 

5* 


3 SR,PC1,PC2,PC3,PC4, 

CD 

O 

OOAOl 

6* 


4 RO, 

30101 

7* 


4 EF1, EF2, EF3 , EF4 , , 

00 

CC101 

8* 


5 PI, P2, P3 , P4 , 

4* 

30101 

9* 


6 PR1, PR2, PR3 , PR4, 

o 

cam 

1C* 


7 Ul, W3, W 4 , 

00 

oaiai 

11* 


8 HP 1 , MP2, MP3, MP4) 

0 

1 

30101 

12* 

C 


ro 

00101 

13* 

C 

PURPOSE. MODEL POWER ACCUMULATOR 

30 

03101 

14* 

C 


a> 

30101 

15* 

C 

METHOD. PRIMARY RECUEST ALLOCATION RESULTING FROM PRIORITY 


30101 

16* 

c 

ASSIGNMENTS. SECONPARY RECUEST ALLOCATION RESULTING 


00101 

17* 

c 

F POM WEIGHT ASSIGNMENTS. 


02101 

18* 

c 

THAT IS, REOUESTS ARE ALLOCATED ACCORDING TO* 


3 0101 

19* 

c 

* PORT PRIORITY (HIGHEST PRIORITY = 11 


00101 

2 0* 

c 

* PORT WEIGHTS (IN CASE OF ECUAL PRIORTIES. ) 


30101 

21* 

c 



30101 

22* 

c 

FORMAL ARGUMENT DEFINITION. 


00101 

23* 

c 

Rl,..., R4 * POWER RECUESTS IN KW (OUTPUTS) 


.»■ A 
»• 


000003 
OOOPOO 
D30000 
000000 
OODOOP 
03CP00 
onopoo 
OPOPOP 
0 "10000 
OCuCCQ 
OPGPDO 
OPDOOO 
010030 
000003 
3P0C00 
D1JCC0 

or>ooao 

LP0030 

030000 

313030 

000030 

030300 

030300 
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0 


00101 

2 * 1 * 

c 

MPO * TOTAL MAXIMUM POWER 

1 OUTPUT) 

000000 

00101 

25* 

c 

SP * SURPLUS REOUEST 

(OUTPUT) 

OOCOOO 

00101 

26* 

c 

pn X TOTAL LOAD IN X H 

•OUTPUT ) 

OnOPQO 

00101 

27* 

c 

SR X SUM OF ENERGY REQUESTED, 

KWH IOUTPUT) 

O 0 O"OO 

0C101 

26* 

c 

PCI « . . . ,PC4 PERCENT OF CUM LOAD DELIVERED ICiUTPUTI 

0-70030 

00101 

29* 

c 

l»P X TOTAL POWER REQUESTED, KW 

IINPJJT) 

OSGOOD 

00101 

30* 

c 

PI,.,,, P4 X INPUT POWER IN AW 

•INPUTS) 

020000 

OCiOl 

31* 

c 

PR1 , . . . , PR4 X PORT PRIORITIES 

•INPUTS) 

o'lonoo 

001 01 x 

32* 

c 

hi,..,, W 4 X PORT WEIGHTS 

•INPUTS) 

o^oinn 

00101 

33* 

c 

MP1, mm MP4 X MAXIMUM POWERS 

1 INPUTS) 

QOUDOO 

00101 

34* 

c 

EF1, .... EF4 X EFFICIENCIES 

1 INPUTS) 

Q 10000 

00101 

35* 

c 

COMMON STORAGE 


002000 

0 C 101 

36* 


COMMON/ CIMPL / IMPL 


01 i)" 0 CI 

0010 * 

37* 


COMMON / CSIMUL / PUK 161, TINC, TM AX 


o^uioo 

00105 

36 * 


PEAL MP0,MP1,MP2,HP3,MP4 


DPC 00 O 

0G1C5 

39* 

c 



O1U0OO 

C 010 S 

40* 

c 

LOCAL VARIABLES 


C 10030 

UC1C3 

4 l* 

c 



CIGDQO 

00105 

42* 

c 



37O n 00 

C 0 1 0 5 

43* 

c 

RIK) IS THE POWER REOUEST AT PORT K 


000000 

00106 

44* 


REAL R(4) 

0 

01C03D 

00106 

45* 

c 



000000 

CC1C6 

46* 

c 

PRIM IS THE PRIORITY ASSIGNED TO PORT 

K 

OOulOO 

omr.7 

47* 


REAL PR 1 4 ) 


010100 

DC1C7 

46* 

c 


j 

01G1D3 

G 0 1 0 7 

49* 

c 

fa IK. ) IS THE WEIGHT ASSIGNED TO PORT K 


OOG n QO 

35110 

£ 0 * 


REAL W 1 4 ) 

• ' 

000000 

00110 

51* 

c 



000000 

00110 

52* 

c 

MPIM IS MAXIMUM POWER TO £iE ALLOCATED 

TO PORT K 

010000 

00111 

53* 


REAL MPIM 


OOC 0 OO 

00111 

54* 

c 



000000 

com 

55* 

c 

SWIIl IS THE SUM OF THE WEIGHTS ASSIGNED TO PRIORITV-I PORTS 

00C100 

com 

56* 


REAL Sfa 1 4 ) 


C 0 U 0 OO 

00112 

57* 

c 



OlcOCO 

001 12 

58* 

c 

SM»III IS THE SUM OF THE MAXIMUM POWER 

AT PRIORITV-I PORTS 

30 COOP 

iiClU 

59* 


REAL SMPI4) 


000300 

00113 

60* • 

c 



Q0GOCO 

00113 

61* 

r 

w 

FPU IS FAIR SHARE UNIT FOR PRIORITY-I 

PORTS 

O 0 G 0 OO 

00113 

62* 

c 



000000 

00113 

63* 

c 

FRIM IS THE COMPUTED FAIR SHAPE REQUEST FOP PORT X. 

000-100 

0C1 14 

64* 


REAL FRIM 


oocroo 

00114 

65* 

c 



ODU'000 

00114 

66 * 

c 

LL IS THE LOAD LEFT AT EACH POINT IN THE ITERATION 

000*30 

00115 

67* 


REAL ll.lolo 


OILIUQ 

001 15 

68 * 

c 



010000 

00115 

69* 

c 

IF IMPL IS 2ER0, THEN ASSIGN DEFAULT VALUES 

030000 

00116 

70* 


IF 1 IMPL .GT. 01 GO TO 40 


O 1 OP 00 

0 C 12 U 

71* 


RC - 0,0 


0 00002 

UC121 

72* 


IF 1 PR 1 .EC. 0.99999) PR1 = l.C 


000007 

00123 

73* 


IF IPR2 .EQ. 0.99999) PR2 = 2,0 


000010 

0 012 5 

74* 


IF IPR 1 .EO. C. 95999) PR3 = 3.U 


010015 

00127 

75* 


IF IPR4 .EO. 0.99999) PR4 = 4.U 


01 C *22 

UC131 

76* 


IF IMP1 .EQ. C. 99999) HP1 = 0 


O 0 CD 2 ? 

00133 

77* 


IF IMP 2 .EQ. 0.99999 ) MP2 = 0 


000133 

C0135 

76* 


IF IMP 3 .EC. 0,9 9999) MP3 = C 


010137 

00137 

7 5* 


IF IMP4 .EQ. C. 99959) MP4 = 0 


0*0043 

UC141 

60 * 


1FIP1 .EQ. .99999) P1=C.0 


Q0CC47 
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00143 

61* 



IRS 8 * ,E0 • ,999991 P2=0.3 



000053 

00145 

62* 



IF*?.? .EC. ,999991 P3s U«0 



000057 

00147 

83* 



1F4P4 «£Q o ,99999) P4=0.P 



000063 

00151 

64* 



SP = 0. 



000067 

00152 

85* 



PCi=0. 



000070 

GC153 

66* 



PC2 = £i. 



003071 

G01 54 

87* 



PC3-0, 



000072 

00155 

88* 



PC4=0e 



OTOC73 

00156 

89* 



TINC1- 0.5*TINC 



0'»CC74 

00157 

90* 


40 

CONTINUE 



C3G100 

00157 

91* 

C 





GOG TO 3 

GC157 

92* 

C 


IF THE TOTAL MAXIMUM POWER IS ,LE . TOTAL POWER 


000103 

00157 

93* 

c 


REQUESTED • THEN SUBMIT REQUESTS AT MAX-POWER, SET REQUEST 


me too 

0C157 

94* 

c 


SURPLUS EQUAL TO THE DIFFERENCE, AND 

RETURN 


030100 

30160 

95* 



PC = PI ♦ P2 ♦ P3 ♦ P4 



010100 

00161 

96* 



1F4PR1 ,LE .0.01 MPI=D. 



0'*3104 

0C163 

97* 



lFiPR2.LE.0.C> HP2=0. 



070110 

00165 

98* 



1 F ( Pk 3 .LL *0 .0 1 MP3=C. 



070114 

00167 

99* 



1F4PR4 ,LE ,0,0 1 MP4 - 0 » 



010123 

00171 

ICC* 



MPG = HP 1 ♦ MP2 ♦ MP3 ♦ MP4 



070124 

GC172 

101* 



IF ( h P L. .GT. RO) GO TO 80 



070131 

0017<t 

102* 



R1 = MPl/EFl 



070134 

00175 

103* 



R2 = HPJ/EF2 



030137 

00176 

1 04* 



fi 3 = HP3/EF3 



0 n 01 42 

GGl 77 

10 5* 


*• 

R4 r MP4/EF4 



010145 

00200 

106* 



SP = RC - MPC 



033153 

00231 

107* 



GO TO SCO 



030153 

0 0232 

1J8* 


80 

CONTINUE 



D7G1S5 

0C2 32 

139* 

c 





070155 

00232 

113* 

c 


PROCEEC WITH ALLOCATION ALGORITHM SINCE THE SUM OF 


000155 

0C232 

111* 

c 


ALL MAXIMUM POWER INPUTS EXCEEDS THE 

TOTAL REQUEST RD 


0 D 1 5 5 

30202 

112* 

c 





C70155 

002,32 

113* 

c 


INITIALIZATION 



030155 

00233 

114* 



LL = RO 



03C155 

30204 

115* 



fil = 0.0 



030156 

00235 

116* 



R 2 = 0.0 



07015? 

00236 

117* ‘ 



R3 r 0.0 



0 3 0 T 6 0 

•30237 

118* 



R4 - 0.0 



0 n u 1 5 l 

Q C2 1 3 

119* 



SP = 0.0 



010162 

00210 

120* 

c 





00U1S2 

00210 

121* 

c 


IF THE TOTAL REQUEST 10R LOAO) IS ZERO, THEN RETURN 


000162 

00211 

122* 



IF IRO .LE, 0.0> GO TO 500 



010163 

30213 

123* 



R • 1 1 -R 1 



C 7 0 1 6 5 

002 1 4 

124* 



R «2)=R2 



030167 

00215 

125* 



R(3)=R3 



010171 

uG2 l 6 

12 6*. 



R«4l=R4 



010173 

30217 

127* 



PR t 1 > : PRI 



070175 

00220 

128* 



PR t 2 J = PR2 



07C177 

30221 

129* 



P R ( 3 ) = PR3 


v 

Q10201 

GG222 

13 0* 



PR 4 4 ) = PR4 



010203 

00223 

131* 



W (1) = W1 



070205 

30224 

132* 



W<2) = W2 



OOC237 

00225 

133* 



W ( 3 1 = W3 



030211 

Li 02 26 

134* 



W (4 ) = U4 



07C213 

00227 

135* 



HP { 1 } = M P 1 



07C715 

00231 

136* 



HP ( 2 1 : HP 2 



070217 

J 02 3 1 

137* 



HP(3J - HP 3 



07C221 


/■' 

* \ .. .. V - 
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00232 

138* 



HPC41 = MP4 

000223 

00232 

139* 

c 



Q7U223 

00232 

140* 

c 



000723 

00232 

141* 

c 


ITERATE ON PRIORITY I FOR I = 1. 2, 3. 4 

000223 

00232 

142* 

c 



000223 

00233 

143* 



00 10C0 I = 1, 4 

000231 

00233 

144* 

c 



07C231 

JC236 

145* 



X I - I . 

0"0?31 

0023a 

146* 

C 


OBTAIN SUM OF MAXIMUM POWER FOR PORTS WITH PRIORITY I 

Q30231 

0 02 3 7 

147* 



SMPII) = 0.0 

070234 

00&4G 

me* 



00 ICO K = 1. 4 

000240 

00243 

149* 



IF 1 PR IK ) ,E0. XI) SMPII) = SMPII) ♦ MPIK) 

000240 

0C245 

150* 


100 

CONTINUE 

000247 

00245 

151* 

c 



G7D247 

00245 

152* 

c 


IF NO PRIORITY-I MAXIMUM POWER EXISTS, THEN PROCEED WITH 

Q 7 0 2 4 7 

00245 

153* 

c 


THE NEXT HIGHER PRIORITY 

005247 

00247 

154* 



IF 1SMPII) .EG « 0.0) GO TO 1C00 

Q3Q747 

00247 

155* 

c 

- * 


000247 

0024 7 

156* 

c 


IF THE SUM OF ALL PRIORITY-I MAXIMUM POWER .6T. LOAD 

Q0G247 

0024 7 

16 7* 

c 


LEFT, THEN GO AROUND 

020247 

00251 

158* 



IF ISMPC1) »G T • l.L) bO TO 4 00 

’ 0707E1 

00251 

159* 

c 



37D251 

00251 

160* 

c 


THE SUM OF A15T PRIORI! Y-I MAXIMUM POWER «LE. LOAD 

07G251 

OC251 

161* 

c 


LEFT, SO SUBMIT EACH PRIORITY-I REGUEST 

07U251 

□ C253 

162* 



CO 200 K — 1 , 4 

07C260 

00256 

163* 



IF IPR IK ) .EG. XI) RIK) = MPIK) 

000260 

00260 

164* 


200 

CONTINUE 

00C266 

0C26 j 

165* 

c 



070266 

00260 

16 6* 

c 


UPDATE LOAO LEFT 

000266 

00262 

167* 



LL = LL - SMPII) 

D7D266 

00262 

168* 

c 



□70266 

JC262 

169* 

c 


IF THE REMAINING LOAD IS ZERO, THEN EXIT THE ITERATION 

070266 

C0263 

'-^7G* 



IF ILL .LE. O.U) GO TO 200Q 

0 ^ 0 2 7 1 

00263 


c 



07C271 

□ 026 3 

172* 

e 


OTHERWISE, PROCEED WITH NEXT HIGHER PRIORITY 

3 n C 2 7 1 

5026 5 

173* 



GO TO 1000 

070273 

□ 0265 

174* - 

c 



000273 

00266 

175* 


4QC 

CONTINUE 

Q30275 

00266 

176* 

c 



070275 

0C26b 

177* 

c 


THE SUM OF THE PRIOPITY-I MAXIMUM POWER EXCEEDS THE 

07G275 

0C266 

178* 

c 


LOAD LEFT, SO COMPUTE AND SUBMIT FAIR SHARE REQUESTS 

070275 

u 0266 

179* 

c 


TO EACH PRIORITY-I PORT 

000275 

00266 

16 0* 

c 



OOG275 

0 0267 

161* 


600 

CONTINUE 

070275 

0 0267 

182* 

c 



070275 

00267 

183* 

c 


SAVE LL FOR LATER REFERENCE 

0702 75 

□ 02 70 

18 4* 



LOLD = LL 

00C275 

002 70 

185* 

c 



C7C275 

00270 

186* 

c 


DETERMINE FAIR SHARE UNITS FOR ALL PRIORZTY-I 

G' 1 Q?75 

□ 02 7 j 

167* 

c 


PORTS TO WHICH NO REQUEST HAS BEEN SUBMITTED 

070275 

00271 

188* 



sum = 0.0 

QDC276 

00272 

189* 



DO 70C- K = 1, 4 

37C302 

00275 

19C* 



IF IRIK) .NE. C.UI GO TO 7C0 

070302 

00277 

191* 



IF IPR (< ) .EQ« XI) SW(I) = SWI1) ♦ WIK) 

070303 

00331 

19 2* 


70C 

CONTINUE 

070315 

00333 

193* 



FRU = i.G / SWII) 

o^om 

00303 

194* 

c 



D70313 



00333 

195* 

C 

FIRST, SUBMIT FAIR SHARE REQUESTS 

TO 

PORTS FOR 

WHICH THE 

coons 


00333 

196* 

C 

FAIR SHARE REQUEST EXCEE0S THE MAXIMUM POWER, 

CONSIDER ONLY 

070312 

ro 

0C333 

197* 

c 

PORTS TO WHICH NO REQUEST HAS BEEN 

SUBMITTED 



030313 

cn 

ro 

CC33H 

196* 


DO 6C0 X : 1, 4 






0'»L , 21 


003 3 7 

199* 


IF (RIM .hF. C.J1 GO 

TO PkO 





0^0 321 


00311 

2CjL* 


IF (PR 1 K 1 .'/L. Xll GO 

TO B jfl 





00U322 


0031 1 

2u l* 

c 







07C322 


0031 1 

20 2* 

c 

compute fair share 






070*22 


003 13 

2u3* 


FR(X) = (WIKI • FRUl 

* LL 





030725 


00313 

2 c 4* 

c 







CC0I25 


0C31 3 

205* 

c 

IF FAIR SHAPE EXCEEDS 

MAXIMUM POWER, 

THEN SUBMIT 

REQUEST 

0 7032 5 


0C314 

206* 


IF (FRIX) c GE »■ “Pm 1 

R ( X ) r APIA) 





030331 


0C3l<* 

207* 

c 

- AND REDUCE LOAD 

LEFT TALLY 





030331 


0031b 

2 ee* 


IF (FR(X) .GE. “?(Xl» 

LL = LL - HP ( K 1 




03D327 


303 20 

239* 

8 00 

CONTINUE 






D 30 750 


oci2e 

210* 

c 







030753 


uC3 2-j 

211* 

c 

IF LL .Tit. LOLO, THEh 

LL WAS REDUCED 

DURING THE 


030750 


0032 0 

212* 

c 

PPOCcSSING IN THE DO 

SEu LOOP ABOVE 

. THIS CHANGES 

070350 


u332o 

2 1 3 * 

c 

THE FAIR SHARE COMPUTATION . IT IS 

THEREFORE 



00c350 


00320 

214* 

c 

NECESSARY TO GO BACK. 

THROUGH THE 

DO 

830 LOOP 

IN 


O^OTtO 


00320 

215* 

c 

ORDER TO RECONSIDER ANY PORT WHICH 

MAY NOW 



070353 


00320 

216* 

c 

SAT ISFY THE REQUIREMENT THAT F ft ( X ) 

.GE. hP ( K ) • 

ONLY 

0C-C350 


003 2 0 

217* 

c 

PRIORITY-! PORTS TO WHICH NO REQUEST 

HAS BEEN 



ono3so 


0D37C 

216* 

c 

MADE ARE ELIGIBLE FOR 

RECONSICEPA T1 ON 




GC3353 


03322 

219* 


IF (LL .LT. LOLD) GO 

TO 





070350 


0C322 

22 0* 

c 







030350 


OC322 

221* 

c 

FINALLY, SUBMIT REQUESTS TO THOSE 

PORTS FOR WHICH 

THE FAIR SHARE 

030*50 


dC 322 

2i 2* 

c 

.lt. than their maximum po»er. consider only 



07C350 


00322 

223* 

c 

PRIOftITV-I POPTS to WHICH NO REQUEST 

HAS BEEN 

SUBMITTED 

30-C3ED 


dG324 

224* 


DO 9l7 R = ), 4 






000357 


00327 

2c5* 


IF (R(M .NE. E* ,0) GO 

TO 9^0 





030757 


0C231 

226* 


IF (PR (X) .NE. X I) GO 

TO 900 





030360 


o 03 3 J 

227* 


ft (x 1 : F ft ( x, ) 






03u*fc3 


aC334 

226* 

9C0 

CONTINUE 






030367 


OC33o 

2^' 9* 


LL=C.n 






Q7C367 


0C337 

230* 


DO TO 2000 






030370 


0 03 37 

231* 

c 







072373 


0 03S 0 

23 2* 

1000 

CONTINUE 






C7C375 


00343 

23 2* 

c 







03077s 

ro 

QC'34 2 

234* 

2000 

CONTINUE 






D 70 775 

o 

co 

c03 4 2 

23 5* 

c 







030375 


u C 3 4 2 

236* 

c 

FINALLY, ASSIGN OUTPUTS TO NON-SUBSCRIPTED FORMAL 

PARAMETERS. 

030375 

o 

30342 

237* 

c 

ALSO, MODIFY ALL REQUESTS ACCORDING TO THE INPUT 

EFFICIENCIES 

070775 


0034 3 

233* 


ft 1 = Rill / EF1 






030 7-7 5 

CD 

o 

C034 <4 

239* 


S 2 = R ( 2 1 / EF2 






C-30377 

f 

&C3 4 b 

24 C* 


ft 3 - R (3) / EF3 






070.40 2 

ro 

2234b 

241* 


ft 4 = R ( 4 ) / EF4 






0-75405 


00347 

24 2* 


SP r LL 






03-04 10 

fD 

< 

0235 0 

243* 

500 

IF(lMPL.LE.l) RETURN 






030413 

* 

00352 

24 4* 


SRO= SR 






DCC421 


OC353 

245* 


SP=Sfi* Pl»TINC1 






070423 


00354 

24 6* 


IFISR.LE.C. ) RETURN 






Q7D427 


3035b 

24 7* 


SRQ -SRC/S R 






030435 


003 5 7 

246* 


SPI= TIi:Cl«lDC-./SP 






0 n 0440 


033 bU 

249* 


PCi= PC1*SP0 ♦ P1*SPI 






0 7044 4 


0726 1 

250* 


PC2- P C 2 * S R 0 * P 2 * S P 1 






C 7 0 4 5 0 


C3362 

251* 


PC3" P C 7* SRO ♦ P3*SRI 






000456 


■. ■ A 

i; f 


A 

v~ 
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OC363 252* PC«: PC«*SRO ♦ P«*SRI 

0036% 253* RTTURN 

QC36S 251* £ND 


00C*6* 

OOQH72 

0307815 


PO 

cn 

to 


7.28 POWER DIVIDER 


PD 



This component allocates power to four ports plus surplus based on priority 
and allocation weights for equal priority ports. Each port is assigned a 
priority sequence from 1 to 4, and a weighting F-j >0, i=l,2,3,4 for propor- 
tional allocation among equal priority ports. If power available exceeds 
the power requested for the ports of highest priority, then the remaining 
power is allocated to ports having the next highest priority. If power 
available is less than the power requested for ports of equal priority, 
then power is allocated among them in proportion to their respective allo- 
cation weights. 

The total power request is the sum of the port requests divided by input 
efficiency. The maximum power outputs MP1,...MP4 are necessary for direct 
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PD 

connections to a power accumulator PA. These variables may be used as 
maximum power inputs to other components, although such connections are 


not required. 

(See 1.2.2 and 7c for further discussion.) 


Inputs* 



Parameter/Port 

Description 

Units 

P 0 

Input power 

kw 

RE 1,2, 3, 4 

Power request of output ports 

kw 

PS 1,2, 3, 4 

Priority sequence (default - 1,2, 3, 4) 

kw 

F 1, 2,3,4 

Allocation weight (for equal priorities) 

- 

MP 

Maximum input power (default = P0) 

kw 

EF 

Input efficiency 

- 

Outputs 



Variable/Port 



P 1,2, 3, 4 

Output power for port i 

kw 

RE 0 

Output power request 

kw 

MP 1,2, 3, 4 

Output maximum power based on MP 

kw 

Statistics 



SP 

Surplus power 

kw 


* No capital costs assigned since this is an allocation component, 
not a physical device. 
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CALCULATION LOGIC 


• WT(k) = 1 k ■ 1,2, 3, 4 

• IF PS(k) = 0, RE ( k ) =0 k = 1,2, 3, 4 

• REO = ERE(k) 


REO£ PO ? 


•INITIALIZE ALLOCATION LOGIC 

P ( k ) - 0 k » 1 ,2,3,4 
1=1, PL = PO, MPA = MP 


r * SUM REQUESTS SR AND WEIGHTS 

WT FOR ALL PORTS WITH PRIORITY I 
MP(k) = MPA*F(k)/WT 
'MPA = MAX(MPA-SR,0.) 


1 = 1+1 


REO-PO OR 


SR < PL ? 


YES , 


• GO TO FAIR SHARE 
ALLOCATION (3) 
(NEXT PAGE) 


FOR EACH PRIORITY I PORT, SET 
P(k) “ RE(k) 


UPDATE POWER AVAILABLE PL 


GO TO 3 
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PD FAIR SHARE ALLOCATION 


FOLD = PL 


• DETERMINE THE SUM SW OF WEIGHTS 

F ( K) ASSIGNED TO PRIORITY - I PORTS 
NOT PREVIOUSLY ALLOCATED 


• DETERMINE FAIR SHARE FR(K) FOR PRIORITY - I, 

NON-ALLOCATED PORTS. 

FR(K) = (F(K)/SW)* PL 

• IF FAIR SHARE FR(K) EXCEEDS REQUEST RE ( K) , 

THEN ALLOCATE AMOUNT REQUESTED AND REDUCE PL 
BY AMOUNT ALLOCATED. 


♦ ALLOCATE FAIR SHARE OF REMAINING 
POWER TO PRIORITY I PORTS NOT 
PREVIOUSLY ALLOCATED 


SP = PL IF IMPL>1 


REO = MIN (REO ,MP)/EF 
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ro 

17 ! 

co 


SUBROUTINE PD 


ENTRY POINT 000*53 


STORAGE USED: CODEtlt 000631; DATA (0) 000105; BLANK COMMON ( 2 ) GOOOOD 
COMMON BLOCKS: 
uOU3 CIMPL 00Q0C1 


EXTERNAL REFERENCES 1BL0CK, NAME) 
LiOCA NERR3S 


STORAGE ASSIGNMENT I BLOCK , TYPE, RELATIVE LOCATION. NAME) 


UOO 1 


C 4)04 02 

IC3DL 

000! 


000166 

2065 

0001 


uOO 1 


QCC333 

275G 

0001 


000370 

315G 

0001 


L-001 


DL03C7 

6C0L 

00G1 


0QC324 

7 3 CL 

0001 


0000 

R 

CuOJ3U 

FR 

on jc 

R 

000C43 

FRU 

CP'30 

i 

gOOO 

I 

00004 1 

K 

coco 

R 

POCO0C 

P 

C'CiiO 

R 

CjCQO 

R 

OOCOIO 

PR 

0000 

R 

000004 

R 

coca 

R 

COCO 

R 

000040 

WT 

coco 

R 

000037 

XI 




000176 

214G 

0001 


000274 

246G 

0001 


000314 

2636 

000033 

4CL 

0001 


000307 

40OL 

0001 


000105 

60L 

000077 

SOL 

00Q1 


000361 

80 CL 

0001 


000377 

9Q0L 

000036 

I 

0003 

I 

DDCOOO 

IMPL 

0000 


0C0052 

INJPS 

PGD034 

PL 

0000 

R 

000335 

PMA 

000 0 

R 

000042 

POLO 

C00024 

SR 

CODO 

R 

000020 

SW 

0300 

R 

Oi;O0l4 

w 



00100 

1* 

CPD 


O0COQD 


00101 

2* 


SUBROUTINE PD« 

O0ODOO 


00101 

3* 


1 PI, P2, P3» P4, 

Q0C0OO 


0C101 

4* 


2 RO, 

O0CO3O 


00101 

5* - 


3 SP,PM1,PH2,PM3,PM4* 

000000 


00101 

6# 


4 PO, 

OOC0OO 


3G1D1 

7* 


5 R 1 , R 2, R3 , R4 , 

onccoo 

CD 

OGiai 

8* 


6 PR 1 , PR2, PR 3 , PR4, 

000000 

n 

c/> 

00101 

9* 


7 Wl, 42, U3, W4 , PM,EF) 

aooooo 

fv 

OC101 

1C* 

C 


onoooo 

o 

00101 

11* 

C 

PURPOSE. MODEL POWER DIVIDER 

DO'Jf'OO 

00 

00101 

12* 

c 


000000 

o 

00101 

13* 

c 

METHOD. PRIMARY FLOW ALLOGATION RESULTING FROM PRIORITY 

000.000 

1 

f\3 

QC101 

14* 

c 

ASSIGNMENTS. SECONDARY FLOW ALLOCATION RESULTING 

O0CCUQ 


caioi 

15* 

c 

FROM WEIGHT ASSIGNMENTS. 

300000 

70 

CD 

G01G1 

16* 

c 

THAT IS, TOTAL AVAILABLE POWER IS ALLOCATED 

O0O"0O 

< 

0 3101 

17* 

c 

ACCORDING TOi 

003000 

9 

OC101 

18* 

c 

* PORT REQUESTS 

O0G0D3 


00101 

19* 

c 

* PORT PRIORITY HIGHEST PRIORITY = 1» 

C 00.000 


0 3101 

20* 

c 

* PORT WEIGHTS 1 1 n case OF EQUAL PRIORTIES) 

oo^nao 


001.01 

21* 

c 


000003 


cnioi 

22* 

c 

ALLOCATION SCHEME. 

oraroo 


acioi 

23* 

c 

IS SUM OF ALL REQUESTS .LT. POWER AVAILABLE P 0\ 

000 "00 


00131 

24* 

c 

YES • 

000030 


uoiai 

25* 

c 

FULFILL EACH REQUEST 

0-0 0030 


& h 

\ Y 



P » 

•,'trr 





« 


_,.«****, 

a 



00101 

26* 

C 

UPDATE POWER AVAILABLE 


00101 

27* 

c 

EXIT 

00 

U0101 

28* 

c 

NO. 

o 

C/5 

G0101 

29* 

c 

IS SUN OF ALL PRIORITY-1 REQUESTS .LT. POX 


00101 

30* 

c 

YES. 

o 

00101 

31* 

c 

FULFILL EACH PRIORITY-1 REQUEST 

co 

0C131 

32* 

c 

UPDATE POWER AVAILABLE {TO PL) 

o 

OOlOi 

33* 

c 

60 ON TO PRIORITY-2 REQUESTS 

f\J 

saini 

34* 

c 

NO. 


0010 1 

35* 

c 

ALLOCATE FAIR SHARE TO EACH PRIORITY-1 PORT 

70 

a> 

00101 

36* 

c 

EXIT. 

< 

00101 

37* 

c 

IS SUN. OF ALL PRIORITY-2 REOUESTS .LT. PL 


00101 

36* 

c 



U01CI 

39* 

c 

AND SO GN AND SO FORTH 


00121 

40* 

c 



U0101 

41* 

c 

FORMAL ARGUMENT DEFINITION. 


00101 

42* 

c 

PI e • * . • P 4 X POWER ALLOCATIONS IN KW (OUTPUTS) 


GC101 

43* 

c 

RC X TOTAL POWER REQUESTED (OUTPUT) 


0C1CI 

44* 

c 

SP X SURPLUS POWER (OUTPUT) 


iiClOl 

45* 

c 

PMI f * . ■ .PM4X PORT MAXIMUM OUTPUT POWER IN KW (OUTPUT) 


G0101 

46* 

c 

PO X TOTAL POWER INPUT IN KW (INPUT) 


00101 

47* 

c 

PM X MAXIMUM INPUT POWER IN KW (INPUT) 


acioi 

46 * 

c 

EF X INPUT EFFICIENCY (INPUT) 


00101 

49* 

c 

R1....I R4 X PORT REQUESTS IN KW (INPUTS) 


C0101 

5C* 

c 

PR1 * * » • » PR4 X PORT PRIORITIES (INPUTS) 


00101 

51* 

c 

hl««... W4 X PORT WEIGHTS (INPUTS) 


00101 

52* 

c 



00101 

53* 

c 

common STORAGE 


00103 

54* 


COMMON/ CIMPL / IMPL 


00103 

55* 

c 



00103 

56* 

c 

LOCAL VARIABLES 


UG103 

57* 

c 



00103 

56* 

c 

P (K ) IS THE POWER ALLOCATED TO PORT K 


00134 

59* 


REAL P ( 4 } 


00104 

bQ* 

c 



0Q1C4 

61* 

c 

R(K) IS THE POWER REQUEST AT PORT K 


00105 

62* • 


REAL R(4) 


00105 

63* 

c 



001P5 

64* 

c 

PR4K) IS THE PRIORITY ASSIGNEO TO PORT K 


00106 

65* 


REAL PR ( 4 ) 


001 Ob 

66* 

c 



OGlOb 

d7* 

c 

k(H) IS THE WEIGHT ASSIGNED TO PORT K 


00107 

66* 


REAL W ( 4 ) 


00107 

69* 

c 



001C7 

70* 

c 

SWiI) IS THE SUM OF THE WEIGHTS ASSIGNED TO PRIORITY-I PORTS 


UOIIG 

71* 


REAL SW ( 4 ) 


00110 

72* 

c 


ro 

03110 

73* 

c 

SR < I ) IS THE SUM OF THE REQUESTS AT PRIORITY-I PORTS 

VO 

30111 

74* 


REAL SR ( 4 ) 


00111 

75* 

c 



00111 

76* 

c 

FRU IS FAIR SHARE UNIT FOR PRIORITY-I PORTS 


ilCUl 

77* 

c 



00111 

78* 

c 

FR4K 1 IS THE COMPUTEO FAIR SHARE ALLOCATION TO PORT K 


QC112 

79* 


REAL FRU) 


00112 

8 0* 

c 



00112 

61* 

c 

PL IS THE POWER LEFT AT EACH POINT IN THE ITERATION 


05113 

82* 


REAL PL 



ootinoo 

ooonoo 

oocooo 

oodpoo 

opcnoa 

CPOPPO 

onoroo 

ononoo 

OPuPOO 

oncnoo 

OPDPOD 

000300 

C0CCC0 

OPCiOPO 

opopjo 

oogpoo 

DTL'PGO 

OOLPOO 

OOU0OO 

OPDOOO 

ooocon 

ccoao 

DOGOOO 

ooccoo 

QCGPOO 

000300 

300000 

OOCCOO 

OOG'OQO 

0*10030 

Q0C0SO 

OPG^GP 

030000 

O0C n CO 

000020 

opohoo 

O0DPCO 

OPUOSO 

D0OPCC 

300030 

opcpna 

OPO COO 
COGPOO 
000030 
O0OPOO 
000000 
opgpoo 
OP CHOP 
CPOPOO 
DPL-nop 
000002 
cianGp 
OPGPOO 

ononoo 

oppnon 

oponoo 

13POPOO 
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00113 

83* 

e 





00113 

84* 

c 


IF IMP). IS ZERO, THEN 

ASSIGN DEFAULT 

VALUES 

00114 

85* 



IF «1HPL .ST, 01 SO TO 

40 


00116 

66* 



R1 = 0.0 



00117 

87* 



R2 s 0.0 



D012C 

68* 



R3 = 0.0 



00121 

89* 



R4 : D.P 



30122 

90* 



IF IPR1 .EQ. 0.99999* 

PR1 = 1«0 


3C124 

91* 



IF <PR2 .EQ • 0.99999) 

PRZ = 2.0 


0012b 

92* 



IF «PR3 .EC. 0.999991 

PR3 = 3o0 


0Cl3u 

93* 



IF IPR4 .FO. 0.99999) 

PR4 = 4.0 


00133 

94* 

c 





DC132 

95* 


40 

CONTINUE 



00132 

96* 

c 





05132 

97* 

c 


IF THE TOTAL POWER REQUESTED IS .LE. 

TOTAL POWER 

00132 

98* 

c 


INPUT, THEN SATISFY REQUESTS, SET POWER SURPLUS 

00132 

99* 

c 


EQUAL TO THE DIFFERENCE, 


05133 

let* 



1FCPR1.LE.0.0I Rl = Cf.U 



00135 

1GI* 



1FIPR2.LE.U.0) R2=0.U 



QQ137 

102* 



IFIPR3.LE.C.0) R3=0.0 



00141 

103* 



1FIPR4.LE.0.0) R4-0.C 



00143 

104* 



RC = R1 ♦ R2 R3 * R4 


00144 

lu5* 



IF ) RC o6T. POI CO TO 

80 


00146 

It 6* 



PI = R 1 



0014? 

It 7* 



P2 = R2 


. 

00150 

IDS* 



P 3 = R 3 



00151 

109* 



P 4 = P 4 



3C152 . 

110* 



PL = PO -R0 



00153 

111* 



GO TO 60 



00154 

112* 


80 

CONTINUE 



00154 

113* 

c 





0 Cl 54 

114* 

c 


PROCEED WITH ALLOCATION ALGORITHM SINCE THE SUM 0 

00154 

115* 

c 


ALL REQUESTS EXCEEDS THE TOTAL AVAILABLE POWER PO 

00154 

116* 

c 





0515 4 

117* 

c 


INITIALIZATION 



00155 

118 * 



PL = Pt 



UG156 

119* ‘ 



pi = c.n 



00157 

12 C* 



P2 = D.O 



0C16O 

121* 



P 3 S 0.0 



QClbl 

00161 

122* 

123* 

c 


P4 = 0.0 


.98 

00162 

124* 


60 

PMA = PM 


•n a 

00163 

125* 



IF) PH «EQ • 0.99999) 

PM A = PO 


00165 

U6* 



Pll) = PI 


O > 

05166 

127* 



P ( 2 ) = P2 


30 r- 

00167 

128* 



Pt3) = P3 


O TS 

c > 

00170 

129* 



P ( 4 ) = P4 


00171 

13C* 



R < 1 ) = R1 


> 0 

00172 

131* 



R 12 t = R2 


n rn 

00173 

132* 



fi(3) = R3 


3 771 

00174 

133* 



RC4) - R4 


•< (A 

uC175 

134* 



PR C 1 ) = PR1 



J 3 1 7 o 

135* 



PR l 2 ) = PR2 



00177 

136* 



ppm = PR3 



00200 

137* 



PR 1 4 ) = PR4 



00231 

138* 



fa 1.1 ) -S Wi 



0 C2 02 

139* 



W (2 ) = U 2 
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00203 

140* 


w<31 r W3 

C001S6 

0020<t 

141* 


414) = W4 

000163 

□ 0234 

142* 

C 


030160 

00234 

143* 

C 


O^O 160 

QC204 

144* 

C 

ITERATE ON PRIORITY I FOR I : 1* 2. 3. 4 

033160 

□ 0234 

145* 

c 


330160 

00235 

146* 


oo lona 1=1,4 

030166 

□ 0235 

147* 

c 


C^C 166 

03210 

143* 


XI = I 

030166 

00210 

149* 

c 

OBTAIN SUM OF REQUESTS FROM PORTS WITH PRIORITY I 

£nC166 

oohi 

150* 


SR C 1 1 = G.O 

030171 

0G212 

151* 


WT=C.Q 

31C172 

0C2 13 

152* 


DO luO K = 1, 4 

Q70176 

0321b 

153* 


IF IPPIK) *EQ • XII SR)I> = SR ( I ) ♦ RIKI 

0 30 176 

0022u 

154* 


IFIPRCh) «EQ. XII WT= NT* ylK) 

030233 

00222 

155* 


100 CONTINUE 

03C213 

GC222 

156* 

c 


03u?l 3 

03224 

157* 


I F (PR 1 .EO. XI) PM 1= PMA*W1/WT 

Q9L213 

0C226 

158* 


1 F i PR2.EQ. XI) PM 2= PMA+W2/WT 

- 03C222 

002 3 J 

159* 


1 F l PR 3 »E0 o XI) PM3= PMA*W3/WT 

02Q221 

002 32 

lfcC* 


IFSPR4 .EO. XI) PM4 = PMA*W4/UT 

030243 

002 34 

161* 


PHA = A KA X 1 ( PM A- SR(I),C.) 

090247 

□ 0235 

162 * 


IFlPt.LE.O.IGO TO 1300 

0702 5 6 

00235 

163* 

c 


030256 

□ 02 35 

164* 

c 

IF NO PRIORITY -I REQUESTS EXIST, THEN PROCEED WITH 

C7C756 

30235 

165* 

c 

THE NEXT HIGHER PRIORITY 

0303 5 6 

□ 02 3 V 

166* 


IF ( SR ( I ) .EQ. 3.0) GO TO 1000 

030261 

0024 1 

167* 


IF(RU.LE.PG) GO TO luOG 

07d?63 

00242 

168* 

c 


030263 

0324 2 

169* 

c 

IF THE SUM OF ALL PRlORlTV-I REQUESTS ,ST. POWER 

030263 

03241 

170* 

c 

AVAILABLE, THEN GO AROUND 

OOC263 

00243 

171* 


IF 1 SR 1 1 ) .GT. PL) GO TO 430 

030265 

00243 

172* 

c 


C 73265 

02243 

173* 

c 

THE SUM OF ALL PRIORITY-I REQUESTS .LE. POWER 

Q9G265 

□ 324 5 

174* 

c 

AVAILABLE, SO FULFILL EACH PRIORITY-I REQUEST 

030265 

00245 

175* 


DO 20C K = 1, 4 

0' , Q?74 

U325J 

176* - 


IF !PR(K} .EQ. XI) PIKJ = RIK) 

030274 

00252 

177* 


200 CONTINUE 

030302 

□ 0252 

178* 

c 


000302 

00252 

179* 

c 

UPDATE POWER AVAILABLE 

030302 

00254 

len* 


PL = PL - SR i I ) 

090332 

00254 

lai* 

c 


03D 322 

00255 

16 2* 


GO TO 10C0 

D33305 

00255 

183* 

c 


090305 

00256 

184* 


400 CONTINUE 

090307 

00256 

185* 

c 


030307 

00256 

186* 

c 

THE SUM OF THE PRIORITY-I REQUESTS EXCEEDS THE 

0-30307 

00256 

167* 

c 

POWER AVAILABLE, SO COMPUTE AND ALLOCATE FAIR 

030307 

00256 

168* 

c 

SHARE TO EACH PRIORITY-I PORT 

030307 

□ 02 5 6 

169* 

c 


030307 

00257 

190* 


69C CONTINUE 

030307 

00257 

191* 

c 


090337 

0C257 

192* 

c 

SAVE PL FOR LATER REFERENCE 

D n 039? 

00263 

193* 


POLD = PL 

C9C3CT 

332S0 

194* 

c 


C-9C307 

30260 

195* 

c 

DETERMINE FAIR SHARE UNITS FOR ALL PRIORITV-I 

07C-3Q7 

00260 

196* 

c 

PORTS FOR WHICH NO ALLOCATION HAS BEEN MADE 

07Q3Q7 



00261 

197* 



sum = o*o 

000310 


00262 

198* 



00 700 K = 1, 4 

070 31 4 

ro 

ns 

GQ26S 

199* 



IF (PCK) .NE. 0.0) 60 TO 700 

033 319 

VJl 

ro 

C0267 

200* 



IF (PR IK ) • EQ • XI) SUIZ) r SUII) ♦ U(K) 

000315 


00271 

2bl* 


700 

CONTINUE 

003325 


00273 

202* 



FPU = 1.0 / sum 

000225 


00273 

203* 

C 



CTO 72 5 


00273 

204* 

C 


FIRST, ALLOCATE FAIR SHARE TO PORTS FOR UHICH THE 

000325 


00273 

205* 

c 


FAIR SHARE EXCEEDS THE REQUEST. CONSIDER ONLY PRIORITY-1 

000725 


00273 

206* 

c 


PORTS, AND CONSIDER ONLY PORTS TO UHICH NO ALLOCATION 

000325 


00273 

2l7* 

c 


HAS YET BEEN MADE 

07D325 


UC274 

2ce* 



DO 6 k f> K = 1, 4 

0*0733 


00277 

209* 



IF I P I K 1 . N E • 0.0) GO TO 80C 

070733 


00301 

210* 



IF I PR ( K ) ,NE. XI) GO TO 800 

07C334 


00301 

211* 

c 



070334 


00301 

212* 

c 


COMPUTE FAIR SHARE 

07C734 


00303 

213* 



FRCK ) = ( U ( K ) * FRU) * PL 

070337 


CC303 

214* 

c 



070337 


00303 

215* 

c 


IF FAIR SHARE EXCEEDS REQUEST, THEN FULFILL REQUEST 

070337 


GQ3 04 

216* 



IF (FRCK) .GE. R (K ) ) P(K) - R ( K ) 

070343 


00334 

217* 

c 


- - - AND REDUCE AVAILABLE POWER 

07U343 


J0306 

218* 



IF ( FP 1 K ) .GE. R (K ) ) PL = PL - P|K) 

0 0 L 3 5 1 


30310 

219* 


800 

CONTINUE 

02D362 


00310 

22 C* 

c 



03U362 


00310 

221* 

c 


IF PL .NE. POLO, THEN PL WAS REDUCED DURING THE 

070362 


00310 

22 2* 

c 


PPOCESSING IN THE DO 800 LOOP ABOVE. THIS CHANGES 

O n D362 


00310 

223* 

c 


THE FAIR SHARE COMPUTATION. IT IS THEREFORE 

030352 


00310 

224* 

c 


NECESSARY TO GO BACK THROUGH THE DO 800 LOOP IN 

070762 


00310 

225* 

c 


ORDER TO RECONSIDER ANY PORT WHICH MAY NOW 

U7U362 


00310 

22 6* 

c 


SATISFY THE REQUIREMENT THAT FRCK) .GE. RCK). ONLY 

070362 


00310 

227* 

c 


PRIORITY-I °0R TS FOR WHICH NO ALLOCATION HAS BEEN 

03D362 


Q 03 1 0 

228* 

c 


MADE ARE ELIGIBLE FOR RECONSIDERATION 

070362 


00312 

229* 



IF (PL .NE. POLD) GO TO 60U 

070362 


00312 

230* 

c 



070362 


00312 

231* 

c 


FINALLY, ALLOCATE POWER TO THOSE PORTS REQUESTING 

030362 


003 12 

232* 

c 


MORE Than .THEIR FAIR SHARE. CONSIDER ONLY 

0 73 362 


00312 

233* • 

c 


PRIORITY-1 PORTS FOR WHICH NO ALLOCATION HAS BEEN MADE 

070362 


0G314 

234* 



DO 960 K = 1, 4 

370370 


00317 

235* 



IF |P(K) .NE. u.ul GO TO 9_0 

0 n 0 370 

00 

30321 

236* 



IF CPRIK) • NE • < I ) GO TO 9J0 ' j 

070371 

o 

t A 

00323 

237* 



P(K) - FR(K) 

070774 

l/i 

0 0324 

238* 


900 

CONTINUE 1 

070430 

o 

00326 

239* 



PL = O.C 

07C4D0 

w 

00326 

240* 

c 


1 

030407 

CD 

o 

003 27 

241* 

1000 

CONTINUE 

000404 

1 

03327 

242* 

c 



C7C4C4 

ro 

00327 

243* 

c 



070404 

70 

00327 

244* 

c 


FINALLY, ASSIGN OUTPUTS TO NON-SUBSCRIPTED 

030404 

(X> 

< 

00327 

245* 

c 


FORMAL PARAMETERS 

070404 

• 

00331 

246* 



PI = P(l> 

070404 


00332 

247* 



P2 = P I 2 ) 

030406 


00333 

248* 



P3 = P ( 3 ) 

070410 


0C334 

249* 



P 4 = P (4) 

070412 


00335 

230* 



I F | IMPL .GT. 1> SP = PL 

0:70414 


00337 

251* 



RG=AMIN1 (RO,PM)/EF 

07Q422 


0 0340 

252* 



RETURN 

070«31 


□ 0341 

253* 



END 

070630 
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7.29 PRIORITY INTERRUPT 




This component is used by the storage components to change priority of the 
power requests when minimum or maximum capacity is approached. 


Inputs 

Parameter/Port 
PS 1 

PS 3 

PMX 
I NT 


Description 

Input priority for PS2 output (0 to 4) 

Input priority for PS4 output (default=PSl) 
Maximum priority for PS2 (default = 1) 
Interrupt flag (0,-1, 1) 


Outputs 

Variable/Port 

PS 2 Output priority for charge cycle 

PS 4 Output priority for discharge cycle 

Equations 


PS2 = 

PS1 

if 

I ITT 

=0 

PS2 = 

PMX 

if 

I NT 

>0 

PS2 = 

0 

if 

I NT 

<0 

PS4 = 

PS3 

if 

I NT 

<0 

PS4 = 

0 

if 

I NT 

>0 
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SUBROUTINE PI ENTRY POINT 000045 

r\> 

CTl 

STORAGE USED* CODE(l) 000072* DATA(Q) 000010* BLANK C0HM0NI2) 000000 
COMMON BLOCKS* 

0003 C2MPL 003001 


EXTERNAL REFERENCES (BLOCK. NAME) 
0004 NERR3S 


STORAGE ASSIGNMENT (BLOCK . TYPE. RELATIVE LOCATION. NAME) 



0001 

occois 

IClL 

0003 I 000000 IMPL 0000 000002 INJPS 




00100 

1* 

CPI 



oooooo 


00101 

2* 


SUBROUTINE PI(PS2,PS4,PS1,PS3,PMX,INT) 


000000 


00131 

3* 

C 



OOOGCO 


00131 

4* 

C 

PURPOSE CHANGE PRIORITY OF POUER ALLOCATION TO ST0RA6E 

COMPONENTS 

DOCOOO 


00101 

5* 

C 



000000 


00131 

6* 

c 

WRITTEN BY A. W. WARREN ‘ VERSION 

1. APRIL 14 £977 

000000 


UC131 

7* 

c 



ooonco 


C0101 

8* 

C 

CALL SEQUENCE 


C0QCOO 


CQ101 

9* 

C 

PS2 - OUTPUT PRIORITY (0 TO 4) 


030000 


00101 

10* 

C 

PS4 - OUTPUT PRIORITY (COMPLEMENT TO PS2 ) 


000000 


00101 

11* 

C 

PS1 - INPUT PRIORITV FOR PS2 


030300 


00101 

12* 

C 

PS3 - INPUT PRIORITY FOR PS4 


030300 


00101 

13* 

c 

PMX - MAXIMUM PRIORITY FOR PS2 


030000 


00101 

in* 

c 

int - interrupt flag 


C3C0CO 


00101 

15* 

c 

0= NO interrupt 


030300 


□ 0101 

16* 

c 

1= INCREASE ALLOCATION PRIORITY 


030000 

ro 

00101 

17* 

c 

-1= DECREASE ALLOCATION PRIORITY 


030000 

o 

00101 

18* 

c 



ooonco 

CO 

00103 

19* 


REAL INT 


oounoo 

o 

OCI04 

20* 


COMMON /C IMPL V IMPL 


030000 


00105 

21* 


IF( IMPL.GT.O) GO TO 10 


0003Q0 

00 

o 

u 0 1 0 7 

22* 


IFIPS3.E0. .99999) PS3=PSl 


OOCOD2 

1 

00111 

23* 


lF(PMX.EQ..99999)PMXri. 


030007 

1 

CGI 1 1 

24* 

c 



030337 

TO 

fD 

00113 

25* 


10 PS2=PS1 


000015 

< 

□ 0114 

26* 


PS4-PS3 


000016 

e 

00115 

27* 


IFUNT.GT.O.) PS2SPMX 


000020 


00117 

28* 


1FUNT.LT.0.) PS 2=0 . 


000025 


00221 

29* 


IFUNT.GT.O) PS4 = 0. 


030031 


□ 0223 

33* 


RETURN 


030035 


□ 0124 

31* 


END 
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PU 


HYDRAULIC PUMP 


INPUT POWER (PI) 

MAX. INPUT POWER 
(MP1) 

INPUT EFFICIENCY 
(EF1) 


The hydraulic pump model is based on a constant speed design. The pump is 
assumed to be designed to a nominal operating point and input power. For 
off-design performance the pump efficiency is assumed to be functionally 
related to the square root of the mass flow rate. 

C---' 

Basic Equations 

The output mass flow rate is based on the equations 
M = P1*EFF/{C1*C2*H1) 

EFF = 1 - (1-EFD)*SQRT(MD/M> I* 

where Cl, C2 are conversion constants 
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PU 
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MASS FLOW RATE (M) 

% 

MAX. OUTPUT POWER (MP2) 

I 

OUTPUT EFFICIENCY (EF2) 
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7.30A SOLAR-PHOTOVOLTAIC ARRAY 


PV 


Cell Temperature (TC) 
Solar Insolation (ST) 
Input Terminal Voltage (VT) 



Cell Efficiency (EF1) 

Power Output (P) 

Power Tracker Efficiency (EF2) 


The photovoltaic cell is modeled by the circuit below. Power is delivered 
at terminal voltage V and is dependent on the cell temperature and insola- 
tion. Default for V is the maximum power point. A square array of solar 
cells is assumed with both parallel and series connections. 


RS 



i ' 

wv J 

il C 

) VD 1 

L U 

r 




BASIC SOLAR CELL MODEL 


Basic Equations 

Output current I as a function of terminal voltage V is given by the 
implicit relation 

I = IL + I0*(1-EXP( (V+I*RS)*QBK/(T+273) ) ) (1) 

where 
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PV 


IL = light current (amps) 

10 = diode reverse saturation current (amps) 

T = temperature (°C) 

RS = internal resistance (ohms) 

QBK = device constant (default = electron charge/Boltzmann's constant) 

The light current IL is computed by a bivariate expansion of insolation and 
cell temperature. It has been reported that this model fits observed solar 
cell characteristics within 5% at high temperatures and insolations and 
within less than 1% under more moderate conditions (ref. 2). The reverse 
saturation current 10 is given by 

I0(T) = KD*A0*( (T+273)**3)EXP(~EG0/(T+273) ) (2) 


where 

KD = a device constant 
AO = a material constant 
EGO = band gap at 0°K/Boltzmann' s constant 


Tables Description Units 

EFF Efficiency of maximum power tracker versus 

fractional load (default table provided) 

0P Optimum cell power versus insolation and kw 

temperature (computed table) 

0V Optimum cell voltage versus insolation and volts 

temperature (computed table) 


272B 
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Inputs/Port 

Description 

Units 

VT 

Array terminal voltage (default = maximum 
power voltage) 

volts 

TC 

Cell temperature 

°C 

TL* 

Low temperature value (default = 28) 

°C 

TH* 

High temperature value (default = 120) 

°c 

TR 

Temperature range (default = TH) 

°c 

ST 

Collector solar insolation 

2 

w/m 

SL* 

Low insolation value (default = 1000) 

w/m 2 

SH* 

High insolation value (default - 25000) 

2 

w/m 

SR 

Insolation range (default = SH) 

2 

w/m 

RC 

Concentration ratio (default = 25) 

- 

AA 

Total illuminated cell area 
(default = .00015*NS*NP) 

2 

m 

NS 

Number of cells in series (default = 300) 

- 

NP 

Number of cells in parallel (default = 500) 

- 

11* 

Cell short circuit current at TL,SL 
(default = .06) 

Amps 

12* 

Cell short circuit current at TL,SH 
(default = 1.5) 

Amps 

13* 

Cell short circuit current at TH,SL 
(default = .06) f 

Amps 


These inputs may be ignored if IL1,DS,DT,DST,KD coefficients are sup- 
plied. 
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Inputs/Port 

(cont'd) 

14* 


VI* 


RS 

AO 

EGO 



KD 


CF 

QBK 

RAP 

CC 

CM 


Description 

Cell short circuit current at TH,SH 
(default = 1.56) 

Cell open circuit voltage at TL,SL 
(default = .6) 

Cell internal resistance (default = .055) 

Material constant (default = 1.54E33 for silicon) 

Band-gap at 0°K normalized by Boltzmann's 
constant (default = 1.4E4 for silicon) 

Coefficients in bivariate expansion for the 

light current IL. If not provided, they 

will be computed from the inputs II,..., 14, 


Units 

Amps 

Volts 

Ohms 


°K 


1/°C 

m 2 /w°C 


Device constant, if not provided will be 
computed from II, VI 

Lens radiation transmission coefficient 

|i 

Device constant (default = 1.161E4) °K/V 

Rated power of maximum power point tracker kw 

(default computed) 

9 

O 

Capital cost/year/unit cell area $/m 

Maintenance cost/year $ 


Note: Minimum input parameters to specify PV are cell area AA, number of 

cells in series NS and in parallel NP, concentration ratio RC, and rated 
power RAP. These parameters must be consistent with those for the collec- 
tor model FO or FP. 

\ 

, v 

\ 

* V 

These inputs may be ignored if IL1,DS,DT,DST,KD coefficients are sup- 
plied. 
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Output/Port 

Description 

Units 

V 

Array terminal voltage 

Volts 

P 

Array output power 

kw 

I 

Array output current 

Amps 

EF1 

Solar cell efficiency 

- 

EF2 

Maximum power tracker efficiency 

- 

Statistics 



SP 

Sum of energy delivered 

kwh 


Calculation Sequence 
First Pass 

1) Compute parameter KD (if not input) 

^ KD = 11/ JaO*( (TL+273 )**3 )*EXP ( -EGO/ (TL+273) )*(EXP(QBK*V1/ 

(TL+273) )-EXP(QBK*Il*RS/ (TL+273) ) )] 

2) Compute coefficients IL1,DS,DT,DST (if not input) in the light 
current bivariate expansion in temperature T and insolation S: 

IL = IL1*S*(1+DS*(S-SL)+DT*(T-TL)+DST*(S-SL)*(T-TL) ) (3) 

Define 

FIL(I,T) = I-I(3(T)*(l-EXP(QBK*I*RS/(T+273) ) ) . 

Then 

III = FIL (I1,TL)/SL 

DS = (FIL(I2,TL) - IL1*SH)/(IL1*SH*(SH“SL) ) 

DT = (FIL(I3,TH) - IL1*SL)/(IL1*SL*(TH-TL) ) 

I ' DST = (FIL(I4,TH) - IL1*SH-IL1*SH*DS*(SH-SL) 

- IL1*SH*DT* (TH-TL ) )/ ( IL1*SH* ( SH-SL )* (TH-TL ) ) . 
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;• { 

3) If a terminal voltage VT is not input, calculate the optimal cell 
voltage V=0V(S,T) with S ranging through 10 values equally spaced 
between 0 and SR, and with T ranging through 10 values equally spaced 
between 0 and TR, resulting in a 10 x 10 matrix 0V(S,T). The calcula- 
tion is as follows: Given S and T, the open circuit voltage VOC is 

given by 

VOC = (T+273)*ALOG(1+IL/I0)/QBK, 
where IL and 10 are computed from (2) and (3). 

A binary search is performed in the range from 0 to VOC. For a value V 

in this range, Newton-Raphson iterations are used to solve for the ‘ 

terminal current I satisfying (1). The corresponding power P (in kw) 

is 

P = I*V/1000 . 

The iterative search process to maximize P is given by 

(i) Take the initial interval [VL , VH] to be [0,V0C] . 

(ii) Compute a numerical derivative of P at the midpoint VM of 

[VL,VH] : 

P' = (P(VM+lE-5)-P(VM) )/lE-5 

(iii) If P' > 0, set VL = VM. 

If P' < 0, set VH = VM. 

.<r v 

f 
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(iv) IfVH-VL > 2E-5 and the number of iterations performed is 10, 
go to (ii). Otherwise P is maximized and 
0V(S,T) = VM 
0P(S,T) = P 

The 10 x 10 matrices 0V(S,T) (optimal cell voltage) and 0P(S,T) (maxi- 
mal cell power) are stored for use in subsequent passes. 

Subsequent Passes 

4) Compute insolation S at the cells 

S = ST*RC*CF 

5) If terminal voltage VT is not input, the cell terminal voltage V and 
power P are obtained by interpolation from the arrays 0V(S,T) and 
0P(S,T) . (A diagnostic is printed if S > SR or TC >TR). 

6) If VT is used as an input voltage, then the cell voltage and power are 
determined using 

V = VT/NS 

I = IL(S,TC) + I0(TC)*(1-EXP(QBK*(V+I*RS)/ (TC+273) )) 

P = I*V/1000 

7) Array outputs prior to maximum power tracker: 

V = V*NS ‘ 

P = P*NS*NP 
I = P*1000/V 

EF1 = P*1000/(S*AA) if S>0 

EF2 = 1. 
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8) If the maximum power tracker is used, 

EF2 = EFF( P/RAP) 

P = P*EF2 

REFERENCES FOR PV 

1. J. K. Linn, "Photovoltaic System Analysis Program-SOLCEL," Sandia 
Laboratories Report SAND77-1268, 1977. 

2. L. H. Goldstein and G. R. Case, "PVSS-A Photovoltaic System Simula- 
tion Program," Sandia Laboratories, 1976. 
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PURPOSE THIS COMPONENT COMPUTES THE POWER AND VOLTAGE 
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OUTPUT OF A PHOTO-VOLTAIC CELL ARRAY GIVEN THE 
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HETHOD NEWTON RALPHSON METHOD TS USED TO CALCULATE CELL 

CURRENT AS FUNCTION OF INSOLATION, TEMPERATURE, AND 
TERMINAL VOLTAGE. IF TERMINAL VOLTAGE IS NOT INPUT, 
POWER IS COMPUTED AT OPTIMAL VOLTAGE. THIS IS DONE 
FOR A RANGE OF 10 VALUES OF TFMPERATURE AND 10 
VALUES OF INSOLATION IN THE FIRST PASS. 

AT SUBSEQUENT PASSES, 

INTERPOLATION IS USED. 

CALL SEQUENCE 
TABLES 

EFF -EFFICIENCY OF MAXIMUM POWER TRACKER 
VS FRACTIONAL LOAD (DEFAULT TABLE! 

OP -OPTIMAL POWER, Kb, VS INSOLATION, W/H2, AND 
TEMPERATURES 

OV -OPTIMAL TERMINAL VOLTAGE, V, VS INSOLATION ,W/H2, AND 
TEMPERATURES 

OUTPUTS 

V -ARRAY TERMINAL VOLTAGE, VOLTS 

P -ARRAY OUTPUT POWER, KW 

I -ARRAY OUTPUT CURRENT, AMPS 

EF1 -SOLAR CELL EFFICIENCY 

EF2 -MAXIMUM POWER TRACKER EFFICIENCY 
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SP -SUM OF POWER DELIVERED, KW 
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POWER VOLTAGE! 
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TL -LOW TEMPERATURE VALUE , C , (DE FAULT=Z8 » 
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NS -NUMBER OF CELLS IN SER IES ( OEF AULT=300I 
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IL1,DS,DT,0ST 

-COEFFICTlNTS IN BIVARIATE EXPANSION FOR THE 
LIGHT CURRENT IL. IF NOT PROVIDEO, THEY WILL 
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BE COMPUTED FROM THE INPUTS 11, ..*,14. 
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KD -DEVICE CONSTANT. IF NOT PROVIDED, IT HILL BE 
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8CLTZMANN*S CONSTANT: 1 .16 1E4 1 

030013 

00101 

76* 

c 
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OF TC AND ST 1 

0 9 C 0 1 3 

30191 

BC* 

c 


CC -CAPITAL COST /YEAR/UNIT CELL AREA, 1/K2 
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00107 

37* 
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30 

00111 

38* 


F0=ei*FIN 


00112 

3 5* 


60 TO ICO 


com 

50* 

C 

POSITIVE SATURATION 


00115 

51* 

10 

F0=Cl*C3-»C2*iriN-C3> 

00115 

52* 


60 TO 130 


00115 

53* 

c 

NEGATIVE SATURATION 


ecus 

55* 

20 

F6=C5*C6*CS*(F!N-C4! 

CC116 

55* 


CO TO 100 


00116 

56* 

c 

NEGATIVE UNSATURATEO 


00117 

57* 

30 

F0=C5*f IN 
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0C121 
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7.33A SOLAR ARRAY ORIENTATION 


SO 


-Tracking Power Required (RE1) 
.Array Beam Insolation (SB1) 
■Array Diffuse Insolation (SD1) 
■Array Total Insolation (ST1 ) 


The Solar Orientation model computes flat plate collector insolation for 
five types of solar tracking: 

• Tilted orientation, facing south 

• Tracking about a horizontal EW axis 

• Tracking about a horizontal NS axis 

• Tilted, tracking about a vertical axis 

• Two axis tracking 

Array insolation is the sum of beam and diffuse components. The beam 
component is the product of normal incidence radiation and a geometry- 
dependent incidence factor. The diffuse component is approximated as the 
product of horizontal diffuse insolation times a geometry factor plus 
ground reflectance. 


Day-of-Year (DY) 


Time-of-Dav (TD) 



SO 

Di rect Beam 


Insolation (SB) 


Global Insolation (ST) — — 
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BASIC EQUATIONS 


SO 


ST1 = SB1 + SD1 + SRI ft 

= SB*IF + SD*RD + ST*RR 
SD = ST - SB*SIN(EL) 

RD = .5*(1 + COS(TLT) ) 

RR = ,5*PR*(1 - COS(TLT) ) , 

where 

IF = solar incidence factor (incidence angle cosine) 
TLT = collector tilt angle from horizontal 
PR = ground reflectance 


Inputs/Port 

Description 

Units 

LA 

•fa 

Collector latitude 

Deg 

DY 

Day-of-the-year (1-365) 

- 

TD 

Time-of-day (0-24) 

hr 

M0 

Tracking mode 

- 


1 = fixed orientation and tilt (default) 

2 = horizontal EW axis tracking 

3 = horizontal NS axis tracking 

4 = tilted, vertical axis tracking 

5 = two axis tracking 


TL 

Collector tilt (M0 = 1, 4 inputs) 

Deg 

SB 

Direct normal beam insolation 

2 

w/m 

ST 

Global insolation on a horizontal surface 

2 

w/m 

PR 

Ground reflectance (default = 0.2) 

- 


For TMY stations, see Table 7.7A of the Environmental Data Component ED. 
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Inputs/Port 

(cont'd) 

0 

Description 

Units 

AA 


Collector array area 

m 2 

m 

SBT 


Insolation threshold for tracking 
(default = 100.) 

2 

w/m 

Outputs/Port 

Description 

Units 

SE 


SIN (Solar Elevation Angle)* 

- 

SA 


* 

SIN (Solar Azimuth Angle) 

— 

IF 


COS (Solar Incidence Angle) 

- 

RE 

1 

Tracking power required 

kw 

SB 

1 

Collector beam insolation 

2 

w/m 

SD 

1 

Collector diffuse insolation 

2 

w/m 

SR 

1 

Collector reflected insolation 

2 

w/m 

ST 

1 

Collector total insolation 

2 

w/m 

TLT 


Collector tilt angle 

Deg 



= Solar Orientation Vector 
= Solar Elevation Angle 
= Solar Azimuth Angle 


* Figure 7.33A Solar Orientation Angles 

% < 
^ 'p' 
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CALCULATION SEQUENCE 


SO 


RPD = tt/ 180 
If SB<0 and M0 >1 return 

1) Solar azimuth and elevation 

U = 15*( 12 - TD)*RPD 
<5 = 23.45*SIN(2tt*(284 + DY)/365)*RPD 
LA' = LA*RPD 

SE = SIN 6 *SIN LA‘ + COS 6 *C0S W*C0S LA' 

CE = (1. - SE*SE) 1/2 

TAN(AZ) = COS <J *SIN W/(C0S W*SIN LA'*COS 6 - SIN 6 *COS LA’) 
CA = 1/(1 + TAN 2 (AZ) ) 1/2 
SA = TAN(AZ)*CA 

2) Horizontal diffuse insolation 

SD = ST - SB*SE 


3) Array geometry and tracking power 
RE1 = 0 
If M0 = 1 then 

TLT' = TL*RPD 

IF = SIN TLT'*CE*CA + COS TLT'*SE 
If M0 = 2 then 

IF = yi* " (CE*SA) 2 
TLT' = MIN (COS” 1 (SE/IF) , tt/2) 

RE1 = 3.75 E-4*AA if SB>SBT 

| 
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CALCULATIONS (contd) 

If M0 = 3 then 

IF = Vl- - (CE*CA)2* 

TLT* = MIN (COS” 1 (SE/IF ) , nj 2) 

RE1 = 3.75 E-4*AA 
If M0 = 4 then 

TLT' = TL*RPD 

IF = SIN TLT'*CE + COS TLT'*SE 
RE1 =3.75 E-4*AA 
If M0 = 5 then 

IF = 1 

TLT' = MIN(COS -1 (SE) , n/2) 

RE1 = 5 . E-4*AA 

4) Insolation components 

SBl = SB*IF 

SDl = SD* . 5* ( 1 + COS (TLT ' ) ) 

SRI = ST* . 5*PR* ( 1 - COS (TLT 1 ) ) 
ST1 = SBl + SDl + SRI 

5) Tilt 

TLT = TLT'/RPD 


SO 


if SB >SBT 


if SB > SBT 


if SB > SBT 
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STORAGE USED: CODEdl 000605; DATaIO! 000060s BLANK C0MH0N(2> 000000 
COMMON BLOCKS ? 

CO 03 CIMPL OOOOul 



EXTERNAL REFERENCES (BLOCK « NAME! 

COG* SIN 
C035 COS 
0006 SCRT 
□007 NERR2S 

0010 A cos 

0011 N£ RR3S 

] ':P . i 
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0001 000023 100L 0001 000052 109L 0001 Q0017S 200L 0001 000220 SOIL 0001 0002*6 302L 
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00100 

1* 

CSO 



000000 


00101 

2* 


SUBROUTINE 

SO 1 SE .SA ,IF»RE 1 ,S8 1 1 SOI i SRI t S Tl t TLT, 

000000 


00101 

3* 


1 LA ,0Y ,T0 

,MO,TLfSB,ST,PR,AA,SBT! 

000000 


00101 

** 

C 



030000 

CO 

O 

00101 

5* 

C 

PURPOSE 

THIS COMPONENT COMPUTES FLAT PLATE COLLECTOR 

000300 

CO . 

ZClGl 

6* 

C 


INSOLATION FOR FIVE MODES OF SOLAR TRACKINGS 

ocuooo 


00101 

7* 

C 


TILTED ORIENTATION, FACING SOUTH 

030000 

o 

00101 

8* 

C 


TRACKING ABOUT A HORIZONTAL EH AXIS 

000000 

CO 

0C1CI 

9* 

C 


TRACKING ABOUT A HORIZONTAL NS AXIS 

DOCOOO 

o 

00101 

10* 

C 


TILTED, TRACKING ABOUT THE VERTICAL AXIS 

coopao 

1 

DO 

□ 0101 

11* 

c 


TWO AXIS TRACKING 

(33COOO 


doiui 

12* 

c 



COO 300 

/U 

(T> 

C0101 

13* 

c 

WRITTEN 

BY Y.K.CHAN, 11-6-78, VERSION 1 

aoonao 

< 

□ 0101 

1** 

C 



C3O0QO 


23101 

15* 

C 

METHOD 

ARRAY INSOLATION IS SUM OF BEAM AND DIFFUSE 

DOODOD 


00101 

16* 

C 

. 

COMPONENTS. THE BEAM COMPONENT IS THE PRODUCT OF 

000300 


coiai 

17* 

c 

t 

NORMAL INCIDENCE INSOLATION AND A GEOMETRY DEPENDENT 

030030 


00101 

19* 

c 


INCIDENCE FACTOR. THE DIFFUSE COMPONENT IS 

0 n G?00 


oaici 

19* 

c 


APPROXIMATED AS THE PRODUCT OF HORIZONTAL DIFFUSE 

DPODJO 


0 3101 

20* 

c 


INSOLATION TIMES A GEOMETRY FACTOR PLUS GROUND REFLECTANCE. 

030003 


□ 0101 

21* 

c 



OQCOOO 





a b, 



GC101 

22* 

C 


CALLING SEOUENCE 


oooooo 

CO 

30101 

23* 

C 


OUTPUTS 


000000 

o 

00101 

24* 

c 


SE -SINE. OF SOLAR ELEVATION ANGLE 


oooooo 

CO 

00101 

25* 

c 


SA -SINE OF SOLAR AZIMUTH ANGLE 


000000 


ucini 

26* 

c 


IF -COSINE OF SOLAR INCIDENCE ANGLE 


000000 

o 

— J 

00101 

2 7* 

c 


RE 1 -TRACKING POWER REQUIRE0»KW 


OOUDOO 

00 

00101 

28* 

c 


S61 -COLLECTOR BEAM INS0LAT30N,W/H2 


D3G0J0 

o 

C01D1 

29* 

c 


SOI -COLLECTOR DIFFUSE INSOLATION ,W/M2 


oooooo 

ro 

03101 

30* 

c 


SRI -COLLECTOR REFLECTED INSOLATION, W/M2 


030000 

sO 

00101 

31* 

c 


ST 1 -COLLECTOR TOTAL INSOLA TION, W/HZ 


030900 

<T> 

< 

oo\oi 

32* 

c 


TLT -COLLECTOR TILL ANGLE .DEGREE S 


iiocnoo 

• 

03101 

33* 

c 


INPUTS 

Cr 

oocooo 


uoini 

34* 

c 


LA -COLLECTOR LATITUDE, DEGREES 


0 TC0S3 


00101 

35* 

c 


OY -0 AY OF YEAR! 1-365 ) 


300000 


00101 

36* 

c 


T3 -TIME OF DAY< 0-24 J .HOUR 


D3G n G3 


00101 

37* 

c 


MO -TRACKING MODE 


, 030903 


0010 1 

38* 

c 


INFIXED ORIENTATION AND TILT (DEFAULT) 


030000 


COiCl 

39* 

c 


^HORIZONTAL EW AXIS TRACKING 


oconoo 


00101 

40* 

c 


s=horizontal ns axis tracking 


onanoo 


00101 

41* 

c 


4=tilteo, vertical axis tracking 


003000 


001 01 

42* 

c 


5-TWO AXIS TRACKIGN 


030000 


00101 

43* 

c 


TL -COLLECTOR TILT (MQ=1,4 INPUTS) .DEGREES 


030309 


ocm 

44* 

c 


SB -DIRECT NORKAL PE A H INSOLATION, W/M2 


D90303 


001 J1 

45* 

c 


ST -GLOBAL INSOLATION ON A HORIZONTAL SURFACE, W/MZ 


09G300 


00101 

46* 

c 


PR -GROUND REFLECTANCE S DEFAULT =0.2 ) 


orroooo 


00131 

47* 

c 


AA -COLLECTOR ARRAY AREA, M2 


oocogo 


00101 

48* 

c 


SeT -INSOLATION THRESHOLD FOR TR ACHING, W/M2 


000030 


uom 

49* 

c 


( DEFAULT=1 DO ) 


033030 


03101 

SC* 

c 




OOCOOO 


00103 

51* 



COMMON /CIMPL/IMPL 


OOiiCOD 


00134 

52* 



REAL IF, LA, MO 


030000 


00135 

53* 



2F(IMPL.NE.C)G0 TO 100 


□3CD00 


00107 

54* 



IF(M0.£0. ,99999)M0= ! • 


03DQ01 


QC121 

55* 



IFCPR. E0..99999IPR=*2 


C3C906 


join 

56* 



IF (SbT.EO,. 99999 )SBT=100 


03C013 


03115 

57* 



R PD =3. 1415926/180. 


030923 


3311b 

58* - 


100 

CONTINUE 


030023 


00117 

59* 



IF) (SB.GT.O.) .0R.(M0»LT.2a))60 TO 109 


030 92 3 


0 3121 

6C* 



SA = 0. 


03CC37 


00122 

61* 



1F = C. 


030040 


03123 

62* 



RE1=L . 


030041 


00124 

63* 



SE I =0 . 


030042 


3C125 

64* 



S01=6. 


030343 


0012b 

65* 



SR1=0. 


090344 


33127 

66* 



STl=u. 


030.345 


00130 

67* 



RETURN 


0"0046 


03131 

68* 


109 

continue 


3-10052 

rv) 

00132 

69* 



RE1=0. 


030052 

oo 

00132 

70* 

c 
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w 

rc 

C0132 

71* 

c 


solar AZIMUTH and elevation 


030352 


00132 

72* 

c 
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001 33 

73* 



W=15.*C12.-TD)*9PD 


03035? 


00134 

74* 



ADEL=23.45*SIN( . 017214 2* ( 284*DY) )*RPO 


030357 


0013b 

75* 



PLA=LA*RPD 


090371 


0013b 

76* 



CLAPzCOS ( PL A J 


090974 


00137 

77* 



SADEL = SIK>(ADELI 


03C100 


(J014G 

78* 



CADFL=COSCADEL) 


030134 
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03141 

79* 



SINPICSIN1PL* i 

00142 

80* 



COSW-COSIUI 

OC143 

81* 



SE=SAD£L*SINPLA+CADEL*COSW*CLAP 

OC144 

82* 



CE=SCRTU.-SE*SE> 

OC145 

83* 
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02146 

84* 
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85* 



S A-I • 
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86* 
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87* 
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00153 

88* 
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S A-T A?*CA 

30155 

90* 
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CONTINUE 

00155 

91* 

C 



00155 

92* 

c 


HORIZONTAL DIFFUSE INSOLATION 

J015S 

93* 

c 



00156 

94* 



SD-ST -SB*SE 

03156 

95* 

c 



00156 

96* 

c 


ARRAY GEOMETRY AND TRACKING POWER 

03156 

97* 

c 



0015 7 

98* 



IMO-MO ♦ • 1 

00160 

99* 



60 TOI301, 332, 303, 304, 3051 ,IM0 

03151 

ion* 


301 

TLTP=TL*RPD 

03162 

101* 



IF=SIN ITLTP»*CE*CA*COSITLTPI*SE 

0316 3 

102* 



60 TO 309 

03164 

103* 


302 

IF=SCPT(1.-CE*CE*SA*SA) 

03165 

1l4* 



BIF=AMIN1C1.,SE/IFJ 

0 01 do 

105* 



TLTP = 1 .5708 

0016 7 

106* 



IFCBIF.GT.0.»TL7P=AC0S«B1P» 

03171 

107* 



I F ( Sfc .6T.SBTIRE1=3.7SE-4*AA 

02173 

108* 
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0D174 

1C 9* 
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1F-SQRT J1.-CE*CE*CA*CA ) 

00175 

110* 



6 IF- AM I VI J 1 » , S£/ IF 1 

00176 
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TLTP-1 .5708 

00177 

112* 



IF(BIF.GT.O.JTLTP=ACOS(BIFJ 

0 0231 

113* 
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J 0 2 0 3 
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115* • 
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00205 

life* 



IF -SIN «TLTP>*CE-»COSiTLTPI*SE 

00236 

117* 



IF ( SB. GT. SB TJ RE1-3* 75E“4*A A 

0C21u 

118* 



GO TO 3^-9 

0 0.2 1 i 

119* 
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1 F- 1 . 

0C2 2 2 
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03216 
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CONTINUE. 

00220 
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INSOLATION COMPONENTS 
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12' 8* 
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0 0222 
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0C223 
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00225 

GC226 

00227 


SPIN 


136* TLTrTLTP/RPD 

137* RETURN 

lie* ENO 


000452 

000455 

000604 


GO 


sw 


7.34 SINGLE POLE SWITCH 


V01 



TCI TIME 
SW1=1 


TC2 


TCI 


TIME 

SW1=0 


TC2 


THE SWITCHING OPERATION MAY BE CONTROLLED BY EITHER 
TIME OR THE INPUT PARAMETER SW1. THE TIME DEPENDENCE 
MAY BE ELIMINATED BY SETTING TCI =*10 36 


Inputs 

Parameter/Port 

VAl 

VB1 

SW1 

TCI 

TC2 

Outputs 

Variable/Port 

V01 


Description 
Input to switch 
Input to switch 
Switch control parameter 
Time for first switching (hours) 
Time for second switching (hours) 


Switch output 


Calculation Sequence 
If SW1 = 0 then 
V01 = 


If SW1 = 1 then 
V01 = 


( 

( 


VAl 

TCI < TIME 

< TC2 

VB1 

otherwise 


VB1 

TCI < TIME 

< TC2 

VAl 

otherwise 
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9.0 SOLAR PHOTOVOLTAIC EXAMPLES 


The solar photovoltaic component model's added to the SIMWEST library are 
briefly described and test case results illustrating their use are sum- 
marized in this section. 

Table 9.0-1 summarizes the characteristics of the solar-photovoltaic com- 
ponents. The environmental data component is designed to read Typical 
Meteorological Year (TMY) data tapes containing hourly insolation and 
weather data at 26 U.S. locations. This component can also be used to read 
other hourly data tapes such as the SOLMET tapes by inputing a user speci- 
fied format to the model generation program. The solar orientation or 
tracking component computes the sum of direct beam and global insolation on 
a flat plate array for fixed orientation and four different beam tracking 
options. The flat plate and focusing lens collector components provide 
detailed thermal analyses for determining average solar cell temperature. 
The collector models, and that of the solar array are based on similar 
models developed at Sandia Laboratories for the SOLCEL program. (Refer- 
ence [4]). The array component model is a. simplified model based on 
scaling the characteristics of a single solar cell. Array voltage can 
either be user specified or determined by a maximum power tracker. It 
should be observed that the above components are coded in SI (metric) 
units, whereas most of the SIMWEST components are coded in English units. 
This is generally not a problem since there are at most only a few 
interconnection variables between the solar-photovoltaic generation 
components and other SIMWEST components, and these are easily converted 
using arithmetic components. 

The TMY data tapes are currently the best environmental data sources avail- 
able for simulating typical yearly solar energy system performance. These 
tapes were extracted from SOLMET data tapes containing rehabilitated 
hourly solar and meteorological observation data over a period of many 
years at each observation site. Each Typical Meteorological Year was 
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created by statistical selection of a typical meteorological month for 
each calendar month in the long term data base and catenating the 12 months 
to form a TMY. All of the TMY data files are available for use by a SIMWEST 
user. He thus has access to a high quality environmental data base for 
solar energy simulations and system analyses. 


TABLE 9.0-1 SOLAR-PHOTOVOLTAIC COMPONENTS 


COMPONENT 

SYMBOL 

PURPOSE 

ENVIRONMENTAL 
DATA (TAPE) 

ED 

READ DOE SOLAR INSOLATION AND WEATHER DATA 
TYPICAL METEOROLOGICAL YEAR TAPE 

SOLAR ORIENTATION 
(TRACKING) 

SO 

SOLAR INSOLATION ON TILTED FLAT PLATE 
ARRAY (FIVE OPTIONS) 

FLAT PLATE 
COLLECTOR 

FP 

FLAT PLATE THERMAL MODEL WITH FLUID AND 
PASSIVE COOLING OPTIONS 

FOCUSING LENS 
COLLECTOR 

FO 

FRESNEL LENS THERMAL MODEL WITH FLUID AND 
PASSIVE COOLING OPTIONS 

PHOTOVOLTAIC ARRAY 

PV 

CONVERTS SOLAR INSOLATION TO D.C. ELEC- 
TRICAL POWER. MAXIMUM POWER TRACKER OR 
USER SPECIFIED VOLTAGE 


t 
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9.1 PHOTOVOLTAIC MODEL TEST CASE 


The input data for the photovoltaic model test case is shown in Figure 
9.1-1. The purpose of this model is to obtain characteristic current 
voltage curves for the default solar array parameters. Fortran statements 
are used in the model generation data to let the terminal voltage range 
between 0 and 204 volts for solar insolation values of 5, 20, and 50 suns 
(1 sun = 1000 w/m^). Cell temperature i;|s jspecif ied at 25°C for the first 
simulation and 55°C for the second. Figure 9.1-2 shows the current voltage 
curves and Figure 9.1-3 shows power voltage cross plots at the lower cell 
temperature and for the three solar insolation levels. These curves verify 
the physical characteristics of the solar cell model. It may be noted in 
these figures that current and output power become negative when the speci- 
fied voltage exceeds the array open circuit voltage. Individual cell 
characteristics may be obtained by dividing voltage by 300 (default number 
of cells in series) and by dividing current by 500 (default number of cells 
in parallel). 

9.2 FLAT PLATE COLLECTOR MODEL 

The input data for the flat plate model test case is shown in Figure 9.2-1. 
The purpose of this model is to illustrate water and wind cooling of the 
collector and to test the tracking options of the orientation component SO. 
There are six 1-1/2 day simulation runs. The first run uses water cooling 
(CM0FP=2) , a single glass cover over the front plate and insolation on the 
back. The second run uses passive cooling (CM0FP=0), no plate insolation 
and fins on the back to cool the collector. In the first two runs, the 
collector is tilted and has a fixed, southward facing orientation 
(M0 S0=1). The last four runs are similar to run 2 except different 
tracking options are utilized. 
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MODEL DESCRIPTION PHOTO-VOLTAIC CURRtNT VOLTAGE CURVES 
LOCATION* 1 i TJ 

4— FORTRAN -STATEMENT'S - • 

ST PVS5000 

X F C D V TI ,GT 1 1 1 5) ST PV*20000 

IF-CDV-^IvG^tSJ-ST-PV-SOOOO- — 

VT PV=8,5*TD TI 
LOC ATIONP53 PV 

— END— OF— MODE! " 

PRINT 


a) Model Generation Input Data 


PARAMETER VALUES 
CVCLESbO# TO TI*0 

— DtrINtSsSO - 

TC PV*2*> 

RC PV*J 

— PRINTER-PLOTS #DI SPLA-Y 1 — — • 

V PV#V8#TIME 
I PV#V3#V PV 

P PV# VS# TIME 

TINCP,5#TMAX=72#PRATE*24, PRINT CONTROL*3# INT M0DE*3,0UTRATE* J 

— T i-TL-E ■ PHO TO- VOL T A I C~CE L~L-C UR R EN T-VOL-T A G E- C UR VE 6 • — — 

SIMULATE 

PARAMETER VALUES 

— TC— PV*55 

SIMULATE 


b) Simulation Program Input Data 


Figure 9.1-1 PV Test Case Input Data 
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MODEL DESCRIPTION FLAT PLATE TEST CASE 

LOCATIONS!! TI 

— fcOCA-TIONs35 -fcD INPUTSeTT 

LOCATIONs53 SO iNPUTSaTI ,EDC XlaSB, X2aST) 

LOCATIONS? FP INPUTSaSO#EDCXfl*WD#X3«TA) 

— e*id-of--model — 

PRINT 


a) Model Generation Input Data 


PARAMETER VALUES 

CYCLESb 2 # 01# TO TIa36,TFIFPslo#TFaPPB50#MFMFPs.02#CMOFP«*2#NG FPf!# 

jl — OL-INESs50 — 

HI FPs,01 

CW FPs l # CL Fps2#NT FPslO#CC FPstOOQ#CM FPs i 0 #CPQFPe # 0 1 #LA SOc29,733# 

— TL~ 80b29 ,-733, AA-S0«2 

PRINTER PLOTS# DISPLAY! 

TLTSO# VS# TIME 

— TO— FP #-V8#TTME 

X2 ED#VS#TIM£ 

PI FP#VS#TIME 

— TINC»,5#TMAX«36#PRATEb6#PRINT--LONTROL-p3#-INT-MODE*3#OUTRATEpI 

TITLEepLAT PLATE COLLECTOR TEST CASE 
SIMULATE 

— P A R A M E-T-E-R V A t-OES 

CHOFPsO , Ml FPsl,E9,FIRFPs« 

SIMULATE 

— P A R AME-TE-R -V ALUE S 

MO 8082 
SIMULATE 

-PAR A M ETER V A LHJ ES 

MO SOP 3 

simulate 

— PARAME-TER— VALUE8 

MO SQ=4 
SIMULATE 

— P-A RA ME-TE-H— VA L-UE 3 

MO S0=5 
SIMULATE 


b) Simulation Program Input Data 

Figure 9.2-1 Flat Plate Collector Model Input Data 
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The model schematic produced by the model generation program is shown in 
Figure 9.2-2. The component TI is used to furnish time of day and day of 
year information to SO and to the TMY read component ED. ED supplies direct 
beam and global insolation to SO, and ambient temperature and wind speed 
to the collector component FP. Based on collector orientation, SO supplies 
solar insolation incident to the array, collector tilt angle, and tracking 
power to FP. 

Typical results of the flat plate model runs are shown in Figures 9.2-3 
through 9.2-5. Figure 9.2-3 shows the global horizontal insolation ob- 
tained from ED during the 36 hour simulation period. The data was for mid- 
winter and the daily peak levels are thus low to moderate. The array tilt 
angle daily pattern for horizontal E-W axis tracking is shown in Figure 
9.2-4. At noon the array is oriented normal to the sun's incident rays and 
thus maximizes the insolation gathered during the mid-day peak. The tilt 
angle approaches 90° as the sun approaches the horizon, and remains fixed 
at 90° overnight. Comparison of the solar insolation peaks with the 
various tracking options showed that horizontal E-W axis tracking gave the 
best results of the single axis tracking systems, and was only slightly 
inferior to two-axis beam tracking. Solar cell temperature for this case 
is shown in Figure 9.2-5. The cell temperature is within a few degrees of 
ambient most of the day and rises in mid-day proportional to the solar 
insolation received. The results with water cooling are quite similar. 

9.3 FRESNEL LENS COLLECTOR MODEL AND INCREMENTAL COSTS 

The input data for the Fresnel Lens test case is shown in Figure 9 . 3-1 . 
The purpose of this model is to illustrate a Fresnel Lens collector model 
with thermal fluid loops for collector cooling and for solar heating. 
Three week-long simulations are used to demonstrate incremental cost cal- 
culations for subsystem economic design. A variable speed pump is assumed 
for the collector fluid loop with the flow rate adjusted so that the outlet 
temperature is 5°C greater than the inlet. The collector consists of a 
rectangular grid of 120 Fresnel lenses each of which focuses solar radia- 
tion on a 5 x 5 array of solar cells. Excess thermal energy is conducted to 
a heat sink surface and then dissipated by natural convection, radiation 
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Figure 9.2-2 Flat Plate Model Schematic 
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MODEL DESCRIPTION 

FRESNEL LENS COLLEC 

LOCATION =11 

TI 


L0CATI0N=71 

ED 

I NP UTS = TI 

LOCAT ICN=45 

MA 

INPUT7S = TS(T = FIN) 

FORTRAN STATEMENTS 

"X- jj 

TFOFO 

= -FO MA +5 

• 

LOCAT I0N=33 

FO 

I NP UTS = E D ( X 1=ST, X 3=T 

LOCATI ON=73 

PV 

INPUTS=ED(X1=ST) *FO 

LOCATI 0N=47 

TS 

INPUTS=FO(P,l=P) * TL 

L OC A T I 0N=27 

TL 

INPUTS=TI ,EDCX3=TA) 

L0CATI0N=77 

LO 

INPUTS=PV (P=P,l,P-LO 

LO CAT I 0N = 79 

CM 


END OF MODEL 



PRINT 

a) Model 

Generation 

Input Data 


TOR WITH THERMAL STORAGE AMD 


A, X4=WD) , MA(FO=TFI) 


, 1 ) 


LOAD 


TITLE=FRESNEL LENS COLLECTOR C INCREMENTAL COST COMPUTATION) 
PARAMETER VALUES X^_ 

CYCLES=4.01 ,T0 T 1=0 * CM0F0=2, CW F0=3.75,CL F0=3.9, DLINES=50 

ML FO=120,NT F0=24,MFMF0=Q.5,CC F0=6.,CM FO=50*HI FO=.01»RC FQ=.06 

TS TS=5»DH TS=. 00879, PD TS=12,LE TS=30,NU TS=.01,NC TL=0.2 

Cl MA=.55556»C2 MA=-17.7778, COPFO=0.5 

CC PV=100,CM PV=50,LE TS=30,CR CM=15,LE CM=20 

AA PV=0.6,NS PV=6Q0,NP P V=5 , RA PP V=l. 3 

VE L 0= • 05 , VE TL=.05 

TABLE, HT TS=4 

.0 08 79,. 025491,. 04737 1, ,064072 

90 ,147,147,204 

TABLE, TL0TL=4 

-10, 0,10,25 

4, 2, 1.5,1 

TABLE, TWTTL=4 

0,6, 13 ,24 

• 4, 1,1, .4 

PRINTER PLOTS, DISPLAY1 

RE TL, VS, TIME 

E TS, VS, TIME 

PI FO, VS, TIME 

FMDFO, VS, TIME 

DISPLAY2 

TC FO, VS, TIME 

P PV, VS, TIME 

FO MA, VS, TIME 

INITIAL CONDITIONS=E TS=30 

TINC=.5,TMAX=168,PRATE=12,PRINT C0NTR0L=3 , IN T M0D£=3 , OUTRATE=l 
SIMULATE 

PARAMETER VALUES, TS TS=5.5 

SIMULATE 

PARAMETER VALUES 

TS TS = 5.,NL F 0 = 126, CW F0=3.94»AA PV = 0.63,NS PV=630 
SIMULATE 

b) Simulation Program Input Data 


Figure 9.3-1 Fresnel Lens Model Input Data 
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and heat exchange to the coolant fluid. The collector parameters are 
chosen for a lens concentration ratio of 25 and series connect it>n *of, the 
output from each array. At maximum output the array col lect s ao0ill^'i0kv/ ; ,pf ( * 
solar radiation and produces about 1.7kw of electricalpower. Th^tiser 

'9» 

should be especially careful in specifying the input parameters to the 
collector and array components FO and PV, since inadvertant parameter 
errors can lead to physically inconsistent configurations, e.g., collec- 
tor area smaller than the total lens area. 

The model schematic produced by the model generation program is shown in 
Figure 9.3-2. The collector thermal loop is formed by the connections 
between the collector FO the thermal storage TS and the multiply and add 
component MA. The MA component is used to convert the thermal storage 
outlet temperature from degrees fahrenheit to degrees centigrade. The 
output temperature from MA is supplied as the inlet temperature to FO. The 
total thermal power gathered by the coolant fluid is computed in FO and 
supplied to TS. Similarly, the thermal load fluid loop is represented by a 
power request from the load component TL to TS and by thermal power de- 
livered from TS to TL. The electrical output of the array is computed by PV 
and supplied to a load component LO which monitors the electrical energy 
col lected. 

Results of the first week simulation run are summarized in Figures 9.3-3 
through 9.3-6. The weather was fairly constant during this run and solar 
insolation was fairly strong all week. Figure 9.^-3 shows that with water 
cooling cell temperature was held to less than 70°C at peak insolation. In 
fact, about 60% of the solar energy incident on the array is exchanged to 
the coolant fluid during peak insolation. The electrical output of the 
array is shown in Figure 9.3-4. The fluid flow rate of the pump and 
thermal energy collected exhibit very similar daily patterns. The thermal 
load for this week is shown in Figure 9.3-5. This load is dependent oh 
both time of day and ambient temperature which yields the complex load 
pattern shown. Figure 9.3-6 shows the temperature of the thermal storage 
vessel resulting from the collector and load thermal loops. The daily 

■ ■ \ 
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Fresnel Lens Model Schematic 
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Figure 9.3-3 Solar Cell Temperature for One Week Simulation 
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Figure 9.3-6 Thermal Storage Temperature for One Week Simulation 





cycles are predominant with the periods of strong insolation providing 
sufficient energy to satisfy the load and compensate for thermal losses. 
Average load is fairly well matched to solar generation during the week 
since the temperature remains within a 15° channel and does not have an 
apparent trend away from this range. 

One of the most important measures of performance for a solar energy system 
is the level ized cost of energy, i.e., the life cycle cost to produce one 
unit of usable energy including generation, storage, transmission and con- 
version subsystems. Energy cost may be used to size components and select 
most promising system alternatives, i.e., minimum energy cost is used as a 
selection or optimization principle. Although SIMWEST does not provide 
user optimization capability, optimal sizing of a few key parameters, such 
as the ratio of solar to utility generation and the size of storage rela- 
tive to generation, is possible and may be accomplished quickly using the 
concept of incremental energy cost. The idea is to compute the incremental 
change in level ized energy cost per incremental change in capital cost, for 
the system parameters of interest. Given an initial system configuration 
and M sizing parameters to be selected, optimization proceeds as follows: 

1) Perform M+l back to back simulations to compute the cost and energy 
performance of the baseline configuration and M incremental configu- 
rations from the baseline. 

2) Calculate the incremental energy costs for each parameter variation. 
Then select a new baseline configuration. Since the incremental 
costs are equal at the minimum cost point, increase or decrease the 
sizing parameters so as to equalize the new baseline Incremental 
costs. 

3) Go to 1) and continue adjusting subsystem parameters until either a 
performance limit is reached or until the incremental costs of the 
remaining parameters are equalized. (If two incremental costs are 
unequal, one can always lower the system energy cost by increasing the 
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subsystem with the smallest incremental cost at the expense of the 
other subsystem.) 

This procedure is recommended as more efficient and economical than using a 
series of parametric trade studies for subsystem optimization. 

The process of computing incremental costs is illustrated for the Fresnel 
Lens model. In the first simulation the baseline system performance and 
costs are computed. The second simulation differs from the first in that 
thermal storage capacity has been increased by 10%, and the third simula- 
tion differs from the first in that the solar collector and photovoltaic 
array area have been increased by 5%. Table 9.3-1 summarizes the incre- 
mental cost and simulation results for these runs. Column 1 shows the 
initial capital cost of the baseline system and the incremental capital 
costs for the thermal storage and solar array increases. (These costs are 
meant to be illustrative rather than representative.) Column 2 shows the 
results of a 20 year levelized cost analysis of the three systems, includ- 
ing maintenance and operating costs, e.g., the change in thermal storage 
increases costs by $9.10 per year. Column 3 shows the energy delivered to 
the loads in a year as estimated from the one week simulations. (Note: the 
change in storage capacity lowers the average coolant temperature, thus 
increasing output power.) Column 4 shows the levelized energy costs of the 
baseline system and of the increments in storage and generation. This 
column shows that the levelized energy cost will decrease as thermal stor- 
age or generation are increased, and that thermal storage is undersized 
relative to generation since a fixed $ increase in storage will lower the 
system energy cost more than the same $ increase in array area. Column 5 
shows the % change in levelized energy cost given a 1% increase in capital 
investment. This column contains the same basic information as column 4 
but provides a better quantitative measlure of the economic value of in- 
creased storage capacity. 
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Table 9.3-1 Incremental Cost Calculations 



CC 

LC 

ED 

EC 

NIC 

Basel ine 

7392. 

1272. 

7829. 

16.2 

. ..... 

10% Inc. in Thermal 

61. 

9.10 

110.5 

8.2 

00 

m 

5% Inc. in Solar 

319. 

47.90 

365.0 

13.1 

-.21 


NOMENCLATURE: 


CC = Initial Capital Cost in $ 

LC = Level ized Total Cost/Yr. in $ 

= Capital Cost*Life Cycle*Charge Rate + 
Maintenance Cost + Operating Cost 

ED = Useful Energy Del ivered/Yr. in KWH 
= Electrical Load + Thermal Load + 

Net Change in Thermal Storage 

EC = Level ized Energy Cost in <£/KWH 
= LC*100/ED 

NIC = Normalized Incremental Costs 

= % Change in EC Per % Change in CC 
= (ALC/LC - AED/ED)/(ACC/CC) 
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APPENDIX: UTILITY SUBROUTINES 


This section provides a short description and source code for the utility 
subroutines called by the SIMWEST library components. These routines are 
also available to the user and may be called by FORTRAN statements in the 
user's manual. (See also page 26 of section 2.1.2 on the use of sub- 
routines TBLU1 and TBLU2.) 

• FUNCTION AINR 

AINR computes the current of a photovoltaic cell given light current AIL, 
cell voltage V, and temperature T. Newton-Raphson iterations are used to 
solve the implicit equation (1) for current I: 

I = AIL + BIO (1. - EXP( (V+I*RS)*QBK/(T+273)) ) (1) 

• SUBROUTINE CNVC 

CNVC computes the convection coefficient HC and Reynolds number RE for air 
blown over a flat plate (ref. 1). 

Inputs: T^ - air temperature in °K 

T p = plate temperature in °K 
CL = length of plate in m 
V = velocity of air in m/s 

Equations: 

T M = < T A + V' 2 
VI = 9.0 x 10' 8 *T m -1.115 x 10" 5 

GR = 1.386 x 10 3 - 2.91 *T m 

CO = 7.25 x 10“ 5 *T m + 4.325 x 10" 3 

RE = V*CL/VI 


(mean temp.) 

(viscosity) 

(Grashof's no.) 
(conductivity) 

( 2 ) 
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H FREE = -116*C0*GR*|T A - Tpl * 333 

H WIND = A 597*C0*REE- 5 /CL 

\.032*C0*(RE* 8 - 23000 )/CL 

HC = H FREE + H WIND 

• SUBROUTINE CUBIC 
CUBIC finds the roots of the cubic equation 

x 3 + AAx + BB = 0 (4) 

and selects the real root x with largest value. 

• SUBROUTINE FLUC 

FLUC computes the heat transfer coefficient HF from a collector plate into 
a fluid coolant. The empirical equations used are for water cooling 
(ref. 1). 

Inputs: NT = number of cooling tubes 

DT = diameter of cooling tubes in m 
CW = collector width in m 
COP = conductivity of mounting plate in w/m-K 
THP = mounting plate thickness in m 
FMD = coolant mass flow rate in kg/s 

3 

DEN = coolant density in kg/m 
TF = mean coolant temperature in K 
COC = coolant conductivity in w/m-K 


RE <5 x 10 5 
otherwise 

(3) 


Ns, > 
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Equations: 


NT1 = NT/CW 

HF1 = 12*NT1 2 *C0P*THP (conduction coeff.) 

VI = (21.7*(TF - 256)" 0,8 - .185) x 10" 6 

(fluid viscosity) 

PR = ( .OQ518*TF - 1.25)**(-1.49) (Prandtl no.) 

RE = 4.*FMD/(tt*DT*NT*DEN*VI) (Reynolds no.) (5) 

If RE< 2100, 

HF2 = 4.36*C0C* tt*NT1 

If RE >10000 

HF2 = .023*C0C*RE ,8 *PR‘ 333 * 7r*NTl 
If 2100 < RE < 10000 \ 

X2 = 36 . 5*PR * 33 
D2 = .0029*PR* 33 

A = (4.36-X2)*1.6 x 10” 8 + D2*1.266 x 10" 4 
B = 02 - A*2.xl0 4 
C = X2 + A*10 8 - D2*10 4 
HF2 = (A*RE 2 +B*RE+C)*C0 C*tt*NT1 

HF = (1/HF1 + 1/HF2)” 1 (6) 

• FUNCTION HTGLAS 

HTGLAS computes the top surface heat loss coefficient H t for a collector 
with 1 to 3 glass covers (ref. 2). 

Inputs: N = number of glass covers (1,2,3) 

T a =■ ambient temperature in °K 

= mean cell temperature in °K 

H r = convection coefficient for air blowing over a 
heated flat plate in w/m -k 

e c ,e g = emittance of cell and glass covers 

TLT = collector tilt from horizontal in degrees 
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Equations: 


H t = ( N (V C)/((T C- T A )/{N+f)) °* 33 + VHc)" 1 

+ or( T 2+ T 2)( T c + T A ) / ( A + (2N+f- 1 )/e g -N) (7) 

with a = 5.688 x 10'® w/m 2 -K^ 

C = 365.9 (l.-.00883*TLT+.0001298*TLT 2 ) 
f = (l.-.04*H c +.0005*H 2 )(l.+.091*N) 

A = l/(e c +.05*N(l-e c )) 

• SUBROUTINE IMPLIC 

IMPLIC controls the iteration logic which determines convergence of im- 
plicit variables in the user's system model, and prints convergence diag- 
nostics. (See section 3.6 for a discussion of the iteration and diagnostic 
control logic.) 

t SUBROUTINE RADC 

RADC computes the infrared radiation coefficient HR between two bodies 
with surface temperatures T^ and T 2 . (See section 7.4 of Duffie and 
Beckman, ref. 3.) 

Inputs: t i»^2 = sur ^ ace temperatures in °K 

epe£ = emittances for surfaces corresponding to T^^ 

H r = 5.688 x 10" 8 (T 2 +T 2 ) (T 1 +T 2 )/(eJ 1 + e" 1 -1) (8) 

• FUNCTIONS TBLU1, TBLU2 

TBLU1 and TBLU2 perform one- and two-dimension linear interpolation. A 
binary search is used to locate the nearest grid points for unequally 
spaced data. See section 2.1.2 for subroutine usage within model genera- 
tion FORTRAN statements. 
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• SUBROUTINE UNIF 


UNIF generates uniformly distributed, pseudo-random number sequences in 
the range [0,1]. This routine may be used to obtain random number 
sequences with a specified distribution function. (See for example the 
coding for WD in section 7.47.) 


REFERENCES 
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3. J. A. Duffie and W. A. Beckman, Solar Energy Thermal Processes , Wiley, 
1974. 
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CAINR 

FUNCTION A1NR CAIL»BIO, V »KS *T ) 

NEWTO&i— RALPHSG'J* TG COMPUTE PHQTO-VGLT AIC CELL CURRENT 

F(A)=A-AlL-aiO^Ci.-EXP(QBR*CV+A^RSI/(T+273n ) 
FP(A)=l.+&IG*EXPCGa**( V+A*R$)/eT+273))*QBK*RS/CT+273) 
A=a. 

DO 1 J=l,10 
AN£W=A~F C A) /FP1 A i 
IF ( ( ANEW— A) •Lc» *00001 )C0 TO 2 
1 A=ANEW 
Z AINR-ANEW 
RETURN 
END 
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CIMVC 


SUBRGU7 1 NE CN VC (HC , RE , TP , 7 A , V, CL > 


C 


c 

COMPUTES 

CONVECTION COEFFICIENT HC AND REYNOLDS 

c 

NUMBER RE 

FOR 

AIR BLOWING OVER A FLAT PLAT n. 

c 

CALLED BY 

COMPONENT FO. 

c 

IttFol S 

TA 

-AIR TEMPERATURE, K 

c 


TP 

-PLATE TEMPERATURE,* 

c 


V 

-VELOCITY OF AIR, M/S 

c 


CL 

-LENGTH OF PLATE, M 

L 

1M = ( TA + 1 P ) ♦ . 3 




Vl=9.B-8*7M-I.libh-5 
GR=13S6.-2.91*TM 
CG-7.25t-:^TM+4.32SE-3 
RE=V *CL/V1 

HFR££=.llo+CC^GR*UADSiTA-TP)) ♦♦4.333) ) 
HWIND^^T+CG^SGRT (RE>/CL 

lF(Rc.GT.B.ES) HWIND= .032*CO* 4 R £♦♦< .a ) -23000. )/CL 

HC=rtFREE+MWlND 

RETURN 

END 



•S 
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o r- o o r r> r. n o coo oon r* o r. 


CUBIC 


SUBRuJl 1NL CUBIC ( aA,S8,ANS) 

TER=AA**3/27. 

TERM=B&**2/A • ♦TER 

IF ( AoS ( 1 cRM) «GT. .0001 )0u TO 10 

THREE REAL ROOTS, TWO EQUAL 

************* ***************************************************** 

AB=2.*Cc.RT|-bB/2.) 

Ab&=-AB/2. 

♦*"*-•*■*■♦♦* V***» ***•■■.♦**» 

SELECT POSITIVE ROOT 

****************************************************************** 
AN^=AHAXi ( AB, aBB ) 

RETURN 

10 1M1 cRM.^1 .0. 100 TO 20 

*********** •**-»*« ******************** ********* ******* *********** 
CWE REAL ROOT* TWC CONJUGATE IMAGINARY ROOTS 
♦*♦*♦**♦**************<'***************************♦**************♦ 

STERM=SQRT( TERM) 

AAA-CbRT C— bo/2.+STeRM) 
c-ji=C3RT 1 — Bd/ 2»— ST cRM) 

♦-H 1 **'-*******-* *♦»***♦♦**♦»*****♦* *********** **************** ******* 
SELECT REAL RCOT 

******************** ************************************ ********** 

aNS=AAA«EBB 

RtTURN 

***************************************************** 4 .*^^*****#:*** 

THREE REAL, UNEQUAL ROOTS 

******************************** ******.********** ** **************** 
20 STcK=SGRT*-TER) 

THETA=AC OSC — 6B/2./STER ) 

TL-2.*SQRT(-AA/3.) 

THETA3=THeTA/3. 

Xi= TELCOS ( THETAS ) 

X2=T£*C0S«THETA3+2.09A39) 

X3=TE*C0S4THBTA3+A.188'/9r‘ V \ 

*»*♦* **■*.**.**♦ ********************************************** ** ***** 

SELECT SMALLEST POSITIVE ROOT 

**-* * ******** ********** *************•:****»******♦*!»****** ********** 
— ANS-AMAXI(XlfX2*X3) 

RETURN 

END 
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no on 


cfluc 


c 

o* 

v _ 
c 
C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBkGUT 1 NE FLUC(HF t RE,NT,DT,CW,CCS,THS,FMU,DEM»TF f COC) 

s-GMPUTES HEA'I TRANSFER COEFFICIENT HF TO FLUID 
AND REYNOLDS NUMBER. 

called by component fc 

LrPuTS *T -dumber of cooling tubes 

DT -ulAMtrEK OF COOLING TUBES 
CW -COLLECTOR WIDTH,M 

COS —CONDUCT 1 V 1 1 Y OF MOUNTING PLAT E» W/M— K 

ThS -MOUNTING PLATc. THICKNESS, M 

FMC -COOLANT HASS FLOW RATE, KG/S 

DEo* -COOLANT DENS IT Y f KG/M 3 

TF — MtAN COOLANT TEMPERATURE, K 

CUC -CuOLAMT CONDUCT I Vi TY,W/M-K 

Re. AL NT, Nil 

WrtITE( 6 * 10 B )FMO,CtM,TF,CCC 
lot FORMAT ( lhD, 3 X ,*FLUC INPUTS *, 4 F 1 G. 2 ) 

PK=C .G 0518 *TF- 1 . 2 >)*;M- 1 . 49 ) 

NU=NT/CW 

HF 1 = 12 .*NT 1 *NT 1 *CCS*THS 

VI = ( 2 i. 7 *(TF- 25 t>. )*♦(-. 6 )-. lt >3 )* 1 . E -6 

RE- 4 .*FM 0 /( 3 . 14 16 *DT*N 1 * DEN^V I 5 

IFIRE «,G! = dGe..)GU TO 1 

HF 2 = 4 . 3 b*CGC* 3 . 14 l 6 *Nll 

bG TO i» 

1 IFCRt.GT.iQOOO. )GC 10 2 
X 2 = 36 . 5 *(Pk**< .BE) I 
D 2 =.G 02 y*(PR**( . 33 ) ) 

A= (4.3b— X2)*i.oE-3*D2*1.266E— 4 

b=D 2 -A* 2 .t ,4 

C-X 2 -»-A* 1 . E 8 - 02 * I • E 4 

Hr 2 = (A*Rt.*R£ + &*Rc+C)*CGC* 3 . 14 lb*NT 1 

GO TO 5 

2 CONTINUE 

HF 2 =.G 2 ' 3 *CGC*IRE**C *&)) *(PR**( . 333 ) ) *b . 14 lo*MTl 
5 CONTINUE 

HF= I ./ ( 1 ./HF 1 + 1 ./HF 2 ) 

WRITt( 6 , 109 )HF,RE 

109 FORMAT! 1 HU,^X,*FcUL OUTPUTS *, 2 Fi 0 . 2 ) 

RETURN 

END 
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CHTGLAS 

FUNCTION HI GLAS(NG,TA,1C,HC1,EC,EG,TL1 ) 

TOP HEaT LOSS COEFFICIENT HT FOR GLAS COVERS v CALLED BY FP 


I 

S5. r 


INPUTS 

NG-NUMbER OF GLASS COVERS (1,2,3) 

Ta-AMBIENT TEMPERATURE, K 
TC=MfcAN CELL TEMPERATURE, K 

HCi=COMVECT ION COEFFICIENT FOR AIR BLOWING OVER 
A HEATED FLAT PLATE, W/M2-K 
cC,EG=LMlTTANCE OF CELL AND GLASS COVERS 
TLT=C0LLfcC10R TILT FROM HORIZONTAL IN DEGREES 


REAL NG 

SlGMA=5.6t>a&-8 

C=3toS.9*( 1.-.00SS3*TLT+.D001296* TLT*TLT ) 

F= ( 1 0 A *HC 1+ - OOO ii*rtC I*HC 1 ) * ( 1 . + . 09 l*MG ) 

A'l./i EC+ ■ 05*NG* ( 1 .-EC ) ) 

G=NG*( IC/C)/ 4 4 I TC— TA )/4NG+F ) )**0.33) + l./HCi 

6=SIGMA*(TC*TC+1A*TA)*(TC+TA)/IA+I2.*NG+F-1.)/EG-NG) 

HTGLAS=i./G+tt 

RETURN 

END 
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C IMP L1C 

SUBROUTINE 1MPLAC (CYCLES , DUNES) 

COMM ON/C IMP L/INP L, ICNT /COLDER/ NQX,NOV /COLD/VGlD 
COMMON /CV/ V /CNAMEV/ NAMEV /CTIME/ TIME 
DIMENSION V ( i ) » N aHcV ( 1) * VOLD (II 
C ****** UNIVAC VERSION COut ONLY 

L IFtCYCLES.LE.O. j GO TO 4(i 

C ****** 

IFUMPL.GT.OGG TO 10 
SP=0 

ITERS=C YCLES 

IT £R3= MAXO ( 1 vMINO ( ITERS * 20 ) ) II 

1LINES- AbSt DLINES J 
1TN0 = 0 
IMPL=I 

DC & 1=1, NOV 
3 VOLO(l) = 0- 
10 CONTINUE 

C ****** CDC VERSION CODE ONLY 

lFtCYCLES.GE.l . > GO TO 13 
IMPL=2 

IF( ICMT.GE.ILINES) 1MPL=3 
RETURN 
C ****** 

13 IF ( IMPL.GT • 1 ) GO 10 20 
ITNO= ITNO+1 

IFtllNO.GE.il £RS) 1MPL=2 
ICGN=I 

DO 30 i=i ,NOV 

IF(AbStVtl) l.LT. l.L-fe) GO TL 30 

IF ( AaStVOLDt D-Vtl) ) -GT. 0.03*ABS( V t I ) ) )IC0N = 0 
VOLOtl ) — Vtl) 

30 CONTINUE 

IF t ICON. tu.l i 1MPL=2 

IFiiMPL.EQ.2 .AND. ICNT.GE.ILINES) IMPL=3 
RETURN 
C 
C 

2 0 ITN0=0 

IF( IMPL.GT. 2) GO 10 40 
IFt ICCN.EQ.l ) GO TO 40 
IFtDLlNES.LT.O.) GO TO 40 

ICK = 0 

DO 30 1=1 ,NOV 

IF ( ABSIVtl )) -LT.l.OE-to) GO TO 50 

I F ft ABSi VOLDtl»-Vt I) ) .LI. 0.05*ABS( V( I ) > )&0 TO 50 
X F HICK. EQ . 0 ) WRITE (b,100) TIME 
100 FORMAT ( lHO,10X,3HHME = ,F9.2i 

WRlT'Et 6 , 200 ) NAMEV 1 1 ) , VOLD (I),V(I) 

200 FORM AT 1 1H ,10X,A6,28H NCNCONVER GENCE- OLD VA,LU£= ,F12.3, 
1 13H NEW VALUE=,F12.3 ) 

ICK = 1 

50 CONTINUE 

IFilCK.EG.I) ICNT= ICN T+I 
40 IMPL=4 
RETURN 

tiND 
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CRADC 

SUBROUTINE RADCtHR,! 1, T2,E1®E2) 

\ 

COMPUTES INFRARED RADIATION COEFFICIENT HR ^ 

CALLtO BY COMPONENT FU o 

INPUTS TI.T2 -SURFACE TEMPERATURES • K 

EI.E2 -CORRESPONDING SURFACE EMITTANCES 

HR-i.oBfcE-8*|Tl*TI+T2*12)*<Tl*12)/ll./El+l.,/£2-i.) 

Re TuKN 
tND 




BCS 40180-2 !Rev. 





o o o o <■' n <-« n o n o o r> o 0000 n n r- o r- 


CTBL Ul 


FUNCTION 1 BLUE iXfXT*FT *NDX» NX ) 

PURPOSE ONE DIMENSION LINEAR INTERPOLATION 

CALL SEQUENCE 

X - VALUE OF INDEPENDENT VARIABLE 
XT - ARRAY OF LENGTH ABStNX) CONTAINING X VALUES 
FT - ARRAY OF TABLE VALUES CORRESPONDING TO XT 
fviDX— INDICATOR FOR STEP SPACING 

IF NDX.EQ.O THEN XT CONTAINS EQUAL SPACED DATA 
IF NDX.NE.G THEN XT CONTAINS UNEQUAL SPACED DATA 
NX - ABStNX) IS THE ARRAY LENGTH 

IF NX.LT.O THEN TRUNCATE OUTSIDE TABLE RANGE 
IF NX.GE.O THEN EXTRAPOLATE OUTSIDE TABLE RANGE 

WRITTEN BY A. W. WARREN VERSION I, APRIL 1977 

DIMENSION XT ( i ) • FT II) 

NA=I ABS (NX) 

IFINA.GT.UGC TO b 
T5LU 1-FT i I) 

RETURN 

5 1FINOX.NE.O) GO TO J.00 


EQUI— SPACED TABLE INTERPOLATION 


XO= XT Cl) 

H- XT C 2 I— XT 1 1 8 
XI = (X-XO/H +1. 

2=XI 

IF(l.GT.O) GO TO 10 
TBLUi= FTti) 

IFCNX.GE.ti) T BLU1 = FT I U + 1X1-1. J*IFTt2 J-FTfl) ) 

RETURN 

10 IF(I.LT.NA) GO Tb 20 
TBLUI=FT(NA) 

IF(NX.GE.O) TBLU1= FTCNA) ♦ <XI-NA)*(FT(NA)-FT(NA-1) ) 

RETURN 

2C TBLUI- FT ( 1 ) + f XI-I ) * « FT II +1 ) -FT 1 1) ) 

RETURN 

UNEQUAL SPACED TABLE INTERPOLATION 

100 IFtX.GE.XTt 1) > GO TO 30 
TBLU1=PT ( 1) 

IF(NX.GE.O) TBLUI= FT(1) + t X-XT II) ) * i FTC 2 J-FTIII ) /( XT< 2 ) -XT ( 1) ) 
RETURN 4 

30 IFIX.LT .XT (NA.)) GO TO AO 
TBLUI- FT(NA) 

IF(NX.GE.O) TBLU1=FT (NA) +( X— XT (NA ) )*(FT(NA )— FT ( NA— 1 ) )/(XT(NA) 

1 - XT (NA— 1 ) ) 

RETURN 
AO 1=1 

I&E- NA 

30 1 I=( IGE+I )/2 

IF ( X.LT .XTC II)) Gb TO 60 
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1= II 
GG TQ 70 

bO I»t- II \- 

70 XF1I+1.LT.IGE) GO TO 50 

TBUii= FTUI + <FT(l*!»-Fim>*(X - XT( I) )/(XT< I+ll-Xl Cl) ) 

RcTURAi 

fcUD | /■ 



< % 

:V t; 

**»%. ' 3fT 
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CTbLUZ 

FUNCTION TBLU2(X,Y,XT,YT,FT,IX,IY»NX»NY,MX,MYi 

PURPOSE TWO DIMENSION LINEAR INTERPOLATION 

MEIHGu BINARY SEARCH TO FIND NEAREST GRID POINTS. 

TBLU1 IS USED TO REDUCE THE INTERPOLATION DIMENSION. 

CALL SEQUENCE 

X, Y ~ POINT AT WHICH INTERPOLATION IS DESIRED 
XT » YT~ ARRAYS CONTAINING INDEPENDENT VARIABLE GRID POINTS 
FT - TWO DIEMSNION ARRAY OF VALUES SUCH THAT FTCI,J) 
CORRESPONDS TO XT(I) f YTCj). 

IXflY- INDICATORS FOR GRID SPACING 

IF IX=0 THEN XT CONTAINS EQUAL SPACED VALUES 
IF IX.NE.0 THEN XT CONTAINS UNEQUAL SPACED VALUES 
NX, NY- ABS(NX) ,A6S(NY) ARE THE ARRAY DIMENSIONS FOR XT , YT 
IF NX.LT.O THEN TRUNCATE OUTSIDE XI RANGE 
IF NX.GT.O THEN EXTRAPOLATE OUTSIDE XT RANGE 
LIKEWISE FOR NY AND YT VALUcS. 

MX, MY- DUMMY ARGUMENTS. SET EQUAL TO ABS(NX), ABS(NY). 

WRITTEN BY A.W. WARREN VERSION 1, JUNE 1977 

DIMENSION XT ( 1) ,YT ( I) ,FT ( 1) 

NA = TABS (NX > 

MX = NA 
NB = IABStNY } 

MY = Nb 

1FCNA.GT .1)GC TO 10 
TBLU2 =- TBlUU Y,YT ,FI , 1 Y,NY ) 

RETURN 

10 IF(NB.GT.i)GO TO 20 

TBLU2 = TBLU1(X,XT ,FT,IX*NX ) 

RETURN 

Y OUTSIDE YT TABLE RANGE 

20 IFi Y.GT. YT Cl) )GD 10 100 
E = ( Y— YT (1 ) J/(YT(2)— YT ( 1 1 ) 

FF1 - T BLUI (X,XT,FT(i) , 1X,NX ) 

TBLU2 =FF1 

IF(NY.GT.0)TBLU2 = FF1 + E*( TBLU1 (X,XT , FTC NA+I ) , IX ,NX ) — FF1) 

RETURN 

100 IFi Y.LT. YT(NB))G0 TO 200 

E = ( YT I NB)— Y ) /( YT ( NB )— YT (NB— 1 I ) 

Nftl - NA*( NB— 1 ) +1 
FFI = TBLiil ( X ,XT ,FT(NBI) ,IX,NX) 

TBLU2 = FFI 

IF (N Y.GT .0)1 BLU2 = FF1+ £*(TBLU1 (X,XT,FTI?fBl-NA ) , IX, NX) -FFI> 

RETURN 

YT GRID SEARCH AND INTERPOLATION 

200 IFUY.NE.G)GQ TO 240 

I = (Y — YT( 1 ) ) / ( YT ( 2 )— YT ( I) ) + I. 

GO TO 300 
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IGE = NB 

250 II = t IGE+I )/2 

1F(Y.LT. YT(I1HGG TO 260 
1= II 
GO TO 270 
260 IGE = II 

270 1FCI+1 .LT. I GE ) GO TO 250 

300 E = tY— YT (1 ) )/f YT( I+D— YK I) ) 

11= Ktt*(I-l)+l 

FF1 = TBLU1 (XflXT*FTlIII* IX? NX ) 

T6LU2 = F FI ♦ E*(TBUJUX?XT,FTfIl+NA)?IX v NX) — FF I) 
RETURN 

END 
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C UNI F 


SUBROUTINE UNIF(U.IX) 

COMMON /C1MPL/ IHPL* ICNT » I TEST 
DATA Y/233967./ 

1FC1MPL.&C.G .AND. ITtST.tG.l) IX=43i4fe9 
IF (1X.LQ.1) IX - 431409 
X= AMODt lX*Y,1677?21o.) 

U- X/16777215. 

IX =X 

KE7JRN 

END 
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